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Abstract 

We show that the spectral gap problem is undecidable. Specifically, we 
construct families of translationally-invariant, nearest-neighbour Hamilto¬ 
nians on a 2D square lattice of d- level quantum systems ( d constant), for 
which determining whether the system is gapped or gapless is an undecid¬ 
able problem. This is true even with the promise that each Hamiltonian is 
either gapped or gapless in the strongest sense: it is promised to either have 
continuous spectrum above the ground state in the thermodynamic limit, or 
its spectral gap is lower-bounded by a constant in the thermodynamic limit. 
Moreover, this constant can be taken equal to the local interaction strength of 
the Hamiltonian. 

This implies that it is logically impossible to say in general whether a 
quantum many-body model is gapped or gapless. Our results imply that 
for any consistent, recursive axiomatisation of mathematics, there exist spe¬ 
cific Hamiltonians for which the presence or absence of a spectral gap is 
independent of the axioms. 

These results have a number of important implications for condensed 
matter and many-body quantum theory. 

*tsc25 @ cam.ac.uk 
^dperezga@ucm.es 
*m. wolf@tum.de 


1 



Contents 


I 7 

1 Introduction 7 

2 Preliminaries 9 

2.1 Gapped versus gapless Hamiltonians. 9 

2.2 Undecidability. 10 

3 Main results 12 

3.1 Implications of the results. 13 

4 Extended Overview 14 

4.1 Ground state energy density. 14 

4.2 Wang Tilings . 17 

4.3 QMA constructions. 19 

4.4 Constant local dimension. 21 

4.5 Translational invariance. 23 

4.6 The thermodynamic limit. 25 

4.7 Structure of the paper. 27 

II 29 

5 Unconstrained local Hilbert space dimension 29 

5.1 Undecidability of the spectral gap via tiling . 29 

5.1.1 Ingredient 1: a tiling Hamiltonian. 29 

5.1.2 Ingredient 2: a gapless frustration-free Hamiltonian .... 30 

5.1.3 Reducing tiling to spectral gap. 30 

5.2 Undecidability of low energy properties. 33 

5.2.1 Reducing tile completion to a ground state energy problem 33 

5.2.2 Reduction of the halting problem to arbitrary low energy 

properties. 36 

III 41 

6 Quantum Phase Estimation Turing Machine 41 

6.1 Quantum Turing Machinery. 42 

6.1.1 Turing Machine Programming. 47 


2 



















6.1.2 Reversible Turing Machine Toolbox. 50 

6.2 Quantum phase estimation overview. 57 

6.3 Preparation stage. 58 

6.4 Control- U v stage . 58 

6.4.1 cU k machine . 59 

6.4.2 Iterating over the control qubits. 62 

6.5 Locating the LSB. 64 

6.6 QFT stage. 68 

6.7 Reset Stage. 70 

6.8 Analysis. 71 

7 Encoding QTMs in local Hamiltonians 73 

7.1 Preliminaries . 74 

7.2 Clock Oscillator. 77 

7.3 Initialisation sweep. 79 

7.4 Clock Counter. 80 

7.4.1 Counter TM construction. 80 

7.4.2 Counter TM Hamiltonian. 89 

7.5 Clock Hamiltonian. 95 

7.6 QTM Hamiltonian. 98 

7.6.1 QTM transition rules. 99 

7.6.2 QTM Unitarity. 100 

7.6.3 QTM initialisation sweep. 107 

7.7 Analysis. 107 

8 Quasi-periodic tilings 115 

8.1 Robinson’s tiling . 115 

8.2 Rigidity of the Robinson tiling. 118 

9 Putting it all together 124 

9.1 Undecidability of the g.s. energy density. 124 

9.2 From g.s. energy density to spectral gap. 133 

9.3 Periodic boundary conditions. 139 

10 Acknowledgements 143 


3 





























List of Theorems 

Definition 1 Gapped. 10 

Definition 2 Gapless. 10 

Theorem 3 Algorithmic undecidability of the spectral gap. 12 

Theorem 4 Axiomatic independence of the spectral gap. 12 

Propositions Decidability of g.s. energy density with promise gap ... 16 

Theorem 6 Undecidability of g.s. energy density. 16 

Theorem 7 Reducing tiling to spectral gap. 30 

Corollary 8 Undecidability of spectral gap for unconstrained dimension 32 
Lemma 9 Relating ground state energy to the halting problem .... 35 

Theorem 10 Relating low energy properties to the halting problem ... 36 

Theorem 11 Phase-estimation QTM. 41 

Definition 12 Turing Machine. 42 

Definition 13 Quantum Turing Machine . 43 

Definition 14 Generalised TM and QTM. 43 

Definition 15 Well-formed QTM. 43 

Definition 16 Normal form QTM. 44 

Definition 17 Unidirectional QTM . 44 

Theorem 18 Local well-formedness. 44 

Theorem 19 Quantum local well-formedness. 45 

Definition 20 Proper QTM . 46 

Lemma 21 Subroutine Lemma. 47 

Lemma 22 Dovetailing Lemma. 48 

Lemma 23 Reversal Lemma. 48 

Lemma 24 Copying machine. 50 

Lemma 25 Shift-right machine. 50 

Lemma 26 Equality machine. 51 

Lemma 27 Increment and decrement machines . 53 

Lemma 28 Looping Lemma . 54 

Lemma 29 Binary adder . 56 

Lemma 30 Unary to binary converter. 56 

Lemma 31 Controlled- U QTM. 59 

Definition 32 Computational history state. 73 

Theorem 33 Local Hamiltonian QTM encoding. 73 

Definition 34 Bracketed state. 75 

Definition 35 Legal and illegal states. 76 

Lemma 36 Regexps. 77 

Lemma 37 Evolve-to-illegal . 83 

Definition 38 Well-formed state. 95 

Lemma 39 Well-formed transitions. 95 


4 





































Lemma 40 Evolve-to-illegal . 97 

Definition 41 Standard-form Hamiltonian. 98 

Lemma 42 Invariant subspaces. 108 

Lemma 43 Clairvoyance Lemma. 108 

Lemma 44 Geometrical Lemma. 109 

Proposition 45 Unique ground state. Ill 

Lemma 46 Cross rigidity. 119 

Lemma 47 Robinson rigifity. 119 

Lemma 48 Segment bound. 122 

Lemma 49 Segment rigidity . 123 

Definition 50 Gottesman-Irani Hamiltonian. 124 

Lemma 51 Tiling + quantum layers . 124 

Lemma 52 Robinson + Gottesman-Irani Hamiltonian. 128 

Proposition 53 Diverging g.s. energy. 129 

Corollary 54 Undecidability of g.s. energy with promise. 133 

Theorem 6 Undecidability of g.s. energy density - restated. 133 

Theorem 3 Undecidability of the spectral gap (restated). 133 

Definition 55 Gapless - generalised. 140 

Lemma 56 Periodic to open b.c. 141 


5 





















List of Figures 

Figure 1 Gapped and gapless. 11 

Figure 2 From ground state energy density to spectral gap. 15 

Figure 3 Valid and invalid tilings. 17 

Figure 4 Undecidability of the tiling completion problem. 34 

Figure 5 Counting argument for proof of Theorem 10 . 39 

Figure 6 Quantum phase estimation, control-phase stage. 59 

Figure 7 Quantum phase estimation, inverse QFT stage. 68 

Figure 8 Evolution of the Track 1 clock oscillator.. 78 

Figure 9 The five basic tiles of Robinson’s tiling. 115 

Figure 10 Parity markings to force the position of the crosses. 115 

Figure 11 A possible Robinson tiling of the plane. 116 

Figure 12 Robinson tiling squares. 117 

Figure 13 The five basic tiles of the modified Robinson tiling. 118 

Figure 14 Gottesman-Irani / Robinson “sandwich”. 125 

Figure 15 Boundary tiles for periodic boundary condition construction. 141 


6 
















Part I 


1 Introduction 

The spectral gap is one of the most important properties of a quantum many-body 
system. It plays a pivotal role in condensed matter, mathematical, and fundamental 
physics, and also in quantum computing. 

One of the main goals of condensed matter theory is to understand phase 
transitions and phase diagrams. The behaviour of the spectral gap is intimately 
related to the phase diagram of a quantum many-body system, with quantum 
phase transitions occurring at critical points where the gap vanishes. The low- 
temperature physics of the system are governed by the spectral gap: gapped 
systems exhibit “non-critical” behaviour, with low-energy excitations that behave 
as massive particles, preventing long-range correlations [MH06]; gapless systems 
exhibit “critical” behaviour, with low-energy excitations that behave as massless 
particles, allowing long-range correlations. 

Many seminal results in mathematical physics prove that specific systems are 
gapped or gapless. Examples include the Lieb-Schultz-Mattis result showing that 
the Heisenberg chain for half-integer spins is gapless [LSM61] (extended to higher 
dimensions by Hastings [ las04j), or the proof of a spectral gap for the ID AKLT 
model [Aff+88]. Other results only apply to gapped systems. Important examples 
include exponential decay of correlations [MH06], the ID entanglement area-law 
[Has07], or the very recent proof that the ground state energy density problem in 
ID is in P [LVV13], 

The spectral gap also arises as a key quantity in quantum computing, via 
adiabatic quantum computation [ ^ar+OO]. In this type of computation, proven 
to be exactly as powerful as the usual circuit model for quantum computation 
[Aha+07], the system is initialised in the ground state of an easy Hamiltonian, that 
is a Hamiltonian for which the ground state can be easily constructed. Then the 
interactions are slowly changed, to end up in a target Hamiltonian whose ground 
state encodes the output of the computation. The key property controlling the 
efficiency of this method is the minimal spectral gap in the path of Hamiltonians 
connecting the initial Hamiltonian to the target. A computational problem has an 
efficient quantum algorithm iff there exists such a Hamiltonian path for which the 
minimal spectral gap is lower-bounded by an inverse-polynomial in the system 
size. 

Thus understanding whether a given system is gapped or not is one of the 
fundamental questions in quantum many-body physics. Indeed, there are many 
famous open problems concerning the spectral gap. A paradigmatic example is the 
antiferromagnetic Heisenberg model. Even in ID, the case of integer spin remains 
open as the “Haldane conjecture” [ Ial83], first formulated in 1983. The same ques- 
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tion in the case of 2D non-bipartite lattices, such as the Kagome lattice, was posed 
by Anderson in 1973 [And73], The spectral gap of such systems remains one of the 
main unsolved questions about the long-sought topological spin liquid phase. This 
latter problem has gained a lot of attention in recent years [BallO] due to materials, 
such as herbertsmithite [Han+12], whose interactions are well-approximated by 
the Heisenberg coupling, and which—according to the latest numerical evidence 
[YHW11] —are potential candidates for topological spin liquids. In the related 
setting of quantum field theory, determining if Yang-Mills theory is gapped is one 
of the Millennium Prize Problems [IWOO], and is closely related to one of the most 
important open problems in high-energy physics: explaining the phenomenon of 
quark confinement. 

All of these problems are specific cases of the general spectral gap problem: 
given a quantum many-body Hamiltonian, is the system it describes gapped or 
gapless? Our main result is a proof that the spectral gap problem is undecidable. 
This is much stronger than merely showing that a problem is computationally hard. 
It implies that it is not merely difficult to determine whether a systems is gapped or 
gapless; it is logically impossible to say in general. 

This has two subtly distinct meanings. First, the spectral gap problem is 
undecidable in precisely the same sense as the Halting Problem is undecidable: 
there cannot exist an algorithm, no matter how inefficient, to determine if a general 
Hamiltonian is gapped or gapless. This algorithmic undecidability relates to solving 
the spectral gap for a family of Hamiltonians, as we might have if we want to map 
out the phase diagram of a system as we vary some external parameters. Second, 
the spectral gap problem is in general an undecidable question: there exist simple 
Hamiltonians for which neither the presence nor the absence of a spectral gap is 
provable from the axioms of mathematics. Theorems 3 and 4 in Section 3 make 
both of these meanings precise. 

The idea that some of the most difficult open problems in physics could be 
mathematically proven to be “impossible to solve” is not new. Undecidability of 
other physical quantities has been shown in many-body systems for the much easier 
cases where either the many-body lattice structure or the translational-invariance is 
removed [Gu+09; Llo93; PER89]. Indeed, proving such impossibility theorems is 
highlighted as one of the main open problems in mathematical physics in the list 
published by the International Association of Mathematical Physics in the late 90’s, 
edited by Aizenman [Aiz9£ ]. Our result, then, can be seen as a major contribution 
to this. 
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2 Preliminaries 

2.1 Gapped versus gapless Hamiltonians 

The following sets the framework for our analysis and defines the objects and terms 
used. Let A(L) := {1,..., L} 2 be the set of sites (or vertices) of a square lattice of 
size L e N, which we assume to be at least 2. By & c A(L) x A(L) we denote the 
set of edges of the square lattice, directed such that O', j) e & implies that j lies 
north or east of i. Two cases will be considered: periodic boundary conditions 
where the outer rows and columns are connected along the same direction as well, 
so that A(L) becomes a square lattice on a torus, and open boundary conditions 
where these connections are not made. 

We assign a Hilbert space 70 (!) =* to each site i e A(L) and the tensor 
product (X) ; . €S . 'H <0 to any subset S c A(L). To every neighbouring pair (i, j) e 8, 
we assign a Hermitian operator h°' J} e £>bH il} <S> f H (}) ) describing the interaction 
between the sites. In addition, we may assign an on-site Hamiltonian given by a 
Hermitian matrix e < BdH (k) ) to every site k e A(L). 

Throughout, we consider Hamiltonians that are built up from such nearest- 
neighbour and possibly on-site terms in a translational invariant way. That is, when 
identifying Hilbert spaces, h\ k) = hf for all k,l e A(L) and = h ihJ) if there is 
a v € Z 2 so that (/', /) = (i + v,j + v). The total Hamiltonian 

H ML) := h (iJ) + ^ hf (1) 

(i,j)eS k€A(L) 

can thus be specified by three Hermitian matrices: a dxd matrix h\ and two d 2 x d 2 
matrices /? I()W and h ct)] , which describe the interactions between neighbouring sites 
within any row and column respectively. Hence, it may alternatively be written as 

hML) = Z £ 'fir 11 + E E h '^ l) + E < 2 > 

rows c columns r ieA(L) 

max{||fi row ||, ||/? C0 |||, p! ||} is called the local interaction strength of the Hamiltonian 
and can be normalised to be 1. 

Let spec H AiL} := jd 0 , ...} denote the spectrum, i.e. the set of eigenvalues 
of H ml> listed in increasing order To < dj < .... For clarity we will sometimes 
write the Hamiltonian in question as an argument of the eigenvalues. Aq(H ML} ) 
will be called ground state energy, the corresponding eigenvector ground state 
and a Hamiltonian H A(L) frustration-free if its ground state energy is zero while 
all h {Ui >,hf are positive semi-definite. That is, a ground state of a frustration-free 
Hamiltonian minimises the energy of each interaction term individually. 

The quantity we are interested in is the asymptotic behaviour of the spectral 

gap 

A(H ml> ) := A\(H AtL> ) - A u (H a,l> ) (3) 
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in the thermodynamic limit, that is, when L —» oo. The main result of this paper is, 
loosely speaking, that it is in general undecidable whether or not this gap closes in 
the thermodynamic limit. To be more precise we need the following definitions: 

Definition 1 (Gapped) We say that a family \H A{L] ) of Hamiltonians, as described 
above, characterises a gapped system if there is a constant y > 0 and a system size 
Lq such that for all L > Lq, Aq(H A{L] ) is non-degenerate and A(H a<l> ) > y. In this 
case, we say that the spectral gap is at least y. 

Definition 2 (Gapless) We say that a family \H A(L) ) of Hamiltonians, as described 
above, characterises a gapless system if there is a constant c > 0 such that for all 
s > 0 there is an Lq e IT so that for all L > Lq any point in \Aq(H A{L} ), Aq(H ml> ) + c] 
is within distance sfrom spec H A(L) . 

Note that gapped is not defined as the negation of gapless; there are systems 
that fall into neither class. The reason for choosing such strong definitions is 
to deliberately avoid ambiguous cases (such as systems with degenerate ground 
states). A Hamiltonian that is gapped or gapless according to these definitions 
should be recognised as such throughout the literature. Our constructions will 
allow us to use these strong definitions, because we are able to guarantee that each 
instance falls into one of the two classes. Indeed, we could further strengthen the 
definition of “gapless” without changing our undecidability results or their proofs, 
below, by demanding that c = c(L) grows with L so that lim t ^ TO c(L) = oo. 

2.2 Undecidability 

There are two different albeit related notions of undecidability, which we will 
briefly and informally recall in the following: one, which can be called axiomatic 
undecidability, is more related to mathematical logic whereas the other, which 
might be called algorithmic undecidability, feels more at home in theoretical 
computer science. For a rigorous introduction and a more detailed discussion of 
the notions used below we refer to [Bar77]. 

A decision problem (i.e. a YES/NO assignment to a each member of a set of 
inputs) is algorithmically undecidable if there is no algorithm that terminates and 
provides the correct answer upon every input. The notion of an ‘algorithm’ used 
here is formally defined in terms of Turing machines, but the Church-Turing thesis 
suggests that any other reasonable formalisation will do as well. A simple but 
crucial point to note here is that decision problems with only finitely many possible 
inputs are trivially decidable by case distinction. In particular, single-instance prob¬ 
lems like "is the spin-1 anti-ferromagnetic Heisenberg model gapped?" are trivially 
decidable in the algorithmic sense since either the algorithm "PRINT:’YES’" or 
the algorithm "PRINT:’NO’" will provide the correct assignment. 
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A i(H) 
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^ jA(ff )> 7 

A() 


( a ) ( b ) 

Figure 1: (a) A gapped system has a unique ground state Ao(H ) and a constant 
lower-bound y on the spectral gap A (H) in the thermodynamic limit, (b) A gapless 
system has continuous spectrum above the ground state in the thermodynamic limit. 
See Definitions 1 and 2 for the precise definitions. 



This differs from the second notion of axiomatic undecidability, also called 
axiomatic independence. Within a formal system, which comprises a set of axioms 
together with rules defining formal proofs, a statement is undecidable in the 
axiomatic sense (or, synonymously, independent) if it can neither be proven nor 
disproven from the axioms. 

There is a simple, but for our purpose important, relation between these two 
notions of undecidability: if a decision problem is algorithmically undecidable, 
then for every consistent and recursive formal system in which the problem can be 
stated, there are infinitely many instances of the problem each of which can neither 
be proven nor disproven from the axioms. The argument is standard [Pool2], so we 
only sketch it here informally. Suppose for contradiction that all but finitely many 
instances are axiomatically decidable. Since the theorems in our formal system are 
recursively enumerable, this would imply that there exists an algorithm which runs 
through all theorems, until either a proof or a disproof is found. The exceptions, 
since there are finitely many, could be dealt with by case distinction. 

In the main part of this paper we prove undecidability of problems concerning 
the spectral gap, the ground state energy density and other low temperature proper¬ 
ties of quantum spin systems in the algorithmic sense. Axiomatic independence 
then follows as a corollary in the way just stated. 
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3 Main results 

For each natural number n, we define <p = (pin) as the rational number whose binary 
decimal expansion contains the digits of n in reverse order after the decimal. We 
also fix one particular Universal Turing Machine and call it UTM. 

Theorem 3 (Algorithmic undecidability of the spectral gap) We construct ex¬ 
plicitly a dimension d, d 2 x d 2 matrices A, B, C, D and a rational number J3 so 
that 

(i) . A is Hermitian and with coefficients in 7L + (57L + -^=Z, 

(ii) . B, C have integer coefficients, 

(iii) . D is Hermitian and with coefficients in {0,!,/?}. 

For each natural number n, define: 

hi(n) = a(n)t (a(n) an algebraic number) 
h mw (n) = D (independent of n) 

h col (n ) = A +p(e i7up B + e^B' + e in2 ~ M C + . 


Then: 

(i) . The local interaction strength is < 1. I.e. all terms h\(n), h row (n), h co i(n) have 

operator norm bounded by 1. 

(ii) . If UTM halts on input n, then the associated family of Hamiltonians \H Ml) (n)) 

is gapped in the strong sense of Definition 1 and, moreover, the gap y > 1. 

( iii) . If UTM does not halt on input n, then the associated family of Hamiltonians 

\H ML] (n)\ is gapless in the strong sense of Definition 2. 

Using the classic result of Turing [Tur36] that the Halting Problem for UTM 
on input n is algorithmically undecidable, we can conclude that the spectral gap 
problem is also undecidable. As noted above, this immediately gives the following: 

Theorem 4 (Axiomatic independence of the spectral gap) Let d € N be a suffi¬ 
ciently large constant. For any consistent formal system with a recursive set of 
axioms, there exists a translationally-invariant nearest-neighbour Hamiltonian on 
a 2D lattice with local dimension d and algebraic entries for which neither the 
presence nor the absence of a spectral gap is provable from the axioms. 

It is straightforward (if tedious) to extract an explicit value for d in Theorems 3 
and 4 from the construction described in this paper. 
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3.1 Implications of the results 

These results have a number of implications for condensed matter physics and 
fundamental physics. They imply that one can write down simple models whose 
phase diagrams are uncomputably complicated. The standard approach of trying to 
gain insight into such models by solving numerically for larger and larger lattice 
sizes is doomed to failure. The system could display all the features of a gapless 
model, with the gap of the finite system decreasing monotonically with increasing 
size. Then, at some threshold size, it may suddenly switch to having a large gap. 
Not only can this threshold be arbitrarily large; the threshold itself is uncomputable. 
In general, we can never know whether a large many-body system is approaching 
the asymptotic behaviour of the thermodynamic limit—one more row of atoms may 
completely change its properties. Our findings also imply that a result showing 
robustness of the spectral gap under perturbations, as for the case of frustration-free 
Hamiltonians in [MP13], cannot hold for general gapped systems. 

Phase diagrams with infinitely many phases are known in quantum systems in 
connection with the quantum Hall effect, where fractal diagrams like the Hoftstadter 
butterfly can be obtained [OAOl; Hof76], Since membership in many fractal sets is 
not decidable (when formulated in the framework of real computation; see [BS93]), 
it would be interesting to see whether quantum Hall systems could provide a 
real-world manifestation of our findings. 

Conjectures about the spectral gap, such as the Haldane conjecture, the 2D 
AKLT conjecture, or the Yang-Mills mass gap conjecture, implicitly assume that 
these questions can be answered one way or the other. Our results prove that the 
answer to the general spectral gap question is not determined by the axioms of 
mathematics. Whilst our results are still a long way from proving that any of these 
specific conjectures are axiomatically undecidable, they at least open the door to 
the possibility that these - or similar - questions about physical models may be 
provably unanswerable. 
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4 Extended Overview 

Since the proof of our spectral gap undecidability result Theorem 3 is quite long, 
in this section we give an overview and discussion of the main ideas in the proof. 

The first and most important step in proving undecidability of the spectral gap 
is to prove undecidability of another relevant quantity: the ground state energy 
density. Once we have this, it is relatively easy to “lift” it to undecidability of the 
spectral gap. (More precisely, we give a reduction from the ground state energy 
density problem to the spectral gap problem.) The intuition behind this is illustrated 
in Figure 2; the proof can be found in Section 9.2. 

In fact, undecidability of the ground state energy density is stronger than we 
really need to prove undecidability of the spectral gap. It is sufficient to prove 
undecidability of the ground state energy with constant promise gap ; undecidability 
of the ground state energy density implies that this holds even with a promise gap 
diverging to infinity. 

4.1 Ground state energy density 

As in Section 2.1, consider the square lattice A(L) with edge length L € N but in 
the general case of v e N spatial dimensions, supporting a translationally-invariant 
nearest-neighbour Hamiltonian 

H a(L> := ^ h (iJ) + hf, (4) 

0'J)e£ k€A(L) 

and let c be its local interaction strength. Assume open boundary conditions. The 
ground state energy density is defined as 

E p := lim E p (L), where E p (L) := L~ v A 0 (H A{l) ). (5) 

L —>oo 

The following simple argument shows that this limit is indeed well defined. 
Consider two lattices of different sizes L,L'e N such that L = nL' for some n e N. 
Assume w.l.o.g. that the interaction terms in the Hamiltonian are all positive semi- 
definite. Then H A(L> is, as an operator, lower bounded by the sum of if translates 
of So we have that 

Ao(H A(L) ) > ifA 0 (N A,L '>). ( 6 ) 

On the other hand, we can use a product of n v copies of the ground state of H A(L ) 
in order to obtain an analogous upper bound on the ground state energy of H AiL} of 
the form 

A 0 (H A(l) ) < n v A 0 (H A(L,) ) + 2vn r L ,(v_1) c. (7) 

Dividing both inequalities by L v we are left with 

E P (L') < E P (L) < E P (L') + -jj-. (8) 
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Figure 2: To relate ground state energy density and spectral gap, we need: (i) A 
Hamiltonian H u {tp ) whose ground state energy density is either strictly positive 
or tends to 0 from below in the thermodynamic limit, but determining which is 
undecidable, and (ii) a gapless Hamiltonian H d with ground state energy 0. We 
combine H u ((p ) and H d to form a new local interaction, h((p), in such a way that 
h((p) has (iii) an additional non-degenerate 0-energy eigenstate, and the continuous 
spectrum of H d is shifted immediately above the ground state energy of H u . (a) If 
the ground state energy density of H u (<p) is strictly positive, its ground state energy 
in the thermodynamic limit must diverge to +oo, and h((p) is gapped, (b) Whereas 
if the ground state energy density of H u ((p ) tends to 0 from below, then its ground 
state energy in the thermodynamic limit must be < 0, and h(<p) is gapless. A proof 
can be found in Section 9.2. 
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Hence, an interval of order 0(1 /L') contains both lim inf and lim sup of E p (L) so 
that both must coincide, which proves that lini/ ^ E P (L) is well defined. 

The ground state energy density is an important physical quantity in its own 
right, as well as being our main stepping stone to the spectral gap results. It is 
therefore worth a brief digression to note that the above argument also shows that 
the ground state energy density can be computed to any precision 6 > 0 by exact 
diagonalisation of H ML ) for any L' > 2 vc/6. This immediately implies that the 
ground state energy density problem is decidable if we provide a finite promise 

gap 6. 

Proposition 5 (Decidability of g.s. energy density with promise gap) Let 6 > 0 

be a computable number and consider translationally-invariant nearest-neighbour 
Hamiltonians on a v-dimensional square lattice with open boundary conditions, fi¬ 
nite local Hilbert space dimensions and algebraic matrix entries. Then determining 
whether E p < 0 or E p > 6 is decidable under the promise that E p £ (0, 8). 

Since the real algebraic numbers form a computably ordered field whose 
cardinality is countably infinite, we can think of the input Hamiltonian as being 
encoded as natural number. 

This is in sharp contrast to the following, which will form the crucial step in 
our proof: 

Theorem 6 (Undecidability of g.s. energy density) Let d eN be sufficiently large 
but fixed, and consider translationally-invariant nearest-neighbour Hamiltonians 
on a 2D square lattice with open boundary conditions, local Hilbert space dimen¬ 
sion d, algebraic matrix entries, and local interaction strengths bounded by 1. 
Then determining whether E p = 0 or E p > 0 is an undecidable problem. 

In the next two sections, we will discuss two approaches to proving Theorem 6 
that do not work. Section 4.2 describes a purely classical construction based on 
Wang tilings. This gives a Hamiltonian with the correct spectral properties, but it 
necessarily requires unbounded local Hilbert space dimension d. In Section 4.3, 
we briefly review the Kitaev-Feynman-style local Hamiltonian constructions used 
in recent QMA-hardness results. By a new and careful application of the quan¬ 
tum phase estimation technique, described in Section 4.4, extending ideas from 
Gottesman and Irani [GI09] in Section 4.5, this approach can give a Hamiltonian 
with constant local dimension. But it necessarily fails to have the required spec¬ 
tral properties. Finally, in Section 4.6, we discuss how combining ideas from 
both these approaches allows us to achieve the required spectral properties whilst 
simultaneously keeping the local dimension constant. 
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4.2 Wang Tilings 

The first approach one might consider to proving undecidability of the ground state 
energy is to note the close relationship between tilings and (classical) Hamiltonians, 
recalling Berger’s [Ber66] classic result that the tiling problem is undecidable. 

We will soon see that this approach is too weak to prove Theorem 6. Nonethe¬ 
less, not only is it helpful to understand why this approach breaks down, the much 
more involved construction required to prove our main result will also make use of 
Wang tilings, albeit in a less direct way. 

A unit square whose edges are coloured with colours chosen from a finite set is 
called a Wang tile. A finite set T of Wang tiles is said to tile the plane Z 2 if there 
is an assignment Z 2 —> T so that abutting edges of adjacent tiles have the same 
colour. The result we will use is the fact that there exists no algorithm which, given 
any set of tiles as input, decides whether or not this set can tile the plane —tiling 
is undecidable [Ber66], Here, rotations or reflections of the tiles are not allowed - 
they would make the problem trivial. 1 

K Z ►* X 

III / u „ 

KJd BB 

Figure 3: Examples of valid and invalid tilings (bottom) of a set of four Wang tiles 
(top). 

A tiling problem can easily be represented as a ground state energy problem for 
a classical Hamiltonian (i.e. one that is diagonal in a product basis). The mapping 
is straightforward: with the identification T - {1,..., T] we assign a Hilbert space 
'H il) - C 7 to each site i of a square lattice, and define the local interactions via 

h (uj) ^ |m><m| (i) <g> \n)(n \ (j) , (9) 

(m,n)eC^ 

! If one slightly modifies the rules of the game and requires complementary rather than matching 
colours for abutting edges, then the problem remains undecidable even if rotations and reflections 
are allowed. 
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where the set of constraints C ( ' ,7) c T x T includes all pairs of tiles (m, n) which 
are incompatible when placed on adjacent sites i and j. As is standard in this 
context, we have used (and will be using throughout) Dirac notation. The overall 
Hamiltonian on the lattice A(L) is then 



( 10 ) 


(U)e£ 


Undecidability of the ground state energy of H c with a promise gap of 1 now 
follows immediately from undecidability of tiling, and this gives undecidability of 
the ground state energy density in the case of open boundary conditions. The full 
proof is described in Section 5. 

However, there is a crucial and fundamental limitation to this approach: there 
is no upper-bound on the local dimension of the Hamiltonian. Rather, the local 
dimension grows with the number of tile types. And we cannot impose any bound 
on the latter, or else the tiling problem is restricted to a finite number of cases and 
is trivially decidable by case enumeration. 

On the other hand, this does already allow us to prove a weaker form of 
our main result: undecidability of the spectral gap for families of Hamiltonians 
with no constraint on the local dimension. From the discussion in Section 2.2, 
this immediately implies existence of Hamiltonians for which the spectral gap is 
axiomatically independent, albeit Hamiltonians whose local dimension is arbitrarily 
large and unknown (in fact uncomputable). Furthermore, this approach can be 
extended to prove that any low-energy property that distinguishes a Hamiltonian 
from a gapped system with unique product ground state is undecidable. Section 5 
of the paper is devoted to proving these results. 

Nonetheless, from a physical perspective e.g. of characterising the phase dia¬ 
gram of a system, it is unreasonable to allow the local Hilbert space dimension 
to grow arbitrarily large - or indeed to change at all - as the parameters of the 
Hamiltonian are varied. So we are still a long way from proving our main result. 

Fundamentally, the problem is that the corresponding Hamiltonians are too 
simple. For a problem to be algorithmically undecidable, it must admit a countably 
infinite number of problem instances. If the local Hilbert space dimension is fixed, 
the Hamiltonian is completely specified by a finite number of matrix elements 
defining its local interactions. The only way to encode a countably infinite number 
of problem instances is to exploit the fact that the matrix elements themselves can 
take a countable infinity of values (e.g. arbitrarily precise rational numbers, or even 
computable numbers). Whereas the above tiling approach is only sensitive to the 
pattern of non-zero matrix elements. 

To overcome this, we will need an inherently quantum approach, which is the 
topic of the next section. 
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4.3 QMA constructions 

There is by now a standard approach to proving complexity-theoretic hardness 
results for local Hamiltonian problems. The idea, which dates back to Feyn¬ 
man [Fey8i ] and was significantly developed by Kitaev [KSV02] and others [KKR06; 
OT08; Aha+09; GI09], is to construct a Hamiltonian whose ground state encodes 
the history of a quantum computation in superposition. If we divide the system into 
two registers, a “clock” register and a “computational” register, then the desired 
ground state is 



( 11 ) 


where \if/ t ) denotes the state of the computation after t steps. This superposition 
over the history of a computation is often called a computational history state. 

It is not difficult to construct a Hamiltonian with this as its unique ground state. 
The difficult part is to implement Feynman’s idea using a local Hamiltonian. This 
was first done by Kitaev [KSV02], who showed how to construct such a Hamilto¬ 
nian out of 5-body terms. Kempe, Kitaev, and Regev [KKR06] improved this to 
2-body, Oliveira and Terhal [OT08] to nearest-neighbour two-qubit interactions on 
a 2D square lattice, and Aharonov et al. [Aha+09] to nearest-neighbour two-body 
interactions on a line. 

All of these constructions exploit the fact that the interactions can differ from 
site to site, in order to encode arbitrary computations. Indeed, for translationally- 
invariant nearest-neighbour interactions on a regular lattice, the entire Hamiltonian 
is specified by a finite number of two-body terms (and possibly one single-body 
term), and it might appear that there are not enough parameters available to encode 
arbitrary quantum computations. However, in a remarkable paper, Gottesman and 
Irani [GI09] showed how to construct a translationally-invariant Hamiltonian which 
has as its ground state a computational history state for an arbitrary computation. 
(In fact, the two-body interaction in the Gottesman and Irani [GI09] construction is 
the same fixed interaction for all problem instances. The input is specified by the 
only remaining free parameter: the length of the chain!) 

The aim of all these local Hamiltonian constructions was to prove QMA- 
hardness of the finite-size ground state energy problem for the corresponding class 
of Hamiltonians, by encoding the quantum computation that verifies the witness for 
a QMA problem and adding a local term to the Hamiltonian that gives an additional 
energy penalty to the “no” output. 

An obvious approach to constructing a Hamiltonian with undecidable ground 
state energy is to use one of these local Hamiltonian constructions to encode the 
evolution of a universal (reversible or quantum) Turing Machine, instead of a QMA 
witness verifier, and give an additional energy penalty to the halting state. If we use 
the Gottesman and Irani [GI09] construction, the resulting Hamiltonian will consist 
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of translationally-invariant, nearest-neighbour interactions on a line. Since the 
Halting Problem is undecidable, and the ground state energy depends on whether 
or not the computation halts, the ground state energy of this Hamiltonian would 
seem to be undecidable. As in the tiling approach, this is certainly too weak to 
prove undecidability of the energy density. But one might hope that it is sufficient 
to prove undecidability of the ground state energy. 

However, this Feynman-Kitaev Hamiltonian approach does not even achieve 
the weaker result of the tiling approach. There are now two crucial problems: 

(i) . The Halting Problem is undecidable for the universal Turing Machine on 

arbitrary input (or for an arbitrary Turing Machine running on fixed input). 
As in the tiling approach of the previous section, it is not at all clear how to 
encode this countably infinite family of problems into the constant number 
of matrix elements describing the nearest-neighbour interaction. 

(ii) . The promise gap 5 in all known local Hamiltonian constructions, and in 

particular the translationally-invariant construction of [GI09], scales inverse- 
polynomially in the system size. Thus (assuming the limits exist) 

|lim^co a(L) non - halting 

lim A () (H a<l> ) = \ w . 6 

L ^°° [lirn^eo a(L) + halting (12) 

= lim a(L). 

l _,—>oo 

I.e. the ground state energy in the thermodynamic limit is identical in both 
the halting and non-halting cases. 

These issues are inherent to the spectral gap problem for many-body quantum 
systems, where the question is only meaningful or interesting in the thermodynamic 
limit of Hamiltonians with regular structure (of which translational-invariance is 
the simplest case). Thus they cannot be side-stepped, and overcoming them is the 
main task in proving the result. 

In the following section, we will see that overcoming (i), whilst challenging, 
can be achieved by exploiting the ability to encode quantum computation. Indeed, 
this will essentially be the only quantum part of our construction. 

However, (ii) presents a more serious obstacle to the history state approach. 
There is an inherent trade-off between run-time and promise gap in Kitaev-style 
local Hamiltonian constructions, and the run-time is directly related to the system 
size when the Hamiltonian is constrained to a lattice. But we are necessarily 
working in the thermodynamic limit of arbitrarily large lattice size. We therefore 
need a constant promise gap, independent of the length of the computation. This 
cannot be achieved by any known local Hamiltonian construction, and may well be 


20 



impossible. Without a constant gap between the halting and non-halting cases, the 
ground state energy problem becomes trivially decidable in the thermodynamic 
limit. We discuss how we overcome this obstacle in Section 4.6. 

4.4 Constant local dimension 

To overcome the unbounded local dimension obstacle we faced in Sections 4.2 
and 4.3, we must find a way of encoding the countably infinite family of Halting 
Problem instances into the finite number of matrix elements describing the local 
interactions of a system with fixed local Hilbert space dimension. 

If we encode the evolution of a quantum Turing Machine into the ground state 
of a local Hamiltonian using a Feynman-Kitaev-style construction, as described in 
the previous section, the local dimension will depend on the number of internal 
states and alphabet size of the QTM. Whichever universal Turing Machine we 
choose to encode, that particular TM will have a fixed state space and alphabet size. 
But to encode the Halting Problem, we need a way to feed any desired input to this 
encoded universal TM. It is not difficult to construct a special-purpose classical TM 
which outputs any given string, starting from a fixed input. But, exactly analogous 
to the Wang tiling constructions of Section 4.2, if there is no upper-bound on 
the number of different strings that we must be able to produce, then either the 
number of internal states or the alphabet size of the Turing Machine is necessarily 
unbounded. This is no use to us, as it would again lead to a family of Hamiltonians 
with unbounded local dimension. 

The only way we can hope to generate arbitrarily long strings using constant 
alphabet size and a constant number of internal states is to use a genuinely quantum 
construction. 1 The transition rules of a QTM can have arbitrarily computable 
numbers as coefficients. (In fact, algebraic numbers will suffice for our purposes.) 
So, whereas for given alphabet size and number of internal states there is only a 
finite number of different classical deterministic TMs, there are a countably infinite 
number of different QTMs. We will show how the string we want to produce can 
be encoded in the transition rule coefficients of a QTM, in such a way that the 
QTM writes out this string and then halts deterministically. 

At first sight, this might appear to violate the Busy Beaver bound on the run¬ 
time of a TM [Rad62], or the Holevo bound on the amount of information that 
can be extracted from a finite-dimensional quantum state [ Iol73], or other results 
that limit the amount of information that can be extracted from a finite-size system. 
However, a little more thought reveals there is no contradiction here. 

Indeed, something similar is already possible for classical probabilistic Turing 
Machines. It is a straightforward exercise to construct a classical probabilistic TM 

indeed, this will essentially be the only point in our construction where we exploit the fact 
that the Hamiltonian is quantum. 
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with fixed alphabet and number of internal states which, given access to a coin with 
bias p, outputs the binary expansion of p with high probability, in expected run¬ 
time that is a function of the length of the binary expansion. What is perhaps more 
surprising is that Quantum Turing Machines allow this to be done deterministically. 

The reason this does not violate the Busy Beaver theorem is that, to simulate a 
probabilistic or quantum TM on a deterministic TM, the alphabet and/or internal 
state size must grow with the precision of the entries in the probabilistic or quantum 
transition function. 

Nor is there any contradiction with the Holevo bound. We are not encoding 
the string in a finite-dimensional quantum state, or even in multiple copies of a 
quantum state. We are encoding the string in the unitary transition rules of a QTM, 
which we get to apply as many times as we like on any quantum state we like. 
Applying the transition rules to a fixed quantum state and performing quantum 
state tomography would already allow us to extract the information encoded in the 
transition rules to arbitrary precision. Again, perhaps the only somewhat surprising 
aspect is that, by exploiting the full power of quantum computation, we can recover 
the encoded string exactly, regardless of how long the string is. 

The idea behind our construction is to use the quantum phase estimation 
algorithm [NCOO] (running on a QTM) to extract a phase which we encode in a 
single-qubit unitary, thereby writing out its binary expansion to the tape. However, 
for technical reasons that appear to be insurmountable, it is crucial to our proof 
that the phase estimation be carried out exactly, not merely with high probability. 
Furthermore, the QTM should halt deterministically after a time that depends only 
on the input. Without these properties, the matrix elements of the Hamiltonians we 
construct will not be computable, and Theorem 3 becomes vacuous. 

It is well-known that the quantum phase estimation circuit can output the 
unknown phase <p exactly if the phase has a finite binary decimal expansion (i.e. 
<p = a / 2 b for some a,b e N), given the appropriate set of two-qubit gates [NCOO]. 
The universal QTM construction of Bernstein and Vazirani [BV97] shows that 
any quantum circuit can be implemented on a QTM. But this only implements a 
quantum circuit up to some error, not exactly. Indeed, it is not at all clear whether 
one can carry out exact phase estimation for arbitrary <p, even in principle, given 
that the precision of <p can be arbitrarily high and is unknown to the quantum phase 
estimation circuit. (Indeed, a no-go result by Nishimura and Ozawa [NO02] shows 
that no quantum Turing Machine can exactly simulate quantum phase estimation 
on an arbitrary number of qubits.) 

We therefore cannot appeal to previous results to show that such a QTM exists. 
Instead, in Section 6 we give a detailed and explicit construction of a phase- 
estimation QTM, that satisfies the very strict requirements of our spectral gap 
undecidability proof. It may seem as if we have achieved nothing, since now we 
have the problem of generating the correct input to the phase-estimation QTM! 
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But in fact this is not an issue because, firstly, the input is encoded in unary, and 
generating a unary string is much simpler than generating an arbitrary string; and 
secondly, it does not matter if the upper-bound in the input is larger than the actual 
number of digits. Our phase-estimation QTM will still write out the phase exactly, 
and halt deterministically in this case. (The run-time of the QTM will depend on 
the size of the upper-bound, but run-time is not relevant here.) So we only need to 
generate a long string consisting entirely of the same symbol! This is very easy, 
especially when we come to encode this QTM into a local Hamiltonian. 

However, it is important to note that, if the input upper-bound is wrong, so that 
the number of binary digits in the phase exceeds the number specified in the input, 
then the only guarantee is that the phase-estimation QTM eventually halts; it could 
leave an arbitrary string (or even quantum state) written to its tape. 

In this way, if we supply the phase-estimation QTM with a suitable input and 
feed its output into a universal TM, the input on which the universal TM runs is 
entirely determined by the transition amplitudes of the phase-estimation QTM. 
Assuming the properties of the Gottesman and Irani [GI09] construction carry 
over, the local Hilbert space dimension will be determined by the alphabet size and 
number of internal states of the QTMs, which is a constant independent of the input. 
Since the alphabet size and number of internal states of any particular universal 
TM are clearly constant, the local Hilbert space dimension of the Hamiltonian 
encoding this sequence of Turing Machines will also be constant. 

4.5 Translational invariance 

Gottesman and Irani [GI09] showed how to construct a fixed Hamiltonian on a ID 
chain, that can encode in its ground state the evolution of a QTM for a number 
of time-steps polynomial in the length L of the chain. The input to the QTM in 
their construction is determined by the chain length. They accomplish this by first 
constructing a translationally-invariant clock to keep track of time, which runs for 
a total of L steps. This clock drives a binary counter TM for L steps, leaving the 
binary representation of L written on the tape. The clock is then reset, and switches 
over to driving the QTM. The binary counting TM and the QTM share the same 
tape, so the input to the QTM is the binary representation of L. It is important 
to note that the local Hilbert space dimension in the Gottesman and Irani [GI09] 
construction depends only on the alphabet size and number of internal states of the 
Turing Machines (plus some constant multiplicative overhead for the clock). 

However, in our case we are interested in the thermodynamic limit. The 
Gottesman and Irani [GI09] result per se does not achieve what we need. We 
cannot use the chain length to encode the input to the QTM, as we are only 
concerned with the limit as the length tends to infinity. Instead, we want to encode 
the input to the QTM in the Hamiltonian itself, and carry out the same computation 
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for any chain length. 1 

Given our quantum phase estimation QTM from the previous section, it is clear 
how we should adapt the Gottesman and Irani [GI09] construction to achieve what 
we need. Instead of the binary counter TM, we first run our phase estimation QTM. 
Provided the chain length is sufficiently large that L > \n\, the phase-estimation 
QTM will write the desired string to the tape and then halt. We then switch to 
driving a universal reversible TM which shares the same tape. 

Whilst this approach does ultimately work, there are a number of technical 
issues to overcome. In particular, the length of the computation in the Gottesman 
and Irani [GI09] construction is limited by the maximum number of time-steps that 
the clock can encode, which is <9(poly L). Whereas our phase-estimation QTM will 
require time 0(poly(L)2 L ) on input of length L. There are a number of ways around 
this. Perhaps the simplest - and the one we adopt - is to modify the Gottesman 
and Irani [GI09] clock construction to count in binary instead of unary, so that the 
clock can encode at least D(2 L ) time-steps. However, this substantially complicates 
the analysis. We also have to address a subtle but important unitarity issue in order 
to carry out the required spectral analysis of our Hamiltonian. 

We give a complete, self-contained construction in Section 7, heavily inspired 
by Gottesman and Irani [GI09], together with a rigorous proof of the required 
spectral properties of the resulting Hamiltonian. 2 

This construction, together with the constant size phase-estimation quantum 
Turing Machine from Section 4.4, allows us to construct a family of translationally- 
invariant ID local Hamiltonians with constant local Hilbert space dimension, 
whose ground states encode the evolution of a universal Turing Machine running 
on an input encoded solely into matrix elements of the Hamiltonian. Supplying 
an appropriate unary input to the phase-estimation QTM (see Section 4.4) can be 
done simply by adding local terms to the Hamiltonian that force the all the qubits 
in the initial state of the computational history superposition to be in the |1) state. 

Thus we have succeeded in overcoming the constant local dimension obstacle 
of the QMA constructions discussed above in Section 4.3. However, simply adding 
a local term to this Hamiltonian that gives an additional energy penalty to the 
halting state does not work, for the reasons discussed in Section 4.3: the energy 
difference between the halting and non-halting cases decreases polynomially with 

'in their paper, Gottesman and Irani [GI09] also provide a construction for infinitely long 
chains. However this works by adding terms to the Hamiltonian which effectively break up the chain 
into segments of length L, and the finite chain-length constmction then goes through independently 
for each segment. This is also not what we want, as it means that, despite the infinitely long chain, 
there is still a finite bound L on the space available for the computation. 

2 For this, we prove a general version of the “Clairvoyance Lemma”, which has been used in a 
number of QMA-hardness proofs in the literature but never formally stated or rigorously proven. 
This may be of independent interest to Hamiltonian complexity theorists. 
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the system size. So all dependence of the ground state energy on the outcome of 
the computation still vanishes in the thermodynamic limit. 

4.6 The thermodynamic limit 

The more challenging obstacle of the thermodynamic limit still remains. 

We first return to tiling problems. However, instead of using these blindly, as 
in the Wang tiling approach described in Section 4.2 where only undecidability of 
tiling was used, we prove and then exploit very particular properties of an aperiodic 
tiling due to Robinson [Rob71], These will allow us, using the ideas discussed in 
Section 4.4 and Section 4.5, to encode in a quantum local Hamiltonian the parallel 
execution of the same universal Turing Machine running on the same chosen input, 
but running on tapes of all possible finite lengths and for every possible finite 
run-time. 

The idea of encoding the parallel evolution of the same Turing Machine on 
tapes of all possible finite lengths, instead of encoding the evolution of a single 
Turing Machine on an infinite tape, dates back to Berger’s [Ber66] original proof 
of undecidability of tiling. Robinson [Rob71] also used the same idea in his 
simplification of Berger’s proof. However, our reason for exploiting this idea 
is somewhat different. In our case, it is the crucial ingredient that allows us to 
decouple the energy dependence of the computational history ground state - a 
purely quantum property - from the overall system size. 

The Robinson tile set can tile the infinite plane, but only aperiodically. The 
aperiodic pattern generated by a Robinson tiling is shown in Figure 11 . We 
have already seen in Section 4.2 that tiling problems can easily be turned into 
translationally-invariant classical Hamiltonians. So we can readily turn this into 
a classical tiling Hamiltonian whose ground state has the same quasi-periodic 
structure as the Robinson tiling. The idea is to add another “layer” on top of 
this tiling Hamiltonian, and use this second “quantum” layer to place copies of 
the Hamiltonian from the previous section along the borders of all the squares in 
Figure 11 . (Of course, this has to be done by adding additional translationally- 
invariant local terms to the Hamiltonian that effectively restrict where the ID 
Hamiltonian acts, not by literally restricting the ID Hamiltonian to the square 
borders, which would break translational-invariance.) 

In this way, we construct a 2D translationally-invariant local Hamiltonian 
whose ground state contains an encoding of the evolution of the universal Turing 
Machine along each edge of the squares in the Robinson tiling. The effective 
tape length of this Turing Machine is limited only by the size of the border it 
“runs” on. But the Robinson tiling contains borders of all finite sizes. 1 So this 
Hamiltonian encodes Turing Machines with all possible tape lengths. All of these 

'Actually, all powers of 4, but this makes no difference to the argument. 
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Turing Machines are running on the same input, encoded in matrix elements of the 
translationally-invariant local interaction. 

If the universal Turing Machine eventually halts on this input, then for all 
borders above a certain size, the effective tape will be sufficiently long for the 
machine to halt before it runs out of tape space. If we add a translationally-invariant 
local term that penalises the halting state, then the ground state will pick up an 
additional energy from the history states encoding Turing Machines that halt. This 
energy still decreases with the size of the system it acts on. But, crucially, this 
size is now the size of the border it is “running” on, not the overall system size. 
We have decoupled the ground state energy from the overall system size; it now 
depends only on the space required for the universal TM to halt. 

If the universal Turing Machine never halts on this input, then it will not pick 
up any additional energy, providing the tape is sufficiently long for the phase- 
estimation QTM to operate correctly. However, if the effective tape length (border 
size) is too small to contain a valid upper-bound on the number of digits in the phase, 
then the quantum state left on the tape by the phase-estimation QTM is arbitrary 
(see Section 4.4). In this case, we cannot assume that the universal Turing Machine 
runs forever, since its input may be corrupted. Thus, even in the non-halting case, 
the ground state will pick up some additional energy from borders that are too 
small. But, crucially, we know the maximum size of these borders: it is determined 
by the number of binary digits in the phase encoded in the matrix elements of the 
Hamiltonian, which we chose. The energy contribution from the small borders is 
therefore a computable number, e.g. brute-force diagonalisation of the Hamiltonian 
on the small borders. Indeed, since this is an eigenvalue computation, this quantity 
is even an algebraic number. The number of borders of any given size grows 
quadratically in the lattice size, so the total energy contribution from the small 
borders can be removed by subtracting an appropriately weighted 1 term from the 
Hamiltonian. 

We have managed to construct a family of Hamiltonians whose ground state 
energies depend on the solution of the corresponding Halting Problem. However, 
this is still not sufficient to prove Theorem 6 and hence prove Theorem 3. The 
difference in energy between the halting and non-halting cases (the “promise 
gap”) depends inverse-polynomially on the space required for the universal Turing 
Machine to halt. Thus, not only does this fail to provide a uniform bound on the 
promise gap as required, this promise gap is uncomputable. 

In fact, this is a limitation of the above analysis, not the true situation. Instead 
of being uncomputably small, as the above argument suggests, the true promise 
gap for this Hamiltonian is infinite! To show this requires a more careful analysis 
of the low-lying eigenvalues (low-energy excitations). It is easy to see that the 
eigenstate consisting of a valid tiling together with computational history states 
along the borders has energy that diverges with the lattice size in the halting case; 
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once the lattice is large enough, the number of borders that are sufficiently large 
for the encoded Turing Machine to halt grows quadratically, and each of them 
contributes a small but non-zero energy. 

The difficulty is that, since its energy diverges, this eigenstate cannot be the 
ground state in the halting case; we can lower the energy by “breaking” some of the 
Turing Machines, by introducing defects in the tiling layer. In Section 8, we prove 
strong rigidity properties of the Robinson tiling, which show that the tiling pattern 
is robust against defects: any defect (non-matching adjacent tiles) in the tiling only 
affects the pattern of borders in a finite ball around the defect. Thus destroying 
n borders requires O(n) defects, each defect contributes 0(1) energy, and Turing 
Machines on intact borders contribute 0(1) energy. Thus, no matter how many 
defects we introduce, the energy will still grow quadratically with the lattice size. 
The promise gap therefore diverges quadratically in the thermodynamic limit. Thus 
we in fact obtain a stronger result than the uniform bound on the promise gap 
required for our main result; since the ground state energy diverges quadratically, 
and our system is on a 2D lattice, this in fact proves undecidability of the ground 
state energy density , as claimed in Theorem 6. 

Making this argument rigorous is somewhat involved, and is the content of 
Section 9.1. Section 9.2 then combines this result with any frustration-free gapless 
quantum Hamiltonian (of which there are many well-known examples) to complete 
the proof of our main result. Theorem 3: undecidability of the spectral gap for 
2D translationally-invariant Hamiltonians with fixed finite local Hilbert space 
dimension and local interactions specified by algebraic matrix elements, even when 
the Hamiltonians are guaranteed to be either gapped with a spectral gap uniformly 
lower-bounded by the local interaction strength (Definition 1), or gapless with 
continuous spectrum above the ground state (Definition 2). 

4.7 Structure of the paper 

We warm up in Section 5 by giving a complete proof of the weaker but substantially 
easier undecidability result for the spectral gap of families of Hamiltonians with un¬ 
bounded local Hilbert space dimension. We also extend this to prove undecidability 
of essentially all low-energy properties of quantum many-body Hamiltonians. 

The bulk of the paper is concerned with proving our main result: undecidability 
of the spectral gap for Hamiltonians with fixed local dimension. The proof is 
broken down into several stages, each described in a separate section. Each of these 
sections is self-contained, in the sense that only the main theorem of the section is 
needed in subsequent sections. 1 Proofs in any individual section can therefore be 
skipped on a first reading, as long as the main result of the section is provisionally 

Except in the case of Section 8, which collects a number of new results on rigidity of the 
Robinson tiling, needed in subsequent sections. 
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accepted. 

In Section 6, we give an explicit construction of an exact quantum phase estima¬ 
tion Turing machine, which will be crucial to achieving constant local dimension. 
In Section 7 we describe a modification of the Gottesman and Irani [GI09] encoding 
of a quantum Turing machine in a ID, local, translationally-invariant Hamiltonian. 
Section 8 contains the analysis of the rigidity of the Robinson quasi-periodic tiling, 
which will allow us to achieve the required spectral properties. 

With this machinery in place, Section 9 proves our main result. This is divided 
into two steps. As discussed in Section 4, the first step is to prove undecidability of 
the ground state energy density problem, which we do in Section 9.1, making use of 
all the results from the previous sections. The second step, in Section 9.2, reduces 
the ground state energy density problem to the spectral gap problem, completing 
the proof of our main result. Section 9.3 extends these results to periodic boundary 
conditions. 
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Part II 


5 Unconstrained local Hilbert space dimension 

In this part we will describe two approaches that exploit known undecidability 
results for tiling and completion problems. Based on these, we can derive undecid¬ 
ability results for the spectral gap and for other low energy properties of transla¬ 
tional invariant nearest-neighbour Hamiltonians on a square lattice. In contrast to 
the approach in the main part of this paper, however, the families of Hamiltonians 
that we construct here, defined by their local interactions j(/? row (n), h co \(n))}„ el ], 
have local Hilbert space dimension d n that differs for different elements n of the 
family, in such a way that {d n } n€ n is unbounded. 

5.1 Undecidability of the spectral gap via tiling 

As in the more sophisticated constructions that will appear later, the idea is to 
reduce an undecidable ground state energy problem to the spectral gap problem. If 
we do not constrain the local Hilbert space dimension, then this reduction can be 
chosen such that it directly exploits the undecidability of a tiling problem. To this 
end, we need two ingredients: 

5.1.1 Ingredient 1: a tiling Hamiltonian 

We saw in Section 4.2 that a tiling problem can easily be represented as a ground 
state energy problem for the local Hamiltonian 


h (iJ, 

:= ^ \m)(m\ (l> <g> |n)(nP , 

(13) 
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where the local basis \m) label tile types, and the set of constraints C <KJ> c 7<" x 7C 
includes all pairs of tiles (m, n) which are incompatible when placed on adjacent 
sites i and j. 

By construction we have h { '' n > 0, H^ (L) is translational invariant and its 
spectrum is contained in No- If there exists a tiling of the plane, then for open 
boundary conditions VL : 0 e spec H^. Similarly, in the case of periodic 
boundary conditions, the existence of a periodic tiling implies that 0 € spec H^ (L> 
holds for an unbounded sequence of L’s. On the other hand, if there is no tiling, 
then there is an L 0 such that 


VL > L 0 : spec H^ (L) > 1. 


(15) 
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Also note that this implies that in the case of open boundary conditions the ground 
state energy density is E p > 1 /(Lq + l) 2 if no tiling exists, compared to E p = 0 if 
there is a tiling of the plane. 

5.1.2 Ingredient 2: a gapless frustration-free Hamiltonian 

As a second ingredient we will use that there are frustration-free two-body Hamil¬ 
tonians 

K L> ■= Z <16) 

(U)e£ 

such that for L —> oo we get spec H q i,J —> R+ in the sense that the spectrum of 
the finite system approaches a dense subset of IR+ [FG+H ]. We will denote the 
corresponding single site Hilbert space by - C D and we can in fact choose 
D = 2 for instance by assigning to each row of the lattice the XY-model with 
transversal field taken at a gapless point with product ground state. Note that in 
this case the entries of the matrices h q are rational. 

5.1.3 Reducing tiling to spectral gap 

In order to fruitfully merge the ingredients, we assign a Hilbert space 'H <0 := 'H^ © 
‘Hf <g> 'Hq’ ^ C 1 © C A ' <g> C D to each site i e A(L). A corresponding orthonormal 
set of basis vectors will be denoted by |0) e C 1 and \k, a ) := | k) © \a) € D , 

respectively. The Hamiltonian is then defined in terms of the two-body interactions 

H (iJ) := |0)<0| ( ° © ig + 1® © |0)<0| 0) (17) 

+ ^ |m)(m|® <g> <S> \n){n^ © (18) 

D 

+ Yj 1 c ) ® \ a )(fl\ {i) ® ® I^Xrl 0) <a, 6\hf j) \/3,y). (19) 

a,j3,y,5=l 

Here l c , t q and t cq denote the identity operators on E( c ,E( q and *H C ® Tf q , re¬ 
spectively. As before, we define the Hamiltonian assigned to a square A(L) as 

H ml> := H<i,J> - ( 20 ) 

(*J)eS 


Theorem 7 (Reducing tiling to spectral gap) Consider any set % of Wang tiles, 
and the corresponding family of two-body Hamiltonians \H MI) \ L on square lattices 
A(L) either with open or with periodic boundary conditions. 

(i). The family of Hamiltonians is frustration-free. 
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(ii). If'K tiles the plane Z 2 , then in the case of open boundary conditions, we 
have that spec H A{L) —> R + for L —> oo, i.e., there is no gap and an excitation 
spectrum that becomes dense in R+. Similarly, if'K tiles any torus, then in the 
case of periodic boundary conditions spec H A(L,) —> R +for a subsequence 
Li -» oo. 

(Hi). If'K does not tile the plane, then there is an Lq such that for all L > Lq H AiL) 
has a unique ground state and a spectral gap of size at least one, i.e. 

spec// A(i) \{0} > 1. (21) 


Proof Frustration-freeness is evident since H >LJ) > 0 and H AiL) |0,..., 0) = 0. In 
order to arrive at the other assertions, we decompose the Hamiltonian as H A(L) =: 
H 0 + H c + H q , where the three terms on the right are defined by taking the sum 
over edges separately for the expressions in (17), (18) and (19), respectively. Let 
us now assign a signature cr e {0,..., K\ L to every state of our computational 
product basis, so that |0) (,) translates to cr, = 0 and \k, a) >n to cr, = k, irrespective of 
a. By collecting computational basis states with the same signature we can then 
decompose the Hilbert space as 




ieA 


( 22 ) 


The reason behind this is that the Hamiltonian is block diagonal w.r.t. this decom¬ 
position, i.e. it can be written as 

H a(L> = (0 H a , and thus spec H A(L) = M spec H a . (23) 

(T cr 


In order to identify the spectra coming from different signatures we will distinguish 
three cases: 


(i) . cr = (0,..., 0): this yields an eigenvalue 0 as already mentioned. 

(ii) . cr ^ (0,..., 0) but cr, = 0 for some i: for any unit vector \f (T ) with such a 

signature we have (t/vlT/olf/v) > 1- We will denote the part of the spectrum 
which corresponds to these signatures by 5 c R. The only property we will 
use is that S > 1. 

(iii) . cr € {1,..., K} l2 \ In the subspace spanned by all states whose signature does 

not contain 0 we have that H 0 = 0 and H c + H q = H AiL> <g> 1 + 1 ® H AiL> . 
Consequently, the spectrum stemming from this subspace equals spec H AiL> + 
spec H A(L) . 
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Putting things together, we obtain 

spec H A(L) = {0} U j spec H AiL) + spec H A<L, \ u S , (24) 

so that the claimed equivalence between the impossibility of tiling with '7C and 
the existence of a spectral gap of size at least 1 follows from the properties of the 
chosen ingredients. Uniqueness of the ground state in cases where no tiling exists 
follows from the above argument when considering the spectra as multisets rather 
than sets. □ 

From the undecidability of the tiling problem we now obtain: 

Corollary 8 (Undecidability of spectral gap for unconstrained dimension) 

Let e > 0. Consider all families \H A(L} ) of frustration-free translational invariant 
nearest neighbour Hamiltonians on 2D square lattices with open boundary condi¬ 
tions that are described by a pair of Hermitian matrices (h row , h co i) x with rational 
entries and operator norm smaller than 1 + 6. There is no algorithm that upon 
input of(h mw , h co i) decides whether these describe a gapped or a gapless system, 
even under the promise that one of them is true and that in the gapped case the 
gap is at least of size 1. 

Here, the norm bound follows from the observation that < 1 + \\hf-’ || where 

hq can be rescaled to arbitrary small norm. 

For the case of periodic boundary conditions we obtain the same statement 
if we replace our strong definition of ‘gapless’ by the weaker requirement that 
3c > OVe > OVL 0 3L > L 0 so that every point in [A () (H AiL> ), Ao(H a<L} ) + c] is within 
distance e from spec H AiL) . The reason for this is, that if the period of the torus does 
not match the required one, then a gap can be generated that, however, disappears 
again if we enlarge the size of the system. 

We emphasise that, so far, there is no constraint on the local Hilbert space 
dimension. Since every instance in the above construction has a finite local Hilbert 
space we can, however, at least conclude axiomatic undecidability for finite local 
Hilbert spaces. That is, for every consistent formal system with a recursive set of 
axioms, there exists an instance of a Hamiltonian with properties as in the corollary 
and finite dimensional local Hilbert space such that neither the presence nor the 
absence of a spectral gap can be proven from the axioms. However, the local 
dimension, whilst finite, is uncomputable. In the main part of the paper, we will 
prove the much stronger Theorem 3, which shows that the spectral gap problem 
remains undecidable even for fixed local dimension. This implies axiomatic 
independence for Hamiltonians where the finite local dimension is known. 

'Here we make no use of an additional on-site term, so h\ = 0. 


32 



In the above construction, gapped cases are related to the impossibility of tiling. 
Since the latter always admits a proof, the axiomatically undecidable cases coming 
from this construction have to be gapless. The following section will provide a 
construction where this asymmetry is reversed. 

5.2 Undecidability of low energy properties 

5.2.1 Reducing tile completion to a ground state energy problem 

We now modify the above construction in order to show that not only the spectral 
gap but many other low energy properties are undecidable as well. The idea, 
illustrated in the case of the spectral gap, is to invert the relation between gap 
and existence of tiling. In the previous construction, the existence of a gap was 
associated with the impossibility of a tiling. Now we want to associate it with the 
existence of a tiling. The drawback is that we loose the frustration free property in 
that case.One of the advantages is that we do not have to rely on the undecidability 
of tiling (which is a rather non-trivial result [Ber66; Rob71]) but can exploit 
the simpler undecidability of the completion problem, already proven by Wang 
[Wan61]. In the completion problem, one fixes one tile at the left-bottom comer 
and asks whether the first quadrant can be tiled with a given set of tiles. 

Undecidability of the completion problem is relatively simple, and explained in 
full detail in Robinson [Rob71, Section 4], by reduction from the halting problem 
for Turing Machines. Any Turing machine is specified by a finite number of states 
Q = [qo,q\,. ..}, a finite alphabet for the tape S = {.s 0 , s\, ...}, where q 0 is the 
initial state and s 0 the blank, and quintuplets of the forms 

(q, s, s',D,q') € Q x S x S x {left, right} x Q 

meaning that if the head of the machine reads 5 and the machine is in state q, then 
it will write s' into the current cell on the tape, transit to state q' and move in the 
direction D. Starting in q, s, there must be at most one valid quintuplet. If there is 
none, then the machine halts. 1 

Following Robinson [Rob71], we define the tiles in Figure 4, where colours are 
identified with labelled arrows, in such a way that matching colours are indicated 
by an arrow head abutting a tail with the same label. We denote the set of tiles by 
7C, which contains at most K := \ ( K\ < 2 + (3|(9| + 1)151 elements. 

We will assume that the Turing machines is constructed such that, when running 
on the initially blank tape, its head never moves to the left of the starting cell. Note 
that this can be assumed without loss of generality, since every Turing machine 
can be simulated by one with this property by standard techniques, which amount 
to ‘folding’ the tape and thereby doubling the tape alphabet. 

'Note that this definition is slightly different from the one used later in Section 6. 
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Figure 4: Tiles for encoding a Turing machine with half-infinite tape into a com¬ 
pletion problem. The tile in the bottom-left corner is fixed and, together with 
the second tile in the third row, starts the Turing machines running on the blank 
half-infinite tape. The tiles in the second row are the ’action tiles’ that are able to 
change the state of the Turing machine and correspond to a left or right moving 
head. For each pair ( q , s) there is at most one action tile. 


If the tile that appears in the bottom-left corner in Figure 4, which we denote 
by is fixed in the bottom-left corner of an L x L lattice, then any valid tiling 
corresponds to the evolution of the Turing Machine on the blank tape up to time L, 
where time goes up in the vertical direction and each row of tiles corresponds to the 
total configuration (including the tape) of the Turing machine in two consecutive 
instants of time (represented by the bottom and upper edges of that row of tiles). 
Hence, 

(1) . There exists a valid tile completion for all L x L squares if and only if the 

Turing Machine, when running on an initially blank tape, does not halt. 

(2) . If such a tiling exists, then it is unique. 

Now, the completion problem can be represented as a ground state energy 
problem in the following straightforward way. Consider a square lattice A(L) with 
open boundary conditions and corresponding edge set S. 

We will assign an energy penalty or energy bonus Waj/m, n) & Z to an edge 
(z, j) e & on which a pair of tiles (m, n) c 7C x 7C is placed. The choice of the 
weights is summarised in the following table. Here □ stands for any tile different 
from E and the numbers in brackets define the weights in cases where the colours 
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(i.e., labels and arrows) match. 



Right tile 


B □ 

Left tile p^ 

4 2 (-1) 

□ 

4 2 (0) 



Tile above 


B □ 

Tile below 

4 2 (-1) 

Lll 

4 2 (0) 


(25) 


Assigning a Hilbert space * - C K to each site i € A(L) whose computational 
basis states are labelled by tiles we define a Hamiltonian H„ (U := Tjuj^s with 

^ := ^ a>vj)(m,n) \m)(m\ m ® |n><n| 0) . (26) 

m,ne9C 

By construction, the Hamiltonian is diagonal in computational basis, its eigen¬ 
vectors correspond to tile configurations A(L) —» and, since all weights are 
integers, its spectrum lies in Z. In order to determine the ground state energy note 
that the tile B is the only one that can lead to negative weights. Assume first 
that a tile B is placed on a site different from the bottom-left corner. Then this 
tile will have at least one neighbour below with an arrow pointing upwards or one 
neighbour to the left with an arrow pointing to the right. In either case there will 
be an energy penalty +4 so that the weights involving this tile will sum up to at 
least 2. If, however, B is placed in the bottom-left corner, then no such energy 
penalty is picked up and the tile can contribute -2 to the energy. If the tiling is 
then completed following the evolution of the encoded Turing machine on a blank 
tape, then no additional penalty will occur on A(L) if until the L’th time step either 
the head of the Turing machines moves left of the origin or the Turing machine 
does not halt. Hence, in both cases Ao(H^) = -2 and the ground state is unique 
and given by a product state. 

If the Turing machine halts, then the completion problem has no solution 
beyond some L 0 and every configuration with B at its bottom-left corner will get 
an additional penalty of at least 2 leading to non-negative energy. Since there is 
always a configuration that achieves zero energy, e.g., one using only the type of 
tiles on the top-left in Figure 4, we obtain: 

Lemma 9 (Relating ground state energy to the halting problem) Consider any 
Turing machine with state set Q and tape alphabet S running on an initially blank 
tape such that its head never moves left of the starting cell. There is a family of 
translational invariant nearest neighbour Hamiltonians [H^ L) } L (specified above) 
on the 2D square lattice with open boundary conditions and local Hilbert space 
dimension at most |S|(3|<2| + 1) + 2 such that 
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(i) . If until time L the Turing machines did not halt, then Ao(H^ L) ) = -2. In this 

case the corresponding ground state is a unique product state and the only 
eigenstate with negative energy. 

(ii) . If the Turing machine halts, then 3Lo (given by the halting time) such that 

VL > L 0 : A 0 (Hu {L) ) = 0. 

Since every Turing machine can be simulated by one with a half-infinite tape, 
which automatically satisfies the constraint that its head never moves to the left 
of the origin, the above construction leads to an undecidable ground state energy 
problem, which we will exploit in the following. 

5.2.2 Reduction of the halting problem to arbitrary low energy properties 

Consider translational invariant Hamiltonians on square lattices with open boundary 
conditions. Let H AiL) describe such a family of Hamiltonians on L x L lattices. Our 
aim is to prove the undecidability of essentially any low energy property displayed 
by this Hamiltonian in the thermodynamic limit that distinguishes it from a gapped 
system with unique product ground state, e.g. the existence of topological order, or 
non-vanishing correlation functions. This is implied by the following theorem: 

Theorem 10 (Relating low energy properties to the halting problem) Let TM 

be a Turing machine with state set Q and alphabet S running on an initially blank 
tape. Consider the class of translational invariant Hamiltonians on square lattices 
with open boundary conditions. Let H A{L) with nearest-neighbour interaction 
hq’^ e &CHq } <g> Llq 1 ) with Li^q — C q and on-site term h^ e describe such 

a Hamiltonian and assume that there is a c e R such that To- cL 2 e 
for all lattice sizes L. Then there exists a family of Hamiltonians H A>U (specified 
in the proof) with the following properties: 

(i) . H A(L> is translational invariant on the square lattice with open boundary 

conditions and local Hilbert spaces LL' 1 - C d of dimension d < 2 (q + 1) + 

I5K3IQI + 1). 

(ii) . The interactions of H A(L) are nearest-neighbour (possibly with on-site terms), 

computable and independent of L. 

(iii) . If up to time L the TM did not halt, then H A(L> has a non-degenerate product 

ground state and a spectral gap A(H ML> ) > |. 

(iv) . IfTM halts, then 3L 0 (given by the halting time) such that VL > L 0 the 

following identity between multisets holds in the interval [0, ^): 

spec H A(L) - A {) (H A(L) ) = spec H ML) - A 0 (H A{L) ). (27) 
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In this case there exist isometries V® : H { q —»77®, V := (^). V® such that 
for this part of the spectrum any eigenstate \<p q ) of H q (t) is mapped to the 
corresponding eigenstate \<f>) of H K(L) via V \4> q ) = | <p). 

Some remarks before the proof. Since, on the Turing machine level, (iii) and 
(iv) cannot be told apart by any effective algorithm, the same has to hold for any 
property that distinguishes a gapped system with product ground state from the low 
energy sector of some frustration free Hamiltonian. With small modifications in 
the proof, the nearest neighbour assumption for h q could be replaced by any fixed 
local interaction geometry. In this case, H A(L) will of course inherit this interaction 
geometry. 

Proof W.l.o.g. c = 0 since we can always compensate for it by adding a multiple 
of the identity to the on-site part of the Hamiltonian. 

For later use we will embed the given frustration free Hamiltonian in a larger 
systems such that its ground state energy gets an offset -1. To this end, define a 
translational invariant nearest neighbour Hamiltonian on an auxiliary system with 
local Hilbert space H ( f - C 2 via 

j]:= Yj »7 (, ' /) . ■= -^I0)<0| ( ° <8> |l)<lp + 21® 0 |0><0| O) . 

(/J)e£ 

Note that Xfrp = -1 and A (rj) = 1. Now introduce Hq := H ( f 0 H q ’ and set 
h { Q ] := T] (, I > 0 1® 7) + 1® 7) 0 ltf ) acting on HqH^q (up to reordering of Hilbert 
spaces). Then Hq L) := Yjpj) e£ Ip 0 h ( ^ has the property that in the 

interval ( - oo, we have that 

sp ecH^ {L) - 1 = spec Hq L> 

holds as an identity between multisets (i.e., including multiplicities) and the corre¬ 
sponding eigenstates are related via a local isometric embedding. 

We will now combine this with the Hamiltonian from the previous subsection 
by enlarging the local Hilbert space to H <1> := Hf © Hq and defining 

h (iJ) ;= h <U) + h (iJ) + 71 (0 ® i0\ ( 28 ) 

where Irt 1 ^ and h ( f’ ] act nontrivially on Hq 0 H ( q and H ( ,f 0 H ( J ) respectively and 
1„ and 1 Q denote the projectors onto the respective subspaces. The Hamiltonian 
on the square A(L) is now 

H ML) = ^ h aj) + ( 29 ) 

(i,j) e£ ke\(L) 
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We decompose it as H A(L} =: H q + H u + H p , where the three terms on the right 
are defined by taking the sum over edges separately for the three interactions in 
(28) and H q contains the on-site terms in addition. Let us now assign a signature 
cr e {0,1} L to every state of our computational product basis, so that \i) (j> translates 
to cr j = 0 for all |i) e Tiff* and cr,- = 1 for all |i) e 9. By collecting computational 
basis states with the same signature, we can then decompose the Hilbert space as 




ieA 


(30) 


The reason behind this is that the Hamiltonian is block diagonal w.r.t. this decom¬ 
position, i.e. it can be written as 

H a<l> = H (r , and thus spec H A<L> = |^J spec H tr . (31) 

cr cr 


In order to identify the spectra and eigenstates coming from different signatures 
we will distinguish three cases: 


Case 1: cr = (0,..., 0). In this sector H p = H q = 0 so that we are left with H A<L) . 

Case 2: cr = (1,..., 1). In this sector H p = H u = 0 so that we are left with H A<L> . 

Case 3: cr ± (0,..., 0) but cr, = 0 for some i. Let us define A (T := {i e A(L)|cr, = 
0}. Decompose A (r into rectangles row by row, as in Figure 5 following the 
following conventions for the rectangles: 

(i) They live on the dual lattice. 

(ii) If one is on top of another of the same size (such as those marked 1 and 2 in 
the figure), we glue them together into a single larger rectangle. 

(iii) For each rectangle, we choose a neighbouring pair (i, j) such that both i and 
j are in the same row, i € A a , j i A tr . This selection is illustrated in the 
figure by black dots (joined by a link). If this is not possible because the 
rectangle is a whole row (as for the rectangle marked 3 in the figure), we 
take j in the next (or previous) possible row after applying convention (ii). 
An example is shown in the figure as a pair or linked purple dots. 

(iv) For each overlapping region between pairs of rectangles (such as those 
marked with dashed red lines in the figure), we also choose a neighbouring 
pair (i, j) with i e A (T , j £ A (T , but now i belongs to one of the rectangles (the 
larger one is a valid option) and j in the row of the other rectangle. This 
selection is illustrated in the figure by linked red dots. 
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Figure 5: Illustration of the counting explained in the main text. 

How can we choose all these pairs ( i , j) with the minimal number of repetitions? 
To repeat a pair, two things may happen. Either we have to choose the same 
vertical pair to select two different overlapping regions (located at its right and 
left) - illustrated by a pair of green dots in the figure - or we have to choose the 
same vertical pair to select one or two overlapping regions and also a rectangle 
made out of a full row. In any case, each pair is chosen not more than three times. 
For instance, the purple joint dots in the figure are chosen three times. With this 
analysis we conclude that if we denote by n Air the number of rectangles, m Ar the 
number of overlaps between them and k A<T the number of chosen different pairs 
(i, j ) with i e A a , j <£ A IT , we have that 

3 k A<r > n A<r + m A<r . (32) 

With these conventions we divide the interactions of H u into those within the 
rectangles // rectang i es , those in the overlaps of rectangles //overlaps and the rest H iest . 
We have that for any unit vector |i/v) with signature cr\ 


(*Acr I//rectangles I *Acr) — — 2 > 

(*Aoj//overlap I *A<r) — — 2/??/\ r , 
(*Acrl//restI*A<r) 0 , 


which gives 


(t/vl/ZJt/v) ^ -2 (n Aa + m Aa ) . 


(33) 
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The reason behind is that, by the same argument as in the previous section, 
MH R U ) > —2 for all rectangles R. Finally, we have trivially 

i'PaWql'Pa) > ~ 1 , (34) 

(l/falHpll/fa) > 7 k Aa . (35) 

Therefore, using (32)-(35), we get {ip (r \H ML) \ij/ (r ) - A {) (H A{L) ) > k A<r - 1 > 0. Hence, 
if we denote the part of the spectrum which corresponds to these mixed signatures 
by 5 c R we have that S > Ao(H^ <L> ). 

Putting things together, we get that the spectrum is decomposed according to 
the signatures as 


spec H A(l) = spec H A{L) U spec H A<L) u 5 (36) 

where 

MHq {L) ) = Ao(H AiL> ) - 1, S > - 1 , (37) 

and according to Lemma 9 either spec H A > 0 or H A(L) has one negative eigen¬ 
value -2, corresponding to a product ground state. Which is the case depends on 
the evolution of the encoded Turing machine (and thus becomes undecidable). □ 
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Part III 


6 Quantum Phase Estimation Turing Machine 


We will see in Section 7 how to encode an arbitrary Quantum Turing Machine 
(QTM) into a translationally-invariant nearest-neighbour Hamiltonian, in such a 
way that the local Hilbert space dimension is a function only of the alphabet size 
and number of internal states. This will allow us to encode any specific universal 
(reversible) Turing Machine in a Hamiltonian with fixed local dimension. However, 
to prove undecidability by reduction from the Halting Problem, we will need a way 
to feed any desired input to this universal TM. 

to this end, the main result of this section is to construct a family of quantum 
Turing Machines, all of which have the same alphabet size and number of internal 
states, but whose transition rules vary. For any given string, we can choose the 
transition rules such that the QTM, when started from the empty input, writes out 
the string to its tape and then halts deterministically. 

More precisely, we will prove the following result (see below for the necessarily 
basic definitions, notations, and facts about QTMs). 

Theorem 11 (Phase-estimation QTM) There exists a family of well-formed, nor¬ 
mal form, unidirectional QTMs P n indexed by n e N with the following properties: 

(i) . Both the alphabet and the set of internal states are identical for all P n ; only 

the transition rules differ. 

(ii) . On input N > \n\ written in unary, P n behaves properly, halts deterministi¬ 

cally after 0(poly(/V)2' v ) steps, uses N + 3 space, and outputs the binary 
expansion ofn (padded to N digits with leading 0’s). (Here, \n\ denotes 
length of the binary expansion ofn.) 

(iii) . For each choice of states p, q, alphabet symbols cr, r and directions D, the 

transition amplitude S(p, cr, r, q, D ) is, independently ofn, one of the elements 
of the set 



(38) 


where <p e Q. 

Part (iii) implies that the only dependence of the transition rules of P n on n is 
that implicit in ip, which is defined as the rational number whose binary decimal 
expansion contains the digits ofn in reverse order after the decimal. 

Note that the input N does not determine the output string that gets written to the 
tape; it only determines the number of binary digits in the output. The number 
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represented by that output is determined (up to padding with leading zeros) by the 
choice of the parameter n for the QTM P n . 

Theorem 11 is essentially the only inherently quantum ingredient in our result, 
and the precise properties asserted there will be absolutely crucial to the proof. For 
example, if instead of writing out the string and halting deterministically, the QTM 
did this only with arbitrarily high probability, our proof would not go through. 

It is therefore not at all clear a priori whether QTMs fulfilling these strict 
requirements exist. Since our proof depends so delicately on the precise properties 
of these QTMs, we cannot appeal to previous results. Instead, in this section we 
give a detailed and explicit construction of a QTM based on the quantum phase 
estimation algorithm, that fulfils all the requirements of Theorem 11. 

6.1 Quantum Turing Machinery 

For completeness we quote the basic definitions of Turing Machines and Quantum 
Turing Machines verbatim from [BV97], 

Definition 12 (Turing Machine - Definition 3.2 in [BV ( ']) 

A (deterministic) Turing Machine (TM) is defined by a triplet (X, Q, 6) where X is 
a finite alphabet with an identified blank symbol #, Q is a finite set of states with 
an identified initial state qo and final state qf ± qo, and 6 - the transition function 
is a function 

8:Qx^^T,xQx{L,R}. (39) 

The TM has a two-way infinite tape of cells indexed by 7L and a single read/write 
tape head that moves along the tape. A configuration of the TM is a complete 
description of the contents of the tape, the location of the tape head and the state 
q € Q of the finite control. At any time, only a finite number of the tape cells may 
contain non-blank symbols. 

For any configuration c of the TM, the successor configuration c' is defined by 
applying the transition function to the current state and the symbol scanned by the 
head, replacing them by those specified in the transition function and moving the 
head left (L) or right (R) according to 5. 

By convention, the initial configuration satisfies the following conditions: the 
head is in cell 0, called the starting cell, and the machine is in state qo. We say 
that an initial configuration has input x e (X\{#})* if x is written on the tape in 
positions 0,1,2, • • • and all other tape cells are blank. The TM halts on input x if 
it eventually enters the final state qf. The number of steps a TM takes to halt on 
input x is its running time on input x. If a TM halts, then its output is the string 
in X* consisting of those tape contents from the leftmost non-blank symbol to the 
rightmost non-blank symbol, or the empty string if the entire tape is blank. A TM is 
called reversible if each configuration has at most one predecessor. 
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Let C be the set consisting of a such that there is a deterministic algorithm that 
computes the real and imaginary parts of a to within 2“” in time polynomial in n. 

Definition 13 (Quantum Hiring Machine - Definition 3.2 in [BV97]) 

A quantum Turing Machine (QTM) is defined by a triplet (£, Q, 6) where X is a 
finite alphabet with an identified blank symbol #, Q is a finite set of states with an 
identified initial state cp and final state q t 4 cp, and 6 - the quantum transition 
function - is a function 

(gSxQxiL,^ (40) 

The QTM has a two-way infinite tape of cells indexed by 7L and a single read/write 
tape head that moves along the tape. We define configurations, initial configurations 
and final configurations exactly as for deterministic TMs. 

Let S be the inner-product space of finite complex linear combinations of 
configurations of M with the Euclidean norm. We call each element <p e S a 
superposition of M. The QTM M defines a linear operator Um '■ S —> S, called the 
time evolution operator ofM, as follows: ifM starts in configuration c with current 
state p and scanned symbol cr, then after one step M will be in superposition of 
configurations if/ = cr,c ; , where each nonzero or, corresponds to the amplitude 
6(p, cr, r, q, d) of\r) \q) \d) in the transition 6(p, cr) and c t is the new configuration 
obtained by writing t, changing the internal state to q and moving the head in the 
direction of d. Extending this map to the entire S through linearity gives the linear 
time evolution operator U M . 

Here, for convenience of programming, we will consider generalised TMs and 
QTMs in which the head can also stay still (No-movement), as well as move Left 
or Right. 

Definition 14 (Generalised TM and QTM) A generalised TM or generalised 
QTM is defined exactly as a standard TM (Definition 12) or standard QTM (Defini¬ 
tion 13) except that the set of head movement directions is { L , R, N) instead of just 
{L,R}. 

Following Bernstein and Vazirani [BV9' ], we define various classes of deter¬ 
ministic and quantum Turing Machines: 

Definition 15 (Well-formed - Definition 3.3 in [BV97]) 

We say that a QTM is well-formed if its time evolution operator is an isometry, 
that is, it preserves the Euclidean norm. 

It is easy to see (Theorem 4.2 in [BV97]) that any reversible TM is also a 
well-formed QTM where the quantum transition function 6(p, cr, q, r,d) = 1 if 
6{p, cr) = (q, r, d) for the reversible TM and 0 otherwise. 
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Definition 16 (Normal form - Definition 3.13 in [BV97]) 

A well-formed QTM or reversible TM M = (E, Q, 6) is in normal form if 

Vo-eE 6(q f ,cr) = \cr)\q 0 )\N). (41) 

In [BV9 ] the normal-form condition is S(qj, cr) = |cr) \qf) \R). However, the 
choice of final direction is just an arbitrary convention, since the machine already 
halted and never carries out those transitions. We replace R with N in the definition 
for technical reasons (see the Reversal Lemma 23, below). 

There is a specific class of QTMs, called unidirectional, that play an important 
role in the general theory developed by Bernstein and Vazirani [BV97], 

Definition 17 (Unidirectional - Definition 3.14 in [BV9 ]) 

A QTM M = (E, Q,5)) is unidirectional if each state can be entered from only 
one direction. In other words, if5(p\, cr \, T \, q, d\) and d(p 2 , cr 2 , r 2 , q, d 2 ) are both 
non-zero, then d\ = d 2 . 

Note that in a unidirectional QTM, the direction component in any transition 
rule triple |r) | q) \D) is fully determined by the state | q). Thus we can recover the 
complete transition rules from the |cr) \q) components alone, without the direction 
component. We call these the reduced transition rules, 5 r : Q x E <D 2xS . 
(Equivalently, there exists an isometry V : C e h-> C Gx|i ' w ' ,V| that maps the reduced 
transition rules to the original transition rules: V6 r (p, r) = 6(p, r).) 

The following theorem gives necessary and sufficient conditions for a (partial) 
transition function to define a reversible Turing Machine. 

Theorem 18 (Local well-formedness - Thm. B.l and Cor. B.3 in [BV9 ]) 

A TM M is reversible iff its transition function satisfies the following conditions: 

Unidirection 

Each state of M can be entered while moving in only one direction. In other 
words, if8{p\,cr\) = (ti, q, df) and d(p 2 , erf) = (r 2 , q, d 2 ) then d\ = d 2 . 

One-to-one 

The transition function is one-to-one when direction is ignored. 

Furthermore, a partial transition function satisfying these conditions can always 
be completed to the transition function of a reversible TM. 

Bernstein and Vazirani [BV97] proved that one can w.l.o.g. restrict to unidirec¬ 
tional QTMs. We therefore restrict the following quantum analogue of Theorem 18 
to the unidirectional case: 
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Theorem 19 (Quantum local well-formedness) 

A unidirectional QTM M = (X, Q, 6) is well-formed iff its quantum transition 
function 5 satisfies the following conditions: 

Normalisation 

Vp,creQxZ \\5{p, <r)\\ = 1. (42a) 


Orthogonality 

VQ?i,o-i) * (p 2 ,cr 2 )e QxT, (5 r (p u o-Q, d r (p 2 ,<x 2 )> = 0. (42b) 

Furthermore, a partial quantum transition function satisfying these conditions can 
always be completed to a well-formed transition function. 

Bernstein and Vazirani [BV9 ] proved this result for standard QTMs, that is, 
where the head must move left or right in each time step. (In fact, they prove it 
without the restriction to unidirectional QTMs; see Theorem 5.2.2 and Lemma 5.3.4 
in [BV97].) We therefore extend the proof of Theorem 19 to generalised QTM, 
which is not difficult - indeed, it is made particularly straightforward by our 
restriction to unidirectional machines. 

Proof (of Theorem 19) Let U be the time evolution operator of the QTM M. By 
definition M is well-formed iff U is an isometry, or equivalently iff the columns 
of U have unit length and are mutually orthogonal. Clearly, the normalisation 
condition specifies exactly that each column has unit length. 

Pairs of configurations whose tapes differ in a cell not under either of the 
heads, or whose tape heads are more than two cells apart, cannot yield the same 
configuration in a single time step. Therefore, all such pairs of columns are 
guaranteed to be orthogonal. The orthogonality condition imposes orthogonality of 
pairs of columns for configurations that differ only in that one is in state p\ reading 
crj while the other is in state p 2 reading cr 2 . 

It remains to consider pairs of configurations whose heads are one or two cells 
apart, differing at most in the symbol written in these cells and in their states. 
However, unidirectionality implies that any update triples that share the same state 
must share the same direction. Thus either the state or the head location necessarily 
differs for all such pairs of columns, hence these too are orthogonal. 

The final claim follows straightforwardly from the fact that the normalisation 
and orthogonality conditions imply that a partial unidirectional reduced transition 
function 6, is well-formed iff it defines an isometry on the space of states and tape 
symbols, and this can always be extended to a unitary. This fills in the undefined 
entries of 5 r by extending 6 r (q, cr) to an orthonormal basis for the space of states 
and symbols. □ 
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We will often only be interested in the behaviour of a QTM (or reversible TM) 
on a particular subset of inputs, since the machine will only be run on those. A 
proper machine is guaranteed to behave appropriately on some subset of inputs, 
but not necessarily on all possible inputs. 

Definition 20 (Proper QTM) A QTM behaves properly (or is proper) on a subset 
X of initial superpositions if whenever initialised in f e X, the QTM halts in a 
final superposition where each configuration has the tape head in the start cell 1 , 
the head never moved to the left of the starting cell, and the QTM never enters a 
configuration in which the head is in a superposition of different locations. We will 
refer to this latter property as deterministic head movement. 

Similarly, we say that a deterministic TM behaves properly (or is proper) on 
X c (E\{#})* if the head never moves to the left of the starting cell, and it halts on 
every input x e X with its tape head back in the starting cell. 

When the set X is clear from the context, we will omit it. 

Note that, for a QTM to have deterministic head movement, it is not sufficient 
that none of its transition rules 5(cr, r) produce a superposition of head directions. 
The head can also end up in a superposition of different locations because the tape 
state is in a superposition, so that two transition rules with different deterministic 
head movement apply in superposition. 

Behaving properly is not a severe restriction on classical Turing Machines. 2 
In fact, given any TM, there is always an equivalent proper TM that computes 
the same function. One way to see this is to recall that all computable functions 
are computable by Turing Machines restricted to one-way infinite tapes (see any 
standard text book on the theory of computation, e.g. Kozen [ Loz97 ]), and these 
clearly behave properly if the tape is extended to be two-way infinite. (Returning 
the head to the starting cell at the end of the computation poses no great difficulty.) 
In particular, this means that there exist proper universal Turing Machines. 

Quantum Turing Machines were originally defined in Bernstein and Vazirani 
[BV97] to have two-way infinite tapes. Indeed, those authors point out that there 
are trivial well-formed machines (such as the always-move-right machine) whose 
evolution would not be unitary on a one-way infinite tape, since the starting 
configuration would have no predecessor. However, when we come to encode 
our QTMs into local Hamiltonians, we will only be able to simulate tapes with a 
boundary. 3 To avoid technical issues, we follow Bernstein and Vazirani [BV97] in 
defining QTMs on two-way infinite tapes, but we will ensure that none of the QTMs 

'This property is called stationarity in [BV97] 

2 Nor is it a severe restriction on QTMs, but we will not need or prove this here. 

Effectively, we can only encode quantum bounded linear automata rather than QTMs, though 
there is no limit on the finite tape length we can encode. 
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(or reversible TMs) that we construct ever move their head before the starting cell. 
Thus, when encoding the QTM in a Hamiltonian, we can ignore all of the tape to 
the left of the starting cell. 

In fact, the local Hamiltonians encoding the QTMs will only be able to simulate 
the evolution on a finite (but arbitrarily large) section of tape. We will therefore 
be interested in keeping very tight control on the space requirements of all the 
reversible and quantum TMs that we construct. By carefully controlling the space 
overhead, it will then be sufficient for our purposes to simulate the evolution of the 
QTM on a finite portion of tape that is essentially no longer than the input. 

6.1.1 Turing Machine Programmin g 

For a multi-track Turing machine with k tracks and alphabet S, on track i, we will 
denote the contents of a tape cell by a tuple of symbols \cr\, cr 2 , &, cr /( ] e Si x 'Ll x 
•••xSt specifying the symbol written on each track. Similarly, the configuration 
of the tape will be specified by a tuple [c\ , c 2 ,..., c„] e 2* x x • • • x £*, where by 
convention all c, are aligned to start in the same cell (which will be the starting cell 
unless otherwise specified). We will use • to stand for an arbitrary track symbol. 
We will often describe Turing machines that act only on a subset of tracks, and 
leave the contents of all other tracks alone. In this case, we will only write out the 
states of the acted-upon tracks in the transition function; this transition function 
should be understood to be extended to the remaining tracks in the obvious way. 

As a shorthand, a transition rule on a tuple containing a • on one or more 
tracks should be understood to stand for the set of transition rules that leave the 
tracks marked • unchanged, and act as indicated on the remaining tracks. We will 
often only specify some of the elements of a transition function when defining 
a reversible or quantum TM. We will call such partial transition functions “well- 
formed” if the elements that are defined satisfy the conditions of Theorem 18 or 
Theorem 19, since by those theorems this is sufficient to guarantee that the partial 
transition functions can be completed to well-formed transition functions. QTM 
(or reversible TM) will sometimes use a finite number of auxiliary tracks. These 
are assumed always to start and finish in the all blank configuration. 

We will have frequent recourse to the following basic TM and QTM program¬ 
ming primitives from Bernstein and Vazirani [BV97], which we slightly generalise 
here to account for additional properties of the resulting QTMs that will be impor¬ 
tant to us later. These primitives allow different QTMs to be combined in various 
ways to build up a more complex QTM. 

Lemma 21 (Subroutine Lemma) 

Let Mi be a two-track, normal form, reversible TM and Mi a two-track normal 
form reversible TM (or well-formed, normal form, unidirectional QTM) with the 
same alphabet and the following properties: 
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(i) . Mi is proper on initial configurations in X\ and M 2 is proper on X 2 . 

(ii) . When started on X\, M\ leaves the second track untouched and when started 

on X 2 , M 2 leaves the first track untouched. 

(Hi). There is a state q on M\ that, when started on X\, can only be entered with 
the head in the starting cell. 

(iv). X 2 contains all the output superpositions (with qf replaced by qo) of k 
consecutive executions of M 2 started from an initial configuration in X \ for 
all 0 < k < r, where r e N U {00} is the maximum number of times that q is 
entered when M\ runs on input in X j. 

Then there is a normal form, reversible TM (or well-formed, normal form, unidi¬ 
rectional QTM) M which behaves properly on X\ and acts as M\ except that each 
time it would enter state q, it instead runs machine M 2 . 

Proof Exactly as Lemma 4.8 in [BV97]. □ 

Lemma 22 (Dovetailing Lemma) Let M\ and M 2 be well-formed normal form 
unidirectional QTMs (resp. normal form reversible TMs) with the same alphabet, so 
that Mi is proper on X\, M 2 is proper on X 2 and X 2 contains all final superpositions 
(resp. configurations) of M\ started on X\. Then, there is a well-formed normal 
form unidirectional QTM (resp. normal form reversible TM) M which carries out 
the computation of M\ followed by the computation of M 2 and that is also proper 
on X\. 


Proof Exactly as Lemma 4.9 in [BV97]. □ 

Lemma 23 (Reversal Lemma) If M is a well-formed, normal form, unidirec¬ 
tional QTM (resp. normal form reversible TM) then there is a well formed, normal 
form, unidirectional QTM (resp. normal form reversible TM) M * that reverses the 
computation of M while taking two extra time steps and using the same amount 
of space. Moreover, if M is proper on X, then M' is proper on the set of final 
superpositions (resp. configurations) of M started on X. 

Proof The proof is very similar to that of Lemma 4.12 in [BV97], We will prove 
it for QTMs, since reversible TMs are a special case of these. Consider an initial 
superposition |0 O > e X, and let \(f>i), ■ ■ • \<p n ) be the evolved sequence obtained 
by \<pi+i) = U M | (pd, where \<p n ) is a final superposition. Since M is normal form, 
|0i), • • • , \4> n ) do not have support on qo. 

Let |0') be the superposition obtained by replacing the state qf in \<p n ) with 
the initial state of the new machine q' Q . Let be the superposition obtained by 
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replacing the state qf in | (f> 0 ) with the new final state q' f . We want to construct a 
QTM M t that, when started from \<p' n ), halts on \<p[ } ) in n + 2 steps. M 1 ' will have the 
same alphabet and set of states as M, together with the new initial and final states 
q' 0 , q' f . We define the transition function S' in the following way: 

(i) . S'(q' 0 ,o-) = \ct) \q f ) \d qf ). 

(ii) For each q e Q\{qo\ and each r e E, 

d'(q, r) = S(p, cr, t, q, d q )* |&) \ p) \d p ). ( 43 ) 

PA 

(iii) S'iqo, cr) = |cr> | q' f ) \\d qo ). 

(iv) S'(q' f , cr) = |cr) \q' Q ) \N >. 

Here, for any state q , d q is the unique direction in which that state can be entered, 
and d q is the opposite direction. Since M is unidirectional. Theorem 19 implies 
that M * is a well-formed, normal form, unidirectional QTM. Given a configuration 
c in state q, n(c) is defined as the configuration derived from c by moving the head 
one step in the direction d q . n can be extended by linearity to S. 

Let us now analyse the behaviour of M t started from |</>',). By (i), AT’ maps 
| (f>' n ) in one step to n{\(p„)). Now consider (i). Since M is normal form, it maps 
superposition | if/\) to superposition \iJ/ 2 ) with no support on state q 0 if and only 
if |«Ai) has no support on state qf. Denote Qo = Q\{c/o}, Qf = Q\{qf}. If M 
takes a configuration c\ with a state from Qf with amplitude a to a configuration 
c 2 (necessarily with a state in Qo), then (i) ensures that AT’ takes configuration 
n(c 2 ) to configuration n(c.\) with amplitude a*. Let S 0 (resp. S /) be the space of 
superpositions using only states in Q 0 (resp. Qf). Since M is well-formed, the 
restriction of the evolution operator Um to S/ is an isometry into So- Hence AT’ 
implements, up to conjugation by n, the inverse of Um restricted to Um(S /). 

As an aside, note that this implies that the evolution operator of M is indeed a 
unitary and not just an isometry. Indeed, by (i) to (i) above, the evolution operator 
of AT’ is also an isometry from S 0 into S /. Since n is trivially a unitary on S, this 
implies that Um restricted to S/ is a unitary onto So. Now, since M is normal form, 
Um achieves any possible configuration with state qo by starting from the same 
configuration but with qo replaced by qf. These two facts together show that Um is 
indeed surjective and hence a unitary. 

Returning to the proof of the lemma, we have seen how (i) implies that 
executes the sequence of n steps /r(|0„_i)), • • • n(\<po)). Finally, by (i), in the (n+2)’th 
step, M’ maps n(\(f>o)) to \<p' {) ) as desired. 

Moreover, it is trivial to see that M t uses exactly the same space as M, and 
behaves properly. □ 
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Note that our way of defining normal form (Definition 16 as opposed to that in 
[BV97]) is crucial to guarantee that the reversal machine is proper. As a corollary 
of the proof, we have also shown that the evolution of a well-formed, normal 
form unidirectional QTM is a unitary operator. (The proof of this fact for non- 
unidirectional QTMs is much more involved, and can be found in [BV97]). 

6.1.2 Reversible Hiring Machine Toolbox 

It will be helpful to have a toolbox of reversible TMs which carry out various 
elementary computations. In order to satisfy the space constraints of Theorem 11, 
we will be interested in keeping tight control of the space requirements of all our 
constructions. 

Lemma 24 (Copying machine) There is a two-track, normal form, reversible TM 
Copy with alphabet 2x1 that, on input s written on the first track, behaves properly, 
copies the input to the second track and runs for time 2\s\ + 1, using l^l + 1 space. 

Proof We simply step the head right, copying the symbol from the first track 
to the second. However, we defer copying the starting cell until the end of the 
computation, so that we can locate the starting cell again. It is straightforward to 
verify that the following normal form transition function implements Copy: 



[ O ',#] 

[#,#] 

[t,t] 

qo 

([ cr,#\,q u R ) 

([#,#], q/,N) 


qi 

([o-,o-],qi,R) 

([#,#], q 2 ,L) 


qi 

([o-,o-],q f ,N) 


([ T , T ], q 2 , L ) 

qf 

([(r,m,qo,N) 

([#,#], q 0 , N) 

([t,t],c/o,/V) 


V<T, T € £/# 


This partial transition function verifies the two conditions of Theorem 18, so it is 
well-formed and can be completed to give a reversible TM. □ 

Lemma 25 (Shift-right machine) There exists a normal form, reversible TM 
Shift with alphabet 2 that, on input s behaves properly, shifts s one cell to the 
right and runs for time 2|s| + 2, using space |si + 2. 

Proof It is straightforward to verify that the following normal-form transition 
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function implements Shift: 



# cr 

40 


41’ 

(v,q 2 ,R) (v, qf R) 

42 


43 

( #,q f ,N ) (cr, < 73 , L) 

4/ 

(#,q 0 ,N) ( cr,q 0 ,N ) 


Vcr, v G £/# 

As this partial transition function verifies the two conditions of Theorem 18, it is 
well-formed and can be completed to give a reversible TM. □ 

Lemma 26 (Equality machine) There is a three-track, normal form, reversible 
TM Eql with alphabet E x E x {#, 0,1} that, on input s; t; b, where s and t are 
arbitrary input strings and b is a single bit, behaves properly and outputs s; t; b', 
where b' = ->b if s = t and b' = b otherwise. Furthermore, Eql runs for time 
2|si + 1 and uses |s| + 1 space. 

Proof Implementing a non-reversible equality-testing machine is trivial: just scan 
the head right checking if the symbols on the two input tracks match. If we 
encounter a non-matching pair, return to the starting cell and flip the output bit. If 
we reach the end of the input without encountering a non-matching pair, return to 
the starting cell and leave the output bit unchanged. 

Doing this reversibly requires more care. The problem is that the computation 
splits into two possible paths, depending on whether a non-match was encountered 
or not, and we must merge these two divergent computations back together again 
reversibly. For example, we cannot simply halt after either flipping the output bit 
or leaving it unchanged, as that would give multiple transitions into the final state 
that write the same symbol to the starting cell. 

The trick is to return to the point at which the computational paths diverged 
(either the first non-matching pair of symbols, or the end of the input) after setting 
the output bit, in order to merge the two computational paths back together, before 
returning to the starting cell again to halt. 

Because we want to avoid ever moving the head before the starting cell, there 
is also the issue of how to identify the starting cell so that we can return to it again. 
If the symbols in the first cell differ, we can set the output bit immediately and halt. 
If they are identical, we temporarily change the symbol on the second input track 
to a # to mark the starting cell. We can recover the original second-track input at 
the end of the computation by copying over the symbol from the first track. 
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The following normal-form transition function carries out this procedure: 



[cr, cr, b] 

[cr, t, b] 

[cr,#,b] 

[#,r,Z>] 

40 

([cr,#, b\,q u R) 

([ cr,T,b],q f ,N ) 

([cr,#, b],q f ,N ) 

([#, r, AO 

41 





42 



([cr,#, -4?],<73,i?) 


q ' 2 



([(r,#,Z>],<z',i?) 


43 





4s 





44 



([cr, cr, &],g/, AT) 


4/ 

([ cr, cr, b ], q 0 , N ) 

([ cr,T,b],q 0 ,N ) 

([cr,#, b],q 0 ,N ) 

([#, r, A], <?o, AO 


[cr, cr, #] 

[cr, T, #] 

[cr,#,#] 

[#,r,#] 

4i 

([cr,cr,#],q u R) 

([cr, r, #], q 2 , L) 

([<t,#,#],<7 ',L) 

([#,T,#], q' 2 ,L) 

42 

([o-,cr,#\,q 2 ,L) 




42 

([cr,cr,#],q 2 ,L) 




43 

([ cr,cr,#],q 3 ,R ) 




4a 

([cr,cr,#],q' 3 ,R) 

([(T,T,#],q 4 ,L) 

([cr,#, #],<? 4 ,L) 

([#,T,#],q 4 ,L) 

44 

([cr,cr,#\,q 4 ,L) 




4/ 

([cr,cr,#],qo,N) 

([cr,T,#],q (h N) 

([cr, #, #],^ 0 ,^V) 

([#,T,#],q 0 ,N) 



[#,#,#] 

4i 

([#,#, #],$ 2 ,L) 

42 


42 


43 

([#,#, #], ? 4,L) 

4s 


44 


4/ 

([#,#, #],q 0 ,N) 


Vcr* T6S/#, G {0, 1} 


One can verify that this partial transition function satisfies the two conditions of 
Theorem 18, so it is well-formed and can be completed to give a reversible TM, as 
required. □ 

It is somewhat easier to construct reversible implementations of basic arithmetic 
operations if the numbers are written on the tape in little-endian order (i.e. least- 
significant bit first), as it avoids any need to shift the entire number to the right 
to accommodate additional digits. We adopt this convention for all the following 
basic arithmetic machines. We do not allow numbers to be padded with leading 
0 ’s as that would allow multiple binary representations of the same number, which 
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is inconvenient when constructing reversible machines. Note that this means the 
number zero is represented by the blank string, not the string “0”. 

Lemma 27 (Increment and decrement machines) There exist normal form, re¬ 
versible TMs Inc and Dec with alphabet {#, 0,1 } that, on little-endian binary input 
n (with n > 1 for Dec), behave properly and output n + 1 or n - 1 respectively. 
Both machines run for time 0(log n) and uses at most \n\ + 2 space. 

Incrementing a binary number on a non-reversible TM is straightforward: 
simply step the head along the number starting from the least significant bit, 
flipping l’s to 0’s to propagate the carry until the first 0 or #, then flip that 0 or # to 
1 and halt. (Little-endian order avoids any need to shift the whole input to the right 
to accommodate an additional digit, should one be required.) However, making 
this procedure reversible is more fiddly. One option of course is to use Bennett’s 
[Ben73] general reversible simulation and uncomputation of a non-reversible TM, 
but this comes at a cost of polynomial space overhead. A more careful construction 
allows us to implement Inc directly, using just two additional tape cells. 

Proof (of Lemma 27) We first use the Shift machine from Lemma 25 to shift 
the entire input one cell to the right, as a convenient way of getting a # in the 
starting cell so that we can return to it later. We then reversibly increment the 
binary number written on the tape, and finish by running Shift ' (the reversal of the 
Shift machine, constructed using Lemma 23) to shift the output back one cell to 
the left. The following normal-form transition function implements the increment 
part of this procedure: 



# 0 1 

40 

(#,4i,#) 

41 

(l,q' 2 ,R) (1,42,#) (0,q u R) 

42 

(0, 43 ,22) f\,q 3 ,L) 

q' 2 

(#,43,#) 

43 

(1,44,#) 

44 

(#,4/,A0 (0,44,#) 

4/ 

(#,qo,N) (0,4o, AO (l,4o, A0 


This partial transition function verifies the two conditions of Theorem 18, so it is 
well-formed and can be completed to give a reversible TM. 

This completes the construction of Inc. To implement Dec, simply construct 
the reversal of Inc using Lemma 23. □ 

We can use the Inc construction to construct a looping primitive. (This gives a 
slightly more general version of the Looping Lemma from [BV97, Lemma 4.13], 
which also has tighter control on the space requirements.) 


53 




Lemma 28 (Looping Lemma) There is a two-track, normal form, reversible TM 
Loop2 with alphabet {#, 0,1}, which has the following properties. On input n; m, 
with n < m both little-endian binary numbers, Loop2 behaves properly, runs 
for time 0((m - n ) log m), uses space \m\ + 2 and halts with its tape unchanged. 
Moreover, Loop2 has a special state q such that on input n; m, it visits state q 
exactly m- n times, each time with its tape head back in the start cell. 

There is also a one-track, normal form, proper, reversible TM Loop which, on 
input m > 1, behaves as Loop2 on input 0; m. 

Proof Our construction closely follows the proof of Bernstein and Vazirani [BV97, 
Lemma 4.2.10], We will use two auxiliary tracks in addition to the two input tracks, 
both with alphabet {#, 0,1} and both initially blank. 

The core of Loop2 machine is a reversible TM M’ constructed out of two proper, 
normal form, reversible TMs M\ and M 2 . M\ has initial and final states q 0 ,q /, and 
transforms input n; m; x; b into n;m;x+ 1; //, where b is a bit and b' = ->b if x = n 
or x = m - 1 but not both, otherwise b' = b. M 2 has new initial and final states 
q a , q w , with qo and qj as its only other normal (non-initial or final) internal states. 
M 2 behaves as follows: 

(i) . If it is in state q a with b - 0, it flips b to 1 and enters state qo. 

(ii) . If it is in state qf with b = 0, it enters state qo. 

(iii) . If it is in state qf with b = 1, it flips b to 0 and halts. 

The transition rules for M' are constructed by deleting the qf to qo transition 
rules from Mi, and adding all the remaining Mi rules to those of M 2 . The initial 
and final states of M' are q a , q 0) . On input n;m;n;0, M' will therefore run M 2 until 
it enters the q 0 state, then run M\ until it halts and re-enters M 2 . It will continue 
to alternate in this way between M 2 and M\ until the former halts. Thus M' goes 
through the following sequence of configurations: 


„ (0 - Ml (ii) Ml (ii) 

n;m;n;0 —> n;m;n;l —> n;m;n+l;0 —» n;m;n+l;0 —> n;m;n+2;0 —> • • • 

(H) , - m „ ..... 

■ ■ • —> n;m;m-l;0 —> n;m;m;l —> n;m;m;0 (48) 


It therefore runs exactly m-n times and enters the state qf once in each run, so we 
can take qf as the special state q. 

To initialise the tracks, we dovetail a proper, reversible TM before M' which 
transforms the input n; m; #; # into the configuration n; m; n; 0. This is easily 
constructed by dovetailing the Copy machine from Lemma 24 (acting on the first 
and third tracks) with a simple reversible TM which changes the first cell of the 
fourth track from # to 0 and halts (implementing the latter is trivial). To return 
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all tracks to their initial configuration at the end, we dovetail another proper, 
reversible TM after M' which transforms the configuration n; m; m; 0 into the final 
output n; m; #; #. This is easily constructed by dovetailing the reversal Copy ' of the 
copying machine from Lemma 24 (acting on the second and third tracks) with a 
simple reversible TM which changes the first cell of the fourth track from 0 back 
to #. From Lemma 29, these initialisation and reset machines run for time 0(1 og n) 
and 0((m - n ) log m), respectively, and use \n\ + 1 and |m| + 1 space. 

It remains to construct Mi and M 2 , and show that combining them as described 
above to form M' gives a proper, normal form, reversible TM. M\ is constructed 
by dovetailing together the Eql machine from Lemma 26 (with the first and third 
tracks as its input tracks and the fourth track as its output), then the Inc machine 
from Lemma 27 (acting only on the third track), and finally another Eql machine 
(this time with the second and third tracks as its input tracks and the fourth track 
as its output). By Lemma 22 and the fact that all the constituent machines are 
proper, normal form and reversible. Mi is proper, normal form and reversible. 
From Lemmas 26 and 27, Mi runs for time 0(log m) and takes at most \m\ + 2 
space. 

The following normal form, partial transition function is essentially the same as 
the corresponding construction in Lemma 4.2.6 of Bernstein and Vazirani [BV97], 
but we can simplify slightly by exploiting the fact that we are allowing generalised 
TMs. It acts only on the fourth track, and implements a machine M 2 that clearly 
satisfies the properties (i) to (iii), above: 
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q a 

q w (#, q a , N ) 
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(hqo,N) 
(0, q 0 , N ) 
(0 ,q a ,N) 
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(0, q w , N) 
(i? qon ^0 


(49) 


M 2 is clearly normal form, and satisfies the two conditions of Theorem 18. 

Since M\ is normal form, it has no transitions into q 0 or out of q t other than the 
ones we deleted before constructing M '. All remaining internal states of M\ and 
M 2 are distinct. Thus, since M\ is reversible, the transition rules for M' also satisfy 
the two conditions of Theorem 18. M' can therefore be be completed to a normal 
form, reversible TM. Since M\ is proper, it is easy to see that M' is also proper. 

M' loops on M| m-n times (with constant time overhead and no space overhead 
coming from M 2 ). Each run of M\ takes time O(logm), so the complete Loop2 
implementation takes time 0(log m). None of the constituent TMs use more than 
\m\ + 2 space, so neither does Loop2. This completes the construction of Loop2.. 

To construct Loop, we simply remove the second input track from the Loop2 
machine, and replace the Eql machine that acts on that track with a trivial machine 
that checks whether the starting cell of the third track contains the # symbol. □ 
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The reversible incrementing, decrementing, and looping machines constructed 
above are sufficient to implement all arithmetic operations. The only ones we will 
need are addition and subtraction. These are now easy to construct. 

Lemma 29 (Binary adder) There exist two-track, normal form, reversible TMs 
Add and Sub with alphabet {#,0, l} x3 , which have the following properties. On 
input n; m with n and m both little-endian binary numbers and m > 1, Add behaves 
properly and outputs n + m; m. On input n; m with n > m, Sub behaves properly and 
outputs n- m;m. Both TMs run for time 0(m log m log n) and use max(|«|, \m\) + 2 
space. 

Proof To construct Add, we simply run the Loop TM of Lemma 28 on the second 
track, inserting for its special state the Inc machine of Lemma 27 acting on the first 
track. Since both Loop and Inc are proper and normal form, the resulting machine 
is also proper and normal form. Sub is simply the reversal of Add, constructed 
using Lemma 23. □ 

We will also need a reversible TM to convert the input from unary representation 
to binary. Again, it is not difficult to construct this using our reversible incrementing 
TM. 

Lemma 30 (Unary to binary converter) There exists a two-track, normal form, 
reversible TM UtoB with alphabet {#, 1} x {#, 0,1} with the following properties. 
On input In;# (n written in unary on the first track), UtoB behaves properly and 
outputs 1 n ; n (n written in little-endian binary on the second track). Furthermore, 
UtoB runs for 0(n 2 log n) steps and uses n + 1 space. 

Proof The basic idea is to step the head right along the unary track until we reach 
the end of the unary input string, running the Inc machine of Lemma 27 once each 
time we step right. However, Inc needs to be started with its head in the starting 
cell. So each time we run it, we need to temporarily mark the current head location, 
move the head back to the starting cell in order to run Inc, then return the head to 
its previous location and continue stepping right. 

Consider the following normal form transition function, acting only on the 
unary track: 
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(50) 




This verifies the two conditions of Theorem 18, so it is well-formed and defines a 
reversible TM. It is also easy to see that this TM is proper. 

The transition function implements a machine that steps right over the 1 ’s on 
the unary track until it reaches the end of the input, at which point it moves its head 
back to the starting cell and halts. However, each time it steps right, it temporarily 
marks the current head location by overwriting the 1 on the unary track with a #, 
returns the head to the starting cell, enters an apparently useless internal state q for 
one time step, then returns the head back to where it was before, restores the 1 on 
the unary track, and continues stepping right from where it left off. The purpose 
of this apparently pointless procedure is that, by Lemma 21, we can substitute the 
Inc machine from Lemma 27 for the state q, with Inc acting on the (initially blank) 
binary track. 

The above TM enters the q state precisely n times. Thus, when we substitute 
Inc for q, Inc will be run precisely n times, leaving the number n written on the 
binary track as required. Each iteration of Inc takes time 0(log n) and at most 
\n\ + 2 space. The step-right TM constructed above adds time overhead 0(n 2 ) and 
uses n + 1 space. Neither machine ever moves its head before the starting cell. 
Thus the overall machine satisfies the time and space claims of the Lemma. □ 

6.2 Quantum phase estimation overview 

The qubits used in the quantum phase estimation circuit can be divided into two sets: 
the “output” qubits which will ultimately contain the binary decimal expansion of 
the phase in the black-box unitary, and the “ancilla” qubits on which the black-box 
unitary is applied. In our case, the black-box unitary will be the single qubit unitary 
U v = (o )’ so ancilla set will be a single qubit. As stated above, ip will refer 
to the number whose binary decimal expansion contains the digits of n in reverse 
order after the decimal. 

The quantum phase estimation algorithm proceeds in five stages: (1) A prepa¬ 
ration stage, in which the qubits are first initialised in the |0> state, then Hadamard 
gates are applied to all the output qubits and the ancilla qubit is prepared in an 
eigenstate of U v (in our case |1)); (2) The controlled-unitary stage, during which 
control- U<p operations are applied between the output qubits and the ancilla qubit; 
(3) A stage in which we locate the least significant bit of the output; (4) The 
quantum Fourier transform stage, in which the inverse quantum Fourier transform 
is applied to the output qubits; (5) A reset stage, which resets all the auxiliary 
systems used during the computation to their initial configuration. 

We will construct QTMs for each of these stages separately, and use the 
Dovetailing Lemma 22 to chain them together. It will be useful to divide the tape 
into multiple tracks. A “quantum track” with alphabet = {#, 0,1} will store the 
qubits involved in the phase estimation algorithm. The input IN will be supplied 
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as a string of TV l’s on the quantum track. All the other tracks will essentially 
be classical; after each stage they will be left in a single standard basis state (i.e. 
neither entangled with other tracks, nor in a superposition of basis states). These 
classical tracks will be used to implement the classical processing needed to control 
the quantum operations applied to the quantum track. 

6.3 Preparation stage 

We will use the first cell of the quantum track as the ancilla qubit, and the following 
N cells as the output qubits of the quantum phase estimation algorithm. 

It will be convenient to first store a copy of the input given in unary on the 
quantum track on an auxiliary “input track”, but written in binary, i.e. to transform 
the configuration 1TV; # of the quantum and input tracks into IN; N. Running the 
unary-to-binary converter UtoB from Lemma 30 on the quantum and input tracks 
carries out the desired transformation. 

We want to initialise the N + 1 qubits that will be used in the phase estimation, 
by preparing the ancilla qubit in the state |1) and the output qubits in the |+) state. 
The first N qubits are already in the |1) state thanks to the input string. We therefore 
prepare the desired state by first temporarily flipping the state of the first qubit on 
the quantum track |#} to mark the starting cell, then stepping the QTM head right 
rotating each qubit into the |+) state, until we reach the first |#) state which we 
again rotate to |+) to initialise the N + l’th qubit. We then move the head back 
to the start location, flip the state of the first cell of the quantum track from |#) to 
|1), and halt. (The contents of all other tracks are ignored.) It is straightforward 
to verify that the following partial transition function satisfies the conditions of 
Theorem 19, so is well-formed, and implements a proper, well-formed, normal 
form, unidirectional QTM that does exactly what we want. 
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6.4 Control- stage 

The second stage of the phase estimation procedure is to apply control- U^' ' 
operations between the n’th output qubit and the ancilla qubit (see Figure 6). 
Constructing a QTM that implements this is more complex. The basic idea is to 
supply the QTM with an internal state q which causes the U v rotation to be applied 
to the quantum track at the current head location. We then construct classical 
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control machinery which iterates over the output qubits, and loops on q to apply 
U v for a total of 2” _1 times. 


10 ) 

10 ) 

10 ) 

10 ) 

ID 



+ e 2^-V)|i > 

+ e 2ni{l2(p) |1) 

+ |i) 

+ e lnilf 11) 


Figure 6: The first stage of the quantum phase estimation circuit for tp (cf. Fig. 5.2 
in [NC00]). 


We will use three auxiliary tracks. One, with alphabet X m = {t, c, m 0 ,mi,#}, 
will be used to mark the position of the current control and target qubits. (The 
m 0> i states will be used to temporarily store an auxiliary qubit on the mark track of 
the starting cell.) The other two auxiliary tracks, with alphabet X/ = {0,1, #}, will 
constitute the work tapes of two different reversible looping TMs from Lemma 28. 

6.4.1 cU k machine 

We will need a QTM which applies the control- U v operation k times. We give 
a construction for an arbitrary controlled single-qubit unitary U, as this will be 
useful later. 

Lemma 31 (Controlled -U QTM) For any single-qubit unitary U, there exists a 
well-formed normal form unidirectional QTM cU k with the following properties. 
The QTM has three tracks: a looping track with alphabet X/ = {0,1,#}, a mark 
track with alphabet X m = {t,c, mo,mi,#}, and a quantum track with alphabet 
X 9 = {0,1,#}. The input consists of a number k > 1 written in binary on the 
looping track in little-endian order, a configuration containing a single t and a 
single c within the first n tape cells on the mark track (and all other cells blank), 
and an arbitrary n-qubit state on the quantum track. 

On such input, the QTM applies the control- U operation k times between the 
control and target qubits on the quantum track marked by c and t, and then halts, 
having run for time 0(kn + k log k), used at most max(«, |k|) + 2 space, behaving 
properly and leaving the configurations of the looping and mark tracks unchanged. 
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Proof It will be convenient to first run the Shift machine from Lemma 25 acting 
only on the mark and quantum tracks, to shift this part of the input one cell to the 
right. This produces a # on the mark and quantum tracks of the start cell, which 
we can use later to return to this cell. The remainder of the construction returns 
the quantum and mark tracks of the start cell to the # state, so we can run the 
corresponding Shift' machine at the very end (where Shift' is the reversal of M\ 
constructed using Lemma 23) to shift these tracks back one cell to the left, so that 
the final output is correctly aligned. These shift operations take 0(n ) time and use 
n + 2 space. 

We construct the core of the cU k QTM out of two simpler machines, M\ and 
M 2 , dovetailed together in the sequence M \, M 2 , M\ (where M\ is the reversal of 
Mi). Machine Mi scans right until it encounters aeon the mark track, at which 
point it applies a CNOT between the quantum track and an internal qubit, returns 
the head to the starting cell, applies a CNOT between the internal qubit and an 
auxiliary qubit on the mark track, and halts. Clearly, Mi is proper, runs for time at 
most 0(n) and uses at most n + 1 space. 

The following well-formed, normal form, partial quantum transition function 
implements Mi (which acts only on the mark and quantum tracks). Note that 
one of the internal states of Mi comes in two varieties, indicated by a superscript 
j = {0,1}, which are used to temporarily store a qubit in the internal state of the 
machine: 
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(52) 

Machine M 2 loops k times, where k is specified by the number written in binary 
on the looping track, applying a control- U operation between the auxiliary qubit 
stored on the mark track of the start cell and the target qubit, and then halts. We 
construct M 2 using the reversible looping TM Loop from Lemma 28, and inserting 
for its special state a QTM M' that is very similar to Mi. The M' machine first 
applies a CNOT between the auxiliary qubit stored in the mark track of the start 
cell and an internal qubit. It then scans right until it finds the t, and applies a single 
control- U operation between the internal qubit and the target qubit. Finally, it 
moves the head back to the starting cell, and applies another CNOT between the 
internal qubit and the auxiliary qubit, before halting. 

If Ujj denotes the i, j’th matrix element of U, then the following well-formed, 
normal form, partial quantum transition function (which acts only on the mark and 
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quantum tracks) implements M'\ 
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This machine is proper, never alters the mark track, and, for given mark track 
configuration, always runs for the same number of time steps. So substituting it 
in the looping machine produces a proper QTM M 2 . Furthermore, since M' takes 
0(n) time and at most n + 1 space, M 2 runs for time Oikn + klogk), uses at most 
max(n, \k\) + 2 space, and never moves the head before the starting cell. 

The only qubits on which the overall QTM acts are the two qubits in the 
quantum tracks of the cells marked by t and c, the auxiliary qubit stored on the 
mark track (let’s label this m), and the two internal qubits stored in the internal 
states qi of M\ and q\ 2 of M 2 (let’s label these auxiliary M\ and M 2 internal qubits 
by 01 , 2 )- Qubits m, ci\ and a 2 are initially in the |0) state. Apart from classical 
processing to move the head into the correct location, M\ just applies a CNOT 
between the c and a\ qubits, followed by a CNOT between the ci\ and m qubits. 

Similarly, M' applies a CNOT between the m and a 2 qubits, a control- U 
operation between the a 2 and t qubits, and a final CNOT between the m and a 2 
qubits. Letting cU denote the control- U gate, and recalling that the a 2 qubit starts 
off in the |0) state, the overall effect of this is: 

CNOT ma2 cU a2t CNOT ma2 \ if,) mt |0) fl2 = cU W) mt |0) fl2 , (54) 

i.e. M' effectively applies a single control- U operation between the m and t qubits. 
M 2 repeats the M' operations k times, so the overall action of M 2 is to apply a 
control- U operation k times between the m and t qubits. 

The time-reversal of Mi simply applies its CNOTs in the reverse order. So 
dovetailing M\,M 2 , M\ carries out the operation 

CNOT cai CNOT aim cU k mt CNOT aim CNOT cai \<jj) ct |0> ai |0> m = cU k ct N> ( , |0> ai |0) m . 

(55) 

Thus the overall action of the entire QTM is to apply a cU k operation between 
the control and target qubits, as required, leaving the looping and mark tracks 
back in the configuration they started in. Furthermore, the QTM runs for time 
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0(kn + k log k), uses at most max(«, |fc|) + 2 space, and is well-formed, normal form, 
unidirectional and proper, as claimed. □ 

6.4.2 Iterating over the control qubits 

We want to apply a cU*" ' operation between the n’th output qubit and the ancilla 
qubit, for each of the N output qubits. (Note that the qubits are stored in reverse 
order on the quantum track, so the 1st output qubit is stored in the N + l’th cell 
of the quantum track and the N’th output qubit is stored in the 2nd cell.) As we 
will need to use a second looping machine to iterate over the qubits, we refer to the 
looping track for the cU k machine as the u cU looping track”, and the track for the 
second looping machine used here as the “outer looping track”. 

We will need a reversible TM which, on suitable mark and cU looping track 
configurations (cf. Lemma 31), moves the c marker on the mark track one cell to 
the left, and doubles the value k written on the cU looping track. It is convenient 
(though slightly less efficient) to divide the implementation into two parts, dove¬ 
tailed together: Step c which shifts the c marker one cell to the left, and Dbl which 
doubles the value on the cU looping track. The following well-formed, normal 
form, partial transition function implements the reversible TM Step c (which only 
acts on the mark track): 
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Vi €{0,1} 

Doubling a number in binary simply appends a 0 onto the binary representation 
of the number. The following well-formed, normal form, partial transition function 
for Dbl (which acts only on the cU looping track) does precisely this, taking 
advantage of the fact that k is always a power of 2 (so, in little-endian order, always 
consists of a string of 0’s followed by a single 1) to simplify the construction: 
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We are now in a position to implement the complete controlled- U v stage of 
the quantum phase estimation algorithm. We initialise the states of the auxiliary 
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tracks, by first running a simple reversible TM that changes the [#, #] in the first 
cell of the mark and cU looping tracks into [to, 1], then returns to the starting 
cell and halts. (Constructing a reversible TM for this is trivial.) We then dovetail 
this with a proper, normal form reversible TM that scans to the end of the output 
qubits, changes the # on the mark track into a c, and halts. The following normal 
form, partial transition function (which acts only on the mark and quantum tracks) 
accomplishes this: 
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On the mark track, this leaves a t 0 and c over the ancilla qubit and the first output 
qubit, 1 respectively, and blanks everywhere else. The configuration prepared on 
the cU looping track corresponds to the number 1 written in binary. 

Next, we copy N from the input track to the outer looping track using the Copy 
machine from Lemma 24. We then run a reversible looping TM from Lemma 28 
which uses this track as its input track (so it will loop N times in total). For the 
special state of this looping machine, we substitute a QTM which dovetails the 
cU k machine from Lemma 31 with the Step c and Dbl machines constructed above. 

On appropriate input, the cU k machine from Lemma 31 runs for a number of 
steps which depends only on the classical part of the input supplied on the mark 
and looping tracks, and not on the quantum state in the quantum track. The Step c 
and Dbl reversible TMs don’t touch the quantum track at all. So as long as the 
mark and cU looping tracks are always in a suitable configuration (cf. Lemma 31) 
before the cU k machine is run in each iteration, the outer looping machine will be 
proper. 

Now, we already initialised the cU looping track to k = 1, above, and marked 
the ancilla qubit as the target and the first output qubit as the control. So the 
initial configuration of the mark and cU looping tracks is suitable input for a cU k 
machine. When cU k is run in the first iteration of the outer looping machine, it 
applies act/ 1 = cU between the first output qubit and the target qubit. Step c and 
Dbl machines then run, which shifts the control marker one cell to the left onto the 
second output qubit and doubles k to 2, ready for the next iteration. In the second 
iteration, the looping machine therefore applies act/ 2 between the second output 

Recall that the first output qubit is the last one on the quantum track 
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qubit and the ancilla, and doubles k to 4. The looping machine goes through a total 
of N iterations, each time feeding suitable input to the cU k machine to make it 
apply a cU 2 " between the n’th output qubit and the ancilla, before halting with 
N written on the outer looping track, a to and c in the first and second cells of the 
mark track, and N + 1 written on the cU looping track. 

Thus the outer looping machine is well-formed, normal form unidirectional and 
proper. It runs for a total of N iterations with overhead NlogN (cf. Lemma 28). In 
the n’th iteration it runs the cU k machine from Lemma 31 once with k = 2” _1 , and 
runs the Step c and Dbl machines once each. The n’th iteration of the cU k machine 
takes time 0(N + 2" _1 log2" _1 ), and the Step c and Dbl machines always take time 
0(n). So the outer looping machine runs for a total time OiN 2 2 N ). 

Finally, we can uncompute (reset) the auxiliary tracks by running a reversible 
TM which transforms the configuration N; N; fac; N + 1 left on the input, outer 
looping, mark, and cU looping tracks, into N\ #; #; #. We do this by running the 
reversal Copy' of the copying machine from Lemma 24 on the input and outer 
looping tracks to erase the outer looping track, decrementing the cU looping track 
and running Copy ' again to erase that track, and running a trivial reversible TM 
that changes t 0 c on the mark track into ##. 

Putting everything together, the control- U stage runs for a total time 0(N 2 2 n ) 
and requires space N + 3 (2 more cells than the number of quits, which is N + 1). 
When dovetailed with the preparation stage from Section 6.3, the overall QTM we 
have constructed implements the control- U circuit of Figure 6. If <p k denotes the 
Cth digit in the binary decimal expansion of <p, then this prepares the state [NCOO] 
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on the N output qubits (ordered as they are on the quantum track). All other tracks 
are blank, except for the input track which still has the input N written on it as a 
little-endian binary number. 


6.5 Locating the LSB 

The final stage of the quantum phase estimation algorithm is to apply the inverse 
quantum Fourier transform to the quantum state generated by the control -U stage. 
The structure of the quantum circuit for the QFT is rather reminiscent of that of 
the control- U stage. It again involves applying a “cascade” of control- U 2 " gates 
between pairs of qubits, which we already know how to implement! Only now 
there is a new cascade starting from each output qubit (see Figure 7). 

However, the phase 2~ N of the control- U 2 - v rotation needed in the QFT circuit 
depends on the total number of qubits N that the QFT is being applied to. If we 


64 


simply implemented the inverse QFT circuit directly on all N output qubits, the 
entries of the transition function in our cU k machine from Section 6.4 would have 
to depend on the input N, which is not allowed. It is not at all obvious whether the 
inverse QFT circuit on N qubits can be implemented exactly when N is given as 
input. 1 

On the other hand, the entries of the transition function are allowed to depend 
on the phase tp that we are estimating - indeed, they necessarily do so, since 
the output of the P n QTM in Theorem 11 depends on n (or equivalently on tp). 
Rather than implementing the QFT on N qubits, we take a different approach. We 
supply our QTM with the U 2 -m gate, and implement the inverse QFT on M qubits, 
independent of the input N. (I.e. just enough qubits to hold all the digits of the 
binary decimal expansion of tp.) However, for this to work, we must first identify 
which output qubit holds the least significant bit (LSB) of tp, so that we know which 
qubits to apply the |</>|-qubit inverse QFT to. The role of the input N is merely to 
provide us with an upper-bound on \tp\, which allows us to identify the LSB in finite 
time and space. 

Recall that the control- U stage prepared the state [NCOO] 
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on the N output qubits (where tp k denotes the &’th digit in the binary decimal 
expansion of tp). 

The least significant bit of tp is the M’th bit, and by assumption (cf. Theorem 11) 
\tp\ < N. So the M’th bit of tp is 1 and the \tp\ + 1... TV’th bits are 0. Thus the last 
N - \tp\ output qubits are in the |+) state, and the M’th qubit, which corresponds to 
the least significant bit of tp, is in the |—) state. (Recall that the output qubits are 
stored on the quantum track in reverse order.) 

Therefore, if we construct a QTM that steps right, applying Hadamard rotations 
to the output qubits in the quantum track until it obtains a |1) and halts, we will be 
able to locate the least significant bit of tp. More precisely, this will rotate the first 
N - \tp\ output qubits into the |0) state, and halt with the N - \tp\ + l’th qubit (which 
corresponds to the M’th bit of tp) rotated into the |1) state. The least significant bit 
of tp is therefore identified by the first |1) on the quantum track after this machine 
has finished running. 

The following well-formed, normal form, partial quantum transition function 
(which acts only on the quantum track) implements a proper QTM that steps along 
the quantum track, applying Hadamards until it obtains a |1). It can be verified that 

1 The universal QTM construction of Bernstein and Vazirani [BV97] can approximate the QFT 
on N qubits to arbitrary precision, but it does not implement it exactly. 
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it obeys the three conditions of Theorem 19, so can be completed to a well-formed 
QTM: 
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(Note that the qubit in the starting cell of the quantum track is still in the state |1) 
from the preparation and control -U stages. As we will not need this qubit again, 
this machine sets it to |#) to mark the starting cell and leaves it in that state when it 
halts, for later convenience.) 

The configuration written on the quantum track after this machine has finished 
is the N + 1-qubit state: 


|#) |0 ) N - •'?' |0> |1> -4ir (|0> + e 2ni0 -^- m ^ |1>) (|0> + |i>)... 

... ^| 0 ) + ( 62 ) 


2 2 


Note that the first N - \(f>\ + 2 cells of the quantum track are not entangled with the 
rest of the track. 

It will be helpful for later to record the number of output bits that we are 
skipping, i.e. to compute N - \ip\ and store the result on a separate auxiliary “count 
track”. We can do this using a construction that is very similar to the unary-to- 
binary converter of Lemma 30. Specifically, we will construct a reversible TM 
Count that steps right over |0)’s until it finds a |1) on the quantum track, running 
the Inc machine once each time it steps right. (This means first returning the 
head to the starting cell, running Inc, then returning the head back to where it was 
previously.) 

Consider the following partial transition function, which satisfies the conditions 
of Theorem 18 so can be completed to a proper, normal form, reversible TM: 
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(63) 







This machine steps right over 0’s until it encounters a 1 (or reaches the end of 
the input), at which point it moves its head back to the starting cell and halts. (It 
assumes that the first cell is marked with a #.) However, each time it is about to 
step right, it first temporarily marks the current head location by overwriting the 0 
with a #, returns the head to the starting cell, enters an internal state q for one time 
step, then returns the head back to where it was before, resets the quantum track 
to 0, and continues stepping right from where it left off. By Lemma 21, we can 
substitute for the state q the Inc machine from Lemma 27 (acting on an initially 
blank “count track”). If we run the above TM on the quantum track, it will step 
right precisely N - \ip\ times before finding the first |1), so it will halt with N - \ip\ 
written to the count track. 

The configuration left on the tape now consists of the little-endian binary 
number N written on the input track, the little-endian binary number N - \(p\ written 
on the count track, and the N + 1-qubit state from (62) written on the quantum 
track. 

As the qubits storing the N - \<p\ trailing 0’s of the TV-digit binary expansion of 
(p now play no further role, it is convenient to shift the starting cell of subsequent 
machines to the tape cell containing the LSB. To this end, we copy the values N 
and N - \<p\ from the input and count tracks onto new auxiliary tracks. We then 
use a Loop machine from Lemma 28, which uses the original count track as its 
input and work track, to run a Shift machine from Lemma 25 acting on the tracks 
holding the copies of the count and input tracks. Thus the Shift machine will run 
a total of N — \<p\ times, thereby shifting the binary strings on the copies so that 
they start in the N - M’th cell - the one containing the LSB of \(f>\ in the quantum 
track. We then run a trivial TM that steps the head right over the initial string of 
|0)s on the quantum track, until it encounters the first |1) identifying the LSB, and 
halts with the head at this LSB cell. When we refer to the input and count tracks in 
the following section, we mean the copies shifted to start at the LSB cell. (At the 
very end of the computation, we will uncompute the shifted copies of the count 
and input tracks and step the head left back to the original starting cell, to ensure 
the overall QTM remains proper.) 

Now, we have been careful to ensure that the head is never moved to a location 
before the starting cell in any of the reversible and quantum TMs that we have 
constructed. Furthermore, the section of the complete tape configuration located 
before the LSB cell is unentangled with the rest (see (62), and note that all tracks 
other than the quantum track are in classical configurations). Thus, if we start 
any of our TMs in the LSB cell, it acts as if its input were the section of the tape 
configuration located after (and including) the LSB cell. We can therefore ignore 
the tape configuration of the first N - \<p\ cells in the following section, and restrict 
our attention to the following \<p\ cells. 
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6.6 QFT stage 

We are now ready to apply the |<£>|-qubit inverse QFT to the |^|-qubit state stored on 
the quantum track. The inverse QFT circuit on |^| qubits consists of cascades of 
cC/*-m gates, very reminiscent of the control-L/ stage we have already implemented 
(see Figure 7). The construction will therefore be similar. In the following, as we 
have shifted the starting cell, we relabel the output qubits and refer to the LSB 
qubit as the 1st qubit, the last one as the M’th qubit. 



hi) 

\h) 

\jN-l) 

\jN-l) 


Figure 7: The inverse QFT stage of the quantum phase estimation circuit 
(cf. Fig. 5.1 in [NCOO]). 


In each cascade, we want to apply a cUim-n = (c[/ 2 -m) 2M+ gate between the 
m’th and n’th qubit (m < n). Once again, we will use mark and cU looping tracks 
to hold the input to a cU k 2 _ w machine from Lemma 31. However, the main loop 
will now consist of two nested loops: an outer loop to iterate the control qubit m 
of the CC/J7 gate over each output qubit, and an inner loop to iterate the target 
qubit n over qubits m + 1 through \< p \. 

We first run a TM that initialises the mark and cU looping tracks, so that 
the mark track contains the configuration ct in the first two cells, and the cU 
looping track contains a string of \ip\ - 1 0’s followed by a 1. Note that this initial 
configuration of the cU looping track is the little-endian binary representation 
of the number 2^ l_1 . (Constructing a proper normal form, reversible TM that 
implements all of this is an easy exercise.) 

We also use the Inc machine from Lemma 27 to increment the number N - \< p \ 
stored on the count track to N - \( p \ + 1. We then run the Sub machine from 
Lemma 29, with the input and count tracks as the input tracks and the inner looping 
track as the output track, to write the number N - (N - \<p\ + 1) = \np\ - 1 to the 
inner looping track. 


Inner loop The inner looping machine is very similar to the main looping ma¬ 
chine from the control-!/ stage. We use a Loop machine from Lemma 28, with the 
inner looping track as its input and work track. In each iteration, we first run a 
TM to divide the value on the cU looping track by 2. (This is simply the reversal 
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Dbl ' of the machine constructed in Section 6.4.2.) We dovetail this with the cC/*_ M 
machine from Lemma 31. Finally, we dovetail this with a TM that moves the target 
marker t on the mark track one cell to the right. (This can be implemented by 
taking the reversal StepJ of the machine from Section 6.4.2 and replacing c with t 
in its transition rules to give a Step,' machine.) 

Assume that the inner looping track is initialised to \ip\ - m, the control and 
target markers on the mark track are initially over the m’th and m + l’th qubits, 
and the cU looping track is initialised to 2 ,<pl . Then the effect of this inner looping 
machine is to apply c[/ 2 m -« operations between the m’th and n’th qubit for all 
m < n < \ip\ - the cascade of cU’s starting from the m’th qubit in Figure 7. 

This inner looping TM leaves the mark track in the configuration with a c 
in the same cell that it started out in, a t in the \np\ + l’th cell (i.e. the cell after 
the last output qubit), and 2 m_1 written on the cU looping track.(The latter is the 
configuration consisting of an initial string of 0’s, followed by a 1 in the cell 
containing the c on the mark track.) 

Outer loop For the outer looping machine, we use a Loop2 machine from 
Lemma 28 running on the input and count tracks (which hold the numbers N 
and N - \<p\ + 1, respectively). By Lemma 28, this machine will therefore loop 
\<p\ - 1 times. 

In each iteration, we first run the inner looping machine. Note that for the first 
iteration, the auxiliary tracks are already initialised as assumed above for the value 
m = 1. We dovetail the this with a simple QTM that applies a Hadamard operation 
to the current control qubit (the one marked by a c on the mark track). 

We then run a TM that changes the 1 on the cU looping track to a 0, shifts 
the c on the mark track one cell to the right, steps right, and changes [#, #] on the 
mark and cU looping tracks to [t, 0]. The machine then steps right, changing # to 
0 on the cU looping track as it goes, until it reaches the end of the output qubits. 
Whereupon it changes # on the cU looping track to 1, steps right, and changes the 
t on the mark track to #. Finally, it returns to its starting cell and halts. (Again, by 
now, constructing a proper, normal form, reversible TM for this is straightforward.) 
We dovetail this with the Dec TM, acting on the inner looping track. 

The effect of all this is to reset the configuration of the auxiliary tracks, ready 
for the next iteration of the inner loop. The control marker c is shifted to the 
next qubit along, qubit m say, and the target marker t is placed over the adjacent 
m + l’th qubit. The cU looping track is reset to 2^ l_1 , and the inner looping track 
is decremented to \<p\ - m, as required. In other words, the auxiliary tracks are 
initialised to the desired configuration for the new value of m. 

Thus the outer looping machine runs the inner looping machine \ip\ - 1 times 
for each of the output qubits m = 1... \(p\ - 1. Each time it runs, the inner looping 
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machine applies the desired cascade of cU operations between the m’th and n’th 
qubits for all m < n < \<p\. The outer looping machine then applies a final Hadamard 
operation to the m’th qubit, before moving onto the next qubit. One can therefore 
see that the overall effect is to apply the inverse quantum Fourier Transform circuit 
to the \ip\ output qubits. 

6.7 Reset Stage 

To reset all the auxiliary tracks, we dovetail the outer looping TM with a sequence 
of TMs that uncompute all the configurations we previously prepared on these 
tracks. 

The first of these TMs changes the ct left on the mark track over the final two 
qubits to ##, and changes the final configuration 1; 2 lv,|_1 of the inner looping and 
cU looping tracks to the blank configuration. (Note that 2^ l_l is the configuration 
consisting of a leading string of 0’s followed by a single 1 on the final output qubit, 
so is straightforward to reset without any arithmetic computations.) 

To reset the final configuration N; N - \<p\ + 1 of the input and count tracks, we 
start by decrementing the count track to N; N - \ip\ (using Dec). Recall that we 
redefined the location of the starting cell before implementing the inverse QFT 
machine, so that the starting cell became the cell containing the LSB of tp. We also 
shifted the copies of the input and count tracks accordingly (Section 6.5). We now 
step the head back to the original starting cell (which is easily located as we left a 
# written there on the quantum track, cf. Section 6.5). We can then run the reversal 
Shift' (constructed using Lemma 23) of the machine used in Section 6.5, to shift 
the copies of the input and count tracks used in the QFT stage back left again to 
the original starting cell, and run the reversal Copy ' of the copy machines to erase 
the copies of the input and count tracks. 

To erase the original count track, which contains N — \<p\, we must run the 
reversal Count' that we constructed in Section 6.5. (The first N - \tp\ qubits remain 
in the |0) state and the LSB qubit in the |1) state, so Count' will indeed uncompute 
N - \ip\.) Finally, to erase the original input track, we first use the Shift' machine 
from Lemma 25 on the quantum track, to shift the entire final state of the output 
qubits left. Note that this will work because the first qubit (the ancilla qubit of the 
control-L stage) was set to # in Section 6.5. We can then run the reversal UtoB' 
of the unary-to-binary machine from Lemma 30, acting on the quantum track and 
input track, to uncompute the binary conversion of the unary input encoded in the 
length of the qubit state. Note that the UtoB machine only checks whether the 
symbols on the unary track are # or non-#. So the fact that the configuration of the 
quantum track is no longer necessarily a string of N |l)’s does not matter; all that 
matters is that it is a string of N non-|#)’s. 

The end result of all this is to reset all the auxiliary tracks to the blank state, 


70 


leaving the output of the inverse QFT circuit stored in the first N cells of the 
quantum track. 

6.8 Analysis 

By dovetailing together the preparation stage (Section 6.3), control- U v stage (Sec¬ 
tion 6.4), inverse-QFT stage (Section 6.6) and reset stage (Section 6.7), we have 
succeeded in constructing a family of well-formed, normal form, unidirectional 
quantum Turing Machines P n that behave properly on input N > \<p\ written in 
unary, and implement the quantum phase estimation algorithm on N qubits for 
phase (f. By construction, P„ satisfies part (i) of Theorem 11. 

Now, we know [NCOO] that the quantum phase estimation algorithm outputs the 
exact binary decimal expansion of the phase, as long as we run it on enough qubits 
to store the entire binary decimal expansion of the phase. Thus for N > \ip\ = \n\, 
the P„ writes out the binary decimal expansion of (p to N bits (in little-endian order, 
padding with 0’s as necessary). We choose <p to be the rational number whose 
binary decimal expansion contains the digits of n in reverse order after the decimal, 
where n indexes P„. So our QTM P„ implements the computation claimed in 
part (ii) of Theorem 11. 

We were careful throughout to keep tight control on the space requirements 
of the reversible and quantum TMs that we used to construct P n . In fact, none of 
them used more than N + 3 space. Finally, all steps of the computation take time 
0(poly N), except the cU k and cU k computations, which take time 0(2 N ). Thus the 
overall run-time is 0(poly(/V)2' v ). Thus P n fulfils the space and time requirements 
in part (ii) of Theorem 11. 

The last thing we must check in order to finish the proof of Theorem 11 is 
part (iii). But the partial transition function defined so far for P n satisfies part (iii). 
It is trivial to check that one can complete this to a full transition function, whilst 
still satisfying part (iii). Indeed, by Theorem 19 the problem is equivalent to 
completing an orthonormal set of vectors with coefficients in 
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to a full orthonormal basis with coefficients in S. Since any normalised vector with 
coefficients in S must be proportional to a vector in the canonical basis {<?/}/ or of 
the form -^(e v ± e k ), j 4 k, the result is immediate. This completes the proof of 
Theorem 11. 

Note that the QTM we have constructed only implements the computation 
correctly when supplied a valid upper-bound on the number of digits in the binary 
decimal expansion of the phase. If the input is not actually a correct upper bound, 
we make no claim about the behaviour of the QTM. (This will be significant later, 
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when we come in Section 9 to bound ground state energies of a Hamiltonian 
constructed out of this QTM.) 
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7 Encoding QTMs in local Hamiltonians 


In this section, we give a self-contained analysis of a local Hamiltonian construction 
that generalises that of Gottesman and Irani [GI09], in order to prove a general 
result about local Hamiltonian encodings of arbitrary quantum computations. 

A full, rigorous analysis of the resulting Hamiltonian will then prove the main 
result of this section: 

Definition 32 (Computational history state) A computational history state \fi) CQ e 
*Hc ® ' Hq is a state of the form 
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where j|f)} is an orthonormal basis for‘H q* an d \fii) = W\ = \ £7; |«Ao ) for some initial 
state |(A 0 ) e 'Hq and set of unitaries Ui € 'B('Hq). 

'He is called the clock register and 'Hq is called the computational register. IfU, 
is the unitary transformation corresponding the t’th step of a quantum computation, 
then \f t ) is the state of the computation after t steps. We say that the history state 
\fi) encodes the evolution of the quantum computation. 

Theorem 33 (Local Hamiltonian QTM encoding) 

Let = (span{@, ©} © <D C ) <g> (span{©, ©} © be the local Hilbert space of 
a 1-dimensional chain of length L. 

For any well-formed unidirectional Quantum Turing Machine M = (2, Q, 6) 
and any constant K > 0, we can construct a two-body interaction h e B(C d ® 
C rf ) such that the 1-dimensional, translationally-invariant, nearest-neighbour 
Hamiltonian H(L ) = Yn=\ e B('H(L)) on the chain of length L> K + 3 has 

the following properties: 

(i) . d depends only on the alphabet size and number of internal states of M. 

(ii) . h > 0, and the overall Hamiltonian H(L) is frustration-free for all L. 

(iii) . If we call 'H(L - 2) := (C c ) 0L_2 ® (C e ) 0L_2 =: 'He ® 'Hq and define S/, r c 

'H = span(|@) ® |©))i ® FKL - 2) ® span(|@) ® |©)) L , the unique ground 
state of H{L)\s hr is a computational history state encoding the evolution of M 
on input corresponding to the unary representation of the number L- K - 3, 
running on a finite tape segment of length L - 3. 

Moreover, if M is proper on input given by the unary representation of L - K - 3, 
then: 
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(iv) . The computational history state always encodes D(|S x Q\ L ) time-steps. If 

M halts in fewer than the number of encoded time steps, exactly one \iff) 
has support on a state |T) that encodes a halting state of the QTM. The 
remaining time steps of the evolution encoded in the history state leave M’s 
tape unaltered, and have zero overlap with |T). 

( v ) . If M runs out of tape within a time T less than the number of encoded time 

steps (i.e. in time-step T + 1 it would move its head before the starting cell or 
beyond cell L- 3), the computational history state only encodes the evolution 
of M up to time T. The remaining steps of the evolution encoded in the 
computational history state leave M’s tape unaltered. 

(vi) . Finally, if M satisfies part (in) of Theorem 11, then h has the following form 

h = A + (e in<p B + e in2 ~ M C + h.c.), (66) 

with B,C e £>(<C d ® C rf ) independent ofn and with coefficients in 7L, and 
A e !B( C rf <8>C rf ) Hermitian independent ofn and with coefficients in 7L + -^=Z. 

Though we will not need it for our purposes, it is not difficult to prove that 
the next-highest eigenstate of H(L)\ Sbr has energy 0(1/T 3 ), where T is the total 
number of time-steps encoded in the computational history state. 

7.1 Preliminaries 

As our construction draws heavily on [GI09], we will follow their notation and 
terminology, which we summarise here. We divide the chain into multiple tracks: 



Track 1: Clock oscillator 



Track 2: Counter TM head and state ... 


© 

Track 3: Counter TM tape 

© 


Track 4: QTM head and state 



Track 5: QTM tape 



Track 6: Time-wasting tape 
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The local Hilbert space at each site is: 

:= (span{|s, p, cr)} © span{|©), |©)}) ® (spanflg, r, y>} © span{|©}, |©>}) 

5 e {©,©,©,©,0,0,®}for/ e {1.ATJ, 

p € P' U {C , ®} where P := P L U P N U P R and P' := P U P' R , 
q e 0' U P U P' L U {r x } U {O, ©} where 0 := Q L U Q N U 0«, 0' := 0 U 0^ 
and * € 0 
<T S X, 

r € E where E := {h,#,0,...,£- 1} with £ = |cr x 0|, 
y e E U {l- ? } where q € 0'. 

(67) 

X is the tape alphabet of our given QTM M. E is the alphabet of the counter TM. 
P L , Pn , Pr are the sets of internal states of the counter TM that can be entered by 
the TM head moving left, not moving, or moving right (respectively). The states 
p' e P' R duplicate the states p e P R , and p' e P' L duplicate those in P L . Similarly 
for the internal states Q L , Q N , Q R of the QTM, with q' e Q' L duplicating the states 
q € Q l . Recall that by the unidirection property of reversible and quantum Turing 
Machines (Theorems 18 and 19), these sets are disjoint. The O and © Track 2 
and 4 symbols are used for cells that do not currently hold the head. 1 The bracket 
states ©, © are the special boundary states |@), |©) from Theorem 33 that will 
only appear at the ends of the chain. The role of the Track 1 states is described 
in Section 7.2. That of the Track 6 “time-wasting” tape, as well as the role of the 
additional Track 4 P U P' L states, will become apparent in Section 7.6.2. 

This set of states defines a standard basis for the single-site Hilbert space ( H. 
The product states over this single-site basis then give a basis for the Hilbert space 
< H® L of the chain. We call these the standard basis states. 

Definition 34 (Bracketed state) We call standard basis states that match the 
regular expression ©O*© bracketed states (where ©, © stand for the states 

I ©^ Tracks 1-3 I Tracks 4-6 aUt ^ I > ) Tracks 1-3 I ® ^ Tracks 4-6 r espeCti\>ely, and ■ Stands for 

any single-site state other than ©, ©). We denote by Sb r the subspace spanned by 
the bracketed states. 

1 There are two blank symbols because we will need different symbols to the left and right of 
the head, in order to enforce the constraint that there is only one head on the track. 
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We will often denote configurations of multiple tracks by writing them verti¬ 
cally, e.g. 


Pi 


<- 


© 

© 


Pi 

0 

1 


denote, respectively, a single-site configuration of two tracks, a configuration on 
two neighbouring sites of one track, and a configuration on two neighbouring sites 
of three tracks. Since we will restrict throughout to the subspace Sb r of bracketed 
states, the © and © states will only ever appear at the left and right ends of the 
chain, across all the tracks. As a shorthand, we will therefore denote configurations 
involving © and © with symbols stretching vertically across all tracks, e.g. 




We will also sometimes need to specify configurations in which the tracks are in 
any configuration other than ©, ©. We denote this with a -i© or -i© symbol 
stretching vertically across all tracks, e.g. 



© 


Pi 

1 



As usual in such constructions, the two-body Hamiltonian will contain two 
types of terms: penalty terms and transition rule terms. Penalty terms have the 
form \ab){ab\ where a, b are standard basis states. This adds a positive energy 
contribution to any configuration containing a to the left of b. We call ab an illegal 
pair, and denote a penalty term \ab){ab\ in the Hamiltonian by its corresponding 
illegal pair. We will sometimes also make use of single-site illegal states a. Note 
that single-site illegal states are easily implemented in terms of illegal pairs, by 
adding penalty terms \ax)(ax\ and \xa)(xa\ for all pairs ax and xa in which the 
single-site state appears. 

Definition 35 (Legal and illegal states) We call a standard basis state legal if it 
does not contain any illegal pairs, and illegal otherwise. 

We will make repeated use of Lemma 5.2 from [ 1109], which lets us use 
penalty terms to restrict to configurations matching a regular expression: 
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Lemma 36 (Regexps) For any regular expression over the standard basis of 
single-site states in which each state appears at most once, we can use penalty 
terms to ensure that any legal standard basis state for the system is a substring of 
a string in the regular set. 

All the regular expressions we will use start with a left-bracket state © and end 
with a right-bracket ©. So whenever we apply this Lemma, the only substrings of 
the regular set that are bracketed states are in fact complete strings in the regular 
set, not substrings. 

Transition rule terms have the form \(\f) - \(p)){{f\ - where \f), \ip) are 
states on the same pair of adjacent sites. We will often take \f) = \ab) and \ f) = \cd) 
to be standard basis states. This forces any zero-energy eigenstate with amplitude 
on a configuration containing ab to also have equal amplitude on the configuration 
in which that ab is replaced by cd. Conversely, if a zero-energy eigenstate has 
amplitude on a configuration containing cd, it must also have equal amplitude on 
the configuration in which that cd is replaced by ab. Thus we can think of these 
terms as implementing the transition rules ab —> cd, which we arbitrarily call 
the “forwards” direction, and the corresponding “backwards” transition cd —» ab. 
Following the notation in [G109], we will denote transition rule Hamiltonian terms 
by their associated forwards transitions ab —> cd or, more generally, \f) —» | ip). 

When a transition rule acts diagonally on a subset of the tracks, we will 
sometimes consider the restriction of the transition rule to those tracks. That is, for 
a transition rule with the general form \ab) T \f) r —» \cd) T \f) T c, where T is some 
subset of the tracks, the restriction of the rule to T is given by \ab) T —» \cd) T . 

When we specify a Hamiltonian term only on a subset of the tracks, we im¬ 
plicitly mean that it acts as identity on the remaining (unspecified) tracks. We will 
assume throughout this section that the ground state subspace of the Hamiltonian 
is restricted to the subspace 5 bl of bracketed states, with © at the left end of the 
chain and © at the right. (We show how to enforce this later, in Section 9.) 

7.2 Clock Oscillator 

Every clock - even one unusual enough to be encoded in superposition into 
the ground state of a quantum many-body Hamiltonian! - needs some form of 
oscillator that oscillates with a fixed period (such as a pendulum), and a counter 
that is incremented after each complete oscillation of the oscillator (such as the 
hands on a clock). For the counter, we will use a reversible counter Turing Machine, 
described in Section 7.4. This section describes the Track 1 clock oscillator. 

The legal configurations of Track 1 are states matching the regular expression 
©O*[©©(§)©©]©*©• By Lemma 36, we can enforce this using penalty 
terms. Independent of the configuration of any other track, a © arrow (either @ 
or @) sweeps right along the chain until it reaches the end, whereupon it turns 
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Figure 8: Evolution of the Track 1 clock oscillator. 


around and becomes a 0. (When a Q) reaches the end and turns around, the label 
on the arrow steps through the sequence (0), (1),..., ( K ). This is used later on to 
correctly initialise the other tracks. The restriction to L > K + 3 in Theorem 33 
ensures that there is enough space for such sequences of transitions.) This 0 
arrow sweeps left until it reaches the beginning of the chain, at which point it turns 
around and becomes a 0 again. We call this entire sequence an oscillator cycle. 
(One complete cycle is illustrated in Figure 8 for a chain of length six.) On a chain 
of length L, one complete oscillator cycle takes 2(L - 2) steps. The following 
transition rules on Track 1 enforce this: 









O0, 0© 


©, O 



When a left-moving 0 arrow returns to the beginning of the chain, it turns 
around and becomes a right-moving © arrow again. However, the label 0 on the 
arrow may change to a 1, depending on the Track 2 state. The @ arrow transitions 
to @ if the Track 2 state is p a (the initial state of the counter TM). Whereas 0 
always transitions to 0 unless the Track 2 state is p a . The following transition 
rules implement this: 



( 69 ) 
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where ->p a here denotes any Track 2 state other than p a . 


7.3 Initialisation sweep 

During the initial sweep of the Track 1 Q) from left to right and back, Track 1 
contains a @ or (^). We call this the initialisation sweep, and we want to use it to 
force Tracks 2 and 3 to be in the initial configurations 

©Po@*© and©#*©, (70) 


respectively. We do this by adding illegal pairs to forbid Track 2 from being 
anything other than © when a @ is over it (except at the beginning of the chain): 



and similarly to forbid Track 3 from being anything other than # when a @ is over 
it (except at the beginning of the chain): 



(72) 


where -i • again denotes anything other than the state • . 

At the beginning of the chain, we use illegal pairs to forbid Track 2 from being 
anything other than p a when @ is over it. We also forbid it from being p a when a 
© is at the beginning of the chain: 



(73) 


Thus any standard basis state containing a @ or (^) on Track 1 that matches 
the regular expressions (70), and is not in the initial configuration on Tracks 2 
and 3, will evolve (either forwards or backwards) under the oscillator transition 
rules into an illegal configuration, in at most L steps. 

If Tracks 2 and 3 are in their initial configurations, then a standard basis 
state containing a @ or @ will evolve backwards until it reaches the initial 
configuration |0 O ), the standard basis state for Tracks 1 to 3 which has the form: 



© 

© • 

•• © 


< 

Pa 

© •• 

• © 

© 


h 

# •• 

• # 



( 74 ) 
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7.4 Clock Counter 

7.4.1 Counter TM construction 

We construct the counter TM using a generalisation of the binary incrementing 
machine Inc from Lemma 27, which will increment integers written in base-£ 
instead of in binary. Constructing a base-f version of Inc is largely a straightforward 
extension of Lemma 27, though ensuring reversibility is slightly subtle. 

The following well-formed, normal-form, partial transition function imple¬ 
ments a reversible TM that increments any little-endian base-£ number written on 
the tape with a leading “start-of-tape” symbol I-, and returns the head to the starting 
cell without ever moving the head before the starting cell: 



h 

# 

0 

1 

i 

f-1 

Po 

(h P\,R) 





Pi 


(1 ,rf,R) 


(2, p 2 ,R) 

(i + l,P2,R) 

(0,Pi,R) 

P2 


(#,P3,L) 

(0, P3,L) 

( 1 5 P3 •> R) 

(i,P3,L) 

(£“ hP3,L) 

P'l 



(0 ,P' V L) 

(1 ,p' v L) 

( i,p' v L ) 

(£-l,p' 3 ,L) 

P2 


(#, p' 3 , L) 





P3 





(i,p 4 ,L) 


P3 




(1» P4t R) 



P4 

(K Pf, N) 

(0, Pa, L) 




Pf 


(#,Po,R) 

(0, po,R) 

(l,Po,R) 

(i, Po,R ) 

(£- hpo,R) 


Vie {2,...,£-2} 

(75) 

We now modify this machine in order to make it loop forever. We must be 
careful to ensure that we reenter the initial state reversibly. To do this, we introduce 
new initial and final states p a and p 0) , and modify the po and pf transitions such that 
the TM reversibly transitions back into the state po instead of halting, exploiting 
the fact that the symbol in the second tape cell will be blank in the first iteration, 
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and non-bla nk thereafter. The following set of transition rules accomplishes this: 



h 

# 

0 

1 

i 

<r-i 

Pa 

(\-,P- U R) 






P- 1 


(#,Po,L) 





Po 

(h Pi,R) 






P\ 


(1 ,P%,R) 

(i ,P' 2 ,R) 

(2,p 2 ,R) 

(J + 1, Pii R) 

(0 , Pl ,R) 

Pi 


(#, Pi,L) 

(0,p 3 ,L) 

(1 ,Pi,L) 

0, Pi,L) 

or-i ,pi,q 

Pi 



(0 ,p' 3 ,L) 

(1 ,P' 3 ,L) 

0, p' 3 ,l ) 

(t-hp' v L) 

Pi 


(#,P' 3 ,L) 





Pi 





0, Pa, L ) 

or-i ,pa,l) 

Pi 




(1,P4,L) 



Pa 

(h p f ,N) 


(0,P4,L) 




Pf 

(i -,p' f ,R) 






P'f 



(0, po, L) 

( l,Po,L ) 

(0 Po, L) 

(£-l,Po,L) 

Poj 

(1 -,Pa,R) 

(#, p a ,R) 

(0, Pa,R) 

(1 ,Pa,R) 

0 i,Pa,R ) 

O' - 1 ,Pa,R) 

V/ € {2, 2} 







(76) 


These transition rules implement a reversible base-£ counter TM. When started 
from the tape configuration consisting of a h symbol in the first cell followed 
by the all-blank tape (representing the number 0), this TM will loop indefinitely, 
incrementing the number written on the tape by 1 in each complete iteration. For 
brevity, we will refer to this particular initial tape configuration as the standard 
input to the counter TM tape. 

We now declare certain configurations of the counter TM to be “illegal”. We 
will choose these illegal configurations to be such that they are never entered 
by the counter TM started from the standard input. Later on, when we come to 
encode the counter TM in a local Hamiltonian, we will use penalty terms to give 
energy penalties to all the configurations we declare “illegal” here, so that the 
corresponding standard basis states of the spin chain are indeed illegal in the sense 
of Definition 35. For now, however, the “illegal configurations” simply define a 
particular subset I of the complete set of counter TM head and tape configurations. 

Anticipating the later use of penalty terms, when defining illegal configurations 
we will often make use of the illegal pair notation introduced above, where the top 
row denotes the counter TM head position and internal state p, and the bottom row 
denotes the tape symbols ab in the section of tape near the head: 

(77) 

A single row always denotes a section of tape: 


p 


a 

b 



p 

a 

b 


( 78 ) 
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Certain configurations that we declare to be illegal will specify that the counter 
TM head is or is not located at the starting cell. Again anticipating later use of 
penalty terms, we denote the starting cell by placing a © symbol on the left. We 
denote tape cells other than the starting cell with -■© on the left: 



(79) 


Note that the counter TM started from the standard input never moves its head 
before the starting cell, never reenters its initial state p a , and never enters its final 
state Pcj. We therefore declare any counter TM configuration which is in state p a 
with the head anywhere other than the starting cell to be illegal: 



(80) 


We also declare any configuration in state p 0J to be illegal: 


(81) 


Similarly, we declare any configuration in state p a with the head adjacent to a 
non-blank symbol to be illegal, i.e. any configuration containing the illegal pair 


Pa 





(82) 


The transition table in (76) is partial; some transitions are never used, so were 
not defined. For each (p, r) for which no transition rule is defined in (76), we 
declare configurations in which the counter TM is in state p and the head is reading 
r to be illegal, i.e. any configuration containing the illegal pair 


P 

T 


(83) 


For each configuration ( p L , r, L), (p R , r, R) or (p N , r, N) which is not entered by 
any transition rule in (76), we declare the corresponding TM head and tape config¬ 
urations to be illegal, i.e. any configurations containing one of the following illegal 
pairs: 



Pr 

T 



Pl 



T 


Pn 


T 



(84) 


Note that during the evolution of the counter TM started from a blank tape, the 
head is never more than one cell to the right of a non-blank tape symbol. In fact, 
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the only moment at which it is over a blank symbol at all is when the incrementer 
TM needs to carry a digit, and moves the head to the next blank cell in order to 
write the carry. We therefore declare configurations in which the head is more than 
one cell to the right of a non-blank tape symbol to be illegal, i.e. all configurations 
containing the illegal pair 


Note also that, when started from the standard input, the counter TM never 
modifies the i- in the first cell, never writes a i- anywhere other than the first cell, 
and never creates an embedded blank symbol on the tape. We therefore define all 
tape configurations without a i- in the starting cell, all tape configuration with a i- 
any where other than the first cell, and all tape configurations with an embedded 
blank, to be illegal. These illegal TM configurations correspond exactly to tape 
configurations matching the regular expression 

hi-#*, i e {0.. 1}. (86) 


The little-endian numbers written by the counter TM are never padded with 
leading 0’s, so the tape never has a 0 to the left of a #. We therefore define this 
combination to be illegal, i.e. all tape configurations containing the illegal pair 


0 # 


(87) 


Later on, we will only be interested in configurations of the first L cells of the 
tape, and will stop the counter TM just before it tries to apply a transition rule that 
maps out of this portion of tape. Therefore, a 0 will never appear in the L’th tape 
cell, and we declare configurations containing a 0 in this cell to be illegal: 


0 © 


( 88 ) 


Lemma 37 (Evolve-to-illegal) Any counter TM configuration that is reached 
starting from the standard input is legal. All other configurations with the head in 
position r e [0, L\ are either illegal, or evolve forwards or backwards to an illegal 
configuration within 0(L ) time steps in such a way that the head never leaves the 
portion of tape [0, L\. 


Proof The first part of the Lemma is true by construction, since the partial transi¬ 
tion rules of (76) implement the counter TM without using any of the undefined 
transitions, and the tape configurations defined to be illegal never occur when the 
counter TM is started from the standard input. 

Recall that the alphabet of the counter TM is S := {h, #} U {0,..., £ - 1}. Any 
tape configuration that does not start with a h, or which contains a h anywhere 
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other than the first cell, or which contains an embedded blank, or contains a 0 to 
the left of a #, or with a 0 in the L-th cell is illegal. So legal tape configurations 
match the regular expression 


h (-#)*-.[(),#]#*. (89) 

Configurations in which the head is more than one cell away from a non-blank 
symbol are illegal by (85), so all legal configurations have the head located in, or 
immediately adjacent to, the non-blank portion of the tape. All other configurations 
are illegal. 

We divide the legal configurations into separate cases, according to the internal 
state of the counter TM: 

State p a : All configurations in state p 0 are illegal by (80), except those with the 
head in the starting cell. The latter are illegal by (82) unless the tape cell to the right 
of the head is blank. The only such tape configuration which is legal is the standard 
input. Thus the only legal configuration in state p a is the initial configuration of 
the counter TM started from the standard input. 

State p-\: There are no transitions out of (p_ \,x) where x # in (76), so by 
(83) the only legal p-\ configurations have their head over a #. Moreover, there is 
only one transition into p-\ in (76). Hence, evolving any such configuration one 
step backwards according to (76) either enters a configuration that is illegal due to 
(83), or steps the head left and transitions to p a . But we have shown that the only 
legal p a configuration is the standard input. Thus any configuration in state p-\ 
not reachable starting from the standard input will evolve backwards to an illegal 
configuration in one time step. 

State p\: Every TM configuration of the form 


Pi 


h 

0 


0 

T 



t e E/ m , n > 0 

with the rest of the tape compatible with (89) is reached when evolving the counter 
TM from the standard input; this configuration is produced whilst the counter TM 
is incrementing the number with little-endian base-£ expansion 


(£-i)...(<r-i )r---. 

n 
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The remaining legal configurations have at least one tape symbol x to the left 
of the head that is neither 0, nor #, nor h: 


Pi 


h 

h 


lm 

X 

0 


0 

T 



x e S/ h, ■ ■ ■, im e S/{ h> #}, r e S/ [h) , m,n> 0. 

Evolving any such configuration backwards according to (76) either enters into an 
illegal configuration due to (83), or rewinds the head left replacing 0’s with £ — l’s, 
eventually reaching the configuration with the head immediately to the right of the 
x: 


Pi 


h 

h 



X 

£-1 


£-1 

r 



x e S/{ K # j0 ,i}, i\,...,i m g H/ {k#) , t g E/{ h> #}, m,n > 0. 

The state pi can only be entered by moving right. But the configuration (p\, x, R) 
is not entered by any transition rule in (76), so this configuration is illegal by (84). 

State p 2 : The state p 2 is entered when incrementing any digit other than 0, # 
or £ - 1, so every TM configuration of the form 


P2 


h 

0 


0 

X 

T 



n 

x € S/jh,#,o,i), t G H/ (1 _), n> 0. 


with the rest of the tape compatible with (89) is reached when evolving the counter 
TM from the standard input. This configuration is produced whilst the counter TM 
is incrementing the number with little-edian base-£ expansion 


(£-!)•••(£— l)( x - l) r ''' 

n 

The state p 2 can only be entered by moving right. Configurations ( p 2 ,x,R ) 
where x g {h, #, 0,1} are not entered by any transition rule in (76), so these 
configurations are illegal by (84). 
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The remaining legal configurations have at least one tape symbol y to the left 
of the head that is neither 0, #, nor k 


P 2 


h 

h 


Im 

y 

0 


0 

X 

T 



x £ E/ {h,#,o,i}» y e H/{H,#,o}i T e S/j,-}, z'i,..., i m £ E/{ h; #}, m,n> 0. 

Evolving any such configuration backwards one time step according to (76) either 
produces an illegal configuration due to (83), or produces the configuration: 


P i 


h 

h 


Im 

y 

0 


0 

x - l 

T 



x £ E/{i- ,#,o,i}, y € E/ { h ,#,o}, t e H/( h ), z'i ,..., i m £ E/{ h ,#}, m,n> 0. 

But we have already shown above that all such pi configurations evolve backwards 
in time to an illegal configuration. 

States p' 2 , p”: The argument for states p' 2 and p" is very similar to that for p 2 , 
except that the legal value of x is now 1. 

State pp. The state p 3 is entered by stepping left from p 2 , so every TM configu¬ 
ration of the form 


Pi 



0 


0 

T 



t € H/{ hj # j0 ,i}, n > 0. 


with the rest of the tape compatible with (89) is reached when evolving the counter 
TM from the standard input. This configuration is produced whilst the counter TM 
is incrementing the number with little-edian base-£ expansion 


(£—!)•••(£— 1)( T - 1) • • • 


There is no transition out of (p 3 ,x) where r e {h,#, 0,1} in (76), so all such 
configurations are illegal by (83). The remaining legal configurations have at least 
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one tape symbol x to the left of the head that is neither 0, #, nor k 


P3 


I- 

h 



X 

0 


0 

T 



T € € S/{i-,#,o} h, ■ ■ ■, i m e S/{ K #j, m,n> 0. 

Evolving any such configuration forwards one time step according to (76) ei¬ 
ther produces an illegal configuration (for example if n = 0), or produces the 
configuration: 


Pa 


h 

h 



X 

0 


0 

T 



T € S/{ h; # ; o,i},x e H/( hj # ; o} z'i,...,z w e E/ {h>#) , m,n> 0. 

We show below that any such p 4 configuration evolves to an illegal configuration. 

Stoic The argument for state p' 3 is very similar to that for pi, except that 
r = 1 in this case. 

State p 4 : In a legal state p 4 configuration, the head must be over a 0 or a h, since 
there is no transition out of (p 4 , v) for x & {0,1-}. We first analyse the case in which 
it is over a 0. 

Any TM configuration of the form 



n > 1 

with the rest of the tape compatible with (89) is reached by evolving the counter 
TM from the standard input. 

The remaining legal configurations have at least one tape symbol x to the left 
of the head that is neither 0, #, nor K 


Pa 


h 

h 



X 

0 


0 



x e E/( hj#j0 ), z'i,...,z m g H/ |hj#) , m>0,n> 1 
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Evolving any such configuration according to (76) steps the head left over the 0’s 
until the head is over the x: 


P\ 


I- 

h 


hn 

X 

0 


0 



x € H/ (k#! o}, h ,..., i m e S/{ K #}, m >0,n> 1. 

But there is no transition out of ( p 4 , x ) in (76), so this configuration is illegal by 
(83). 

We now analyse the case in which p 4 is over K Evolving one step backwards 
according to (76) steps the head right and either enters an illegal configuration, 
transitions to state p 4 with the had over a 0, or transitions to states p 3 or p' y In all 
three cases, we have already shown that all such configurations not reachable from 
the standard input evolve to illegal configurations. 

State pf.‘ There are no transitions out of (pj, x) where x ^i- in (76), so the only 
legal pf configurations have the head over a K The latter can only appear in the 
starting cell in legal configurations. Evolving any such configuration backwards 
according to (76) for one time step transitions to the state p 4 , and we have already 
shown than all p 4 configurations not reachable from the standard input evolve to 
illegal configurations. 

State p'p There are no transitions out of (p'j, x) where x e {i-, #} in (76), so the 
only legal p^ configurations have the head over a symbol in the set {0,..., £ - 1}. 
Evolving any such configuration backwards according to (76) for one time step 
transitions to the state pf, and we have already shown than all p/ configurations 
not reachable from the standard input evolve to illegal configurations. 

State /?o: There are no transitions out of (p 0 , *) where x in (76), so the only 
legal po configurations have the head over a K The latter can only appear at the 
beginning of the tape in legal configurations. Since the only transitions in (76) 
into po are from p\ or p'j, evolving any such configuration one step backwards 
according to (76) either enters an illegal configuration, or steps the head right 
into state p-\ or p' f . In both cases we have already shown that configurations not 
reachable from the standard input evolve to illegal ones. 

State p, 0 : The counter TM never enters the state p 0J , and all such configurations 
are illegal. □ 
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7.4.2 Counter TM Hamiltonian 


We want all legal Track 2 configurations to contain a single state p e P' marking 
the location of the counter TM head, and blanks everywhere else. I.e. the Track 2 
configuration should match the regular expression ©( *p §)'(> where p e F. 
By Lemma 36, we can enforce this using penalty terms. This Track 2 regular 
expression allows us to identify Tracks 2 and 3 with the configurations of the 
counter TM (with internal state p, its position being the head location, and Track 3 
the tape). This will be explicitly or implicitly exploited in all that follows. 

To encode the evolution of this counter TM in our local Hamiltonian, we encode 
its transition rules in transition terms. However, the transition rules will only apply 
when a @ arrow on Track 1 sweeps past the counter TM head encoded on Track 2, 
so that the counter TM advances exactly one step for each complete left-to-right 
sweep of the @ arrow. Since we are restricting to two-body interactions and the 
arrow is moving to the right, transitions in which the head moves left must be 
triggered when the @ is to the left of the TM head, in order to both move the arrow 
and update the head location. This way of implementing the left transitions means 


that 



P 


cr 


(90) 


never occurs during the evolution of the counter TM if 6(p, cr) = (r, p L , L). We add 
a penalty term to make such configurations illegal. 

Transitions in which the head moves to the right or stays still will be triggered 
when the @ is on top of the head. In order to avoid two transitions in which the 
head moves right from being triggered during the same sweep of the @, we must 
split these latter transitions into two stages, in which we first transition into an 
auxiliary state p' R e P' R , and then transition into the correct state p R e P R in the 
following time step (adapting the construction of [GI09]). Thus, when the @ is on 
top of the the head and the head will move right, we will update the tape, move the 
head, and step the @ to the right. But we transition into the auxiliary state p' R e P' R 
instead of p R e P R . In the next step, with the @ now on top of the p' R state, we 
transition to the correct state p R , and step the @ to the right once more so that it is 
no longer on top of the head. This way of carrying out the right-moving transitions 
means that p' R only ever appears below a @, so we add a penalty term to make all 
other p' R configurations illegal: 



(91) 


For TM transition rules 5(p, cr) = (r, p N , N ), 6(p, cr) = (r, p L , L ) and 5(p, cr) = 
(T,p R ,R) (where p e P\ {Pa] , p L ,\>j< e Pl,n,r and cr,r e S), the following local 
transition terms on Tracks 1 to 3 implement the desired transformations (cells 
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can be in any state, and are left unchanged by the transition): 
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© 
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P'r 

© 

-> 

Pr 


9 

Pr 

© 

-> 

Pr 
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(92a) 


(92b) 


(92c) 


Transitions in which the head moves left from the end of the chain are already 
covered by the left-moving transition in (92b). Transitions in which the head would 
move right ofF the end of the chain are not implemented. (Transitions in which the 
head would move left off the beginning of the chain are not implemented either, 
but in fact the counter TM construction of Section 7.4 will never attempt such a 
transition when initialised properly. Moreover, we just declared such configurations 
illegal in (90).) 

When started from the blank input (which represents the number 0), the encoded 
counter TM will loop g L ~ 3 times before it exceeds the L — 2 tape space available 
on Track 2 (the space overhead of 1 is due to the way the incrementer TM in 
implemented in (75)). The incrementer machine takes 0(log x) steps to increment 
the number x (written on the tape in base-4). Thus the counter TM will run for at 
least £l(£ L ) (and at most 0(log L^ L )) time-steps before it exceeds the available tape 
space. 

The transition function for the counter TM defines a transition rule for each pair 
(p, r) e Px S. So for any legal configuration containing a p € P on Track 2, exactly 
one of the transition rules in (92a) and (92b) will apply during each © sweep. 
Similarly, for any configuration containing a p' R e P' R , a (92c) rule will apply. It 
is easy to verify that if any of the transition rules from (92a)-(92c) is applied to 
any legal configuration, none of the (92a) and (92b) rules can apply again to the 
resulting configuration. If the final (92b) rule applies, the (92c) rule will apply 
exactly once in the following step, after which none of the rules can apply again to 
the resulting configuration. The transition rules in (92a)-(92c) therefore implement 
a single step of the TM during each left-to-right sweep of the @ state on Track 1, 
as required. We call this phase of the evolution, in which Track 1 contains a © or 
©, the computation phase. 
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We also modify the clock oscillator rules from Section 7.2 involving a right- 
moving © arrow to only apply when there is no TM head under the ©, since the 
© movement is already taken care of by the above transition rules when there is a 
TM head present: 
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We must also enforce illegality of all the configurations declared to be illegal in 
Section 7.4.1. To do this, we simply introduce local penalty terms corresponding 
to the combinations declared illegal in (80)-(85), (87) and (88), and enforce the 
Track 3 regular expression of (86) using penalty terms and Lemma 36. 

The complete set of Track 1 to 3 transition rules is summarised in Table 1 . The 
regular expressions on Tracks 1 to 3 enforced by penalty terms, together with all 
the additional illegal pairs defined so far, are summarised in Table 2. 

Table 1: All transition rules for Tracks 1 to 3. 


©,© rules 


Track 1 rules 
© rules 


© rules 
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transition rules marked with an (*) will be replaced in Section 7.6 by a set of rules that also 
act on Tracks 4, 5 and 6. 
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Track 1, 2 and 3 rules 
® rules 


© rules 
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Table 2: All illegal pairs and regular expressions enforced by 
illegal pairs for Tracks 1 to 3. 


Regular expressions 

Track 1 Track 2 Track 3 
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Illegal pairs 
Tracks 1 and 2 



Illegal pairs 

Tracks 1 and 3 Track 2 Tracks 2 and 3 Track 3 
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Illegal pairs 

Tracks 2 and 3 for undefined transitions in (76) 
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Illegal pairs 

Tracks 1, 2 and 3 if 5{p, cr ) = (p L , r, L ) (90) 
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cr 


7.5 Clock Hamiltonian 

Before considering the remaining tracks, it is helpful to analyse in depth the 
Hamiltonian defined so far on the first three tracks. For that we start with the 
definition of well-formed states. 

Definition 38 (Well-formed state) We say that a standard basis state on Tracks 1 
to 3 is well-formed if it is a bracketed state, its Track 1 configuration matches 
the regular expression ©O* [©(§)©©©] @*@, and its Track 2 configuration 
matches the regular expression © 0*p where p e P'. 

The penalty terms defined previously give an energy penalty of at least 1 to any 
standard basis state in S br that does not match the desired regular expressions on 
Tracks 1-3. Thus only well-formed states can have zero energy. The following 
result shows that only one transition rule can apply to each well-formed state. 

Lemma 39 (Well-formed transitions) For any well-formed standard basis state, 
at most one transition rule applies in the forward direction, and at most one in the 
backwards direction. Furthermore, the set of well-formed states is closed under 
the transition rules. 

Proof A well-formed standard basis state contains exactly one of Q), @ 

or © on Track 1. Thus clearly only transition rules from one of the sets in Table 1 
(the set of (^) rules, @ rules, or © rules) can apply in the forward direction. 
The left hand sides of the rules within the ©, © and © sets are manifestly 
mutually exclusive. Since the counter TM is deterministic, there is a unique TM 
transition that applies at any step, so the left-hand-sides of the © rules are also 
mutually exclusive. 

The counter TM is reversible, so there is also a unique backwards TM transition 
at every step. Thus the same argument applies in the backwards direction to the 
right hand sides of the rules in each set, with the exception of the @ turning rule 
that changes a © into a © (in the backwards direction). But this rule clearly 
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cannot apply at the same time as the right hand side of any rule from the © set, 
which concludes the proof of the first part of the lemma. 

It is straightforward to verify that all the transition rules in Table 1 preserve 
well-formedness, implying the second part. □ 

Recall that the state \<po) is the standard basis state with ©©<§)*© on Track 1, 
©Pa ®*@ on Track 2, and © h #* © on Track 3. This corresponds to having 
the clock oscillator in the first state of the Track 1 sequence (see Figure 8), and the 
counter TM Tracks 1 and 2 in their initial configurations. There is no backwards 
transition out of \<Pq), so we will refer to it as the initial clock state. Let \<p t ) be the 
standard basis state obtained by applying t transitions (in the forwards direction) to 
\<po), which is well-defined thanks to Lemma 39. 

Let us consider the form of the states \<p t ). Note that | <p 0 ) contains a Q) on 
Track 1. The only transition rules that apply when Track 1 contains a @ or @ 
(see Table 1) are those that sweep the arrow all the way to the right and back again, 
without affecting Tracks 2 or 3. We saw in Section 7.2 that one complete oscillator 
cycle takes 2 (L - 2) steps. Thus t < 2(L - 2) corresponds to the initialisation sweep, 
in which the state | <p t ) contains a Q), @ or (§) on Track 1, and Tracks 2 and 3 
remain in their initial configurations. 

At t = 2 (L - 2) + 1, the turns around and (because Track 2 contains a p 0 
at the beginning of the chain) becomes a @. This @ is now over a p a on Track 2. 
Thus the next step implements the first step of the counter TM on Tracks 2 and 3, 
causing Track 2 to transition to p a . 

Recall from Section 7.4 that the transition rules implement one step of the 
counter TM during each left-to-right sweep of the @ on Track 1, and then move 
the © back to the beginning of the chain in the second half of the oscillator cycle. 
Note that the counter TM construction of Section 7.4 never transitions back to its 
initial state p a . So when the © reaches the beginning of the chain, it will never find 
a p a on Track 3. Thus it will always transition into ©, never into @ (see Table 1). 
Therefore, for all t > 2(L - 2) the state \<p t ) contains a @ or © on Track 1. At 
t = 2 (L - 2)(n + 1) the @ in \<p t ) is back at the beginning of the chain, and Tracks 2 
and 3 are in the configuration corresponding to the n’th step of the counter TM. 

Eventually, the number written on Track 3 will be incremented until it reaches 
4 /-_4 , the maximum number that can be written in base-4 in L - 4 digits (the -4 
accounts for the space-overhead of 2 required by the Inc TM implementation, plus 
the two bracket states at the ends of the chain). When the counter TM tries to 
increment this number, the head will move right until it reaches the cell adjacent 
to the © at the end of the chain. At this point, it will be in an internal state that 
would normally step right in order to write the carry. But the site to the right of the 
head contains © instead of a blank Tape 2 cell. When the @ on Track 1 reaches 
the head at the end of on Track 2, there is no further forwards transition out of 
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this configuration (see Table 1). Let T be the number of transitions required to 
reach this configuration \<p T ) starting from the initial state \4>o). We will refer to 
\(j) T ) as the final clock state. Since this configuration occurs when the counter TM 
is about to exceed L tape space, which occurs after at most <9(log Lf L ) steps of 
the TM (see Section 7.4) each of which requires 2(L - 2) transitions, we have that 
T = 0(L\ogL^ L ). (Though this bound on the run-time will not be important for 
our purposes.) 

The following result will be important later when we come to analyse the 
Hamiltonian. 

Lemma 40 (Evolve-to-illegal) Evolving any \fi t ) forwards or backwards in time 
according to the transition rules will never reach an illegal configuration. All other 
well-formed standard basis states will evolve either forwards or backwards to an 
illegal configuration after 0(L 2 ) transitions. 

Proof As in Lemma 37 the first part of the Lemma is true by construction. For the 
second part, let us analyse all possible well-formed standard basis states, which we 
divide into two cases depending on the type of arrow state we have on Track 1: 

Case 1: Track 1 contains (J) or (§). If Tracks 2 and 3 are in the initial 
configuration, then we are in one of the \4>,) states with 0 < t < 2(L - 2). 

If Tracks 2 and 3 are not in the initial configuration, then within at most 0{L) 
steps the @ or (^) will move forwards or backwards along the chain due to the 
initialisation sweep clock oscillator transition rules, until it is over a site containing 
the wrong initial Track 2 or Track 3 state. But this is illegal due to the Track 1-3 
illegal pairs from Table 2. 

Case 2: Track 1 contains @ or ©. Recall that there is a one-to-one correspon¬ 
dence between configurations of the counter TM and well-formed standard basis 
states that do not contain a p' R e P' R on Track 2. 

Well-formed configurations containing a p' R e P' R on Track 2 evolve by a single 
backwards transition from (92b) to a configuration containing a p i P' R on Track 2, 
and evolve by a single forwards transition from (92c) to a configuration containing 
a p R e P R . A configuration containing p' R e P' R is therefore reachable (resp. 
unreachable) by transition rules from the initial configuration iff the corresponding 
counter TM configuration with p' R replaced by the corresponding p R is reachable 
(resp. unreachable). 

For any well-formed configuration, each sweep of @, the transition rules im¬ 
plements exactly one step of the counter TM encoded in Tracks 2 and 3. The only 
counter TM transitions that are not implemented are those that would move the 
head beyond the L sites between the © and ©. Therefore, if the configuration 
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of Tracks 2 and 3 is not reachable from the initial configuration, then Lemma 37 
guarantees that evolving it forwards or backwards by the transition rules will reach 
an illegal configuration in 0(L ) steps of the counter TM, which takes 0(L 2 ) transi¬ 
tions. (Note that it is crucial for this argument that to reach an illegal configuration 
in Lemma 37 never requires moving the head outside the portion of tape [0, L].) 

It only remains to consider well-formed states with Track 2 and 3 configura¬ 
tions that are reachable from the initial configuration. For the initial Tracks 2-3 
configuration itself, if © is in the starting cell then we have the state \<p t ) for 
t = 2(L - 2) + 1. If not, either we have © to the right of the head or © on Track 1, 
with p a at the start of Track 2. In both cases, evolving such a state forwards in time 
will reach a state with a © over the p a within O(L) transitions, which is illegal by 
the second Track 1 and 2 illegal pair from Table 2. 

If the Track 2-3 configuration is not the initial one, then all positions of @ and 
© give rise to elements in the set {| <p t )} t except exactly those declared illegal in 
(90) and (91). 

This completes the proof of the Lemma. □ 

7.6 QTM Hamiltonian 

With the clock in place, encoding an arbitrary quantum Turing Machine in the 
Hamiltonian is now possible using ideas from [GI09] (and is similar to the construc¬ 
tion used in Section 7.4 to encode the reversible counter TM). The configurations 
are now quantum states, so there can be arbitrary superpositions over standard 
basis states. But it suffices to verify that the construction does the right thing on 
standard basis states; this then extends to arbitrary superpositions by linearity and 
well-formedness. 

When we later come to analyse the spectrum, we will require our Hamiltonian 
to be standard-form in the following sense: 

Definition 41 (Standard-form Hamiltonian) We say that a Hamiltonian H = 
H tm ns + H pen acting on a Hilbert space ( H = (<D C ® C e )® L = (C c )® L ® (C~)* L =: 
‘He ® ‘H q is of standard form ifH tranSiPen = X,ti' h^^, and h tmns>pen satisfy the 
following conditions: 

(i). h tmns e ( B ((C c ® <D^)® 2 ) is a sum of transition rule terms, where all the 
transition rules act diagonally on C c <g> C c in the following sense. Given 
standard basis states a, b,c,d e C c , exactly one of the following holds: 

• there is no transition from ab to cd at all; or 

• a,b,c,d e C c and there exists a unitary U abcd acting on C e <g> C e 
together with an orthonormal basis (I ^‘abed^i f or ^ ® C 2 , both de¬ 
pending only on a, b, c, d, such that the transition rules from ab to cd 
appearing in h tmns are exactly \ ab) W abcd ) -» I cd) U abcd W abcd ) for all i. 
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(ii). h pen e ( B ((C c <S> C e )® 2 ) is a sum of penalty terms which act non-trivially 
only on (C c ) 02 and are diagonal in the standard basis. 

In our case, <D C corresponds to Tracks 1-3 and <D e to Tracks 4-6. Transitions 
from ab to cd will correspond exactly to the transitions shown for Tracks 1-3 
in Table 1. While adding new transitions involving Tracks 4-6, we will need 
to remove some of the transitions in Tracks 1-3, but they will be recovered as 
restrictions of the new rules to those tracks. (In Table 1, transitions that will 
be replaced later are marked with an asterisk.) We will define the unitaries U a bcd 
partially, and then complete them to a full unitary. All we have to take care of is that 
(i) the new transition involving the rest of the tracks, when restricted to Tracks 1-3, 
recover exactly the transitions removed from Table 1 ; and that (ii) orthogonality is 
preserved in the construction. 

By the standard form, 

/*trans = £ (I ab) - \cd))((ab\ - ( cd\ ) ® (21 - U ahcd - U^J . (94) 

ab^cd 

Now, if we start with a family of QTM P n which satisfies part (iii) of Theorem 11, 
it will be immediate from our construction that the partial definition of all U abcd 
will only involve elements in the set 

5 = jo, 1, ±-^=, e m f e in7rW |. (95) 

As argued above, the same will hold for the full U abcd , which gives part (vi) in 
Theorem 33. 

7.6.1 QTM transition rules 

We first analyse how to incorporate the QTM transition rules into the Hamiltonian. 
We will simulate the QTM transition rules 6(p, cr) = Y^ T , q ,n d( p, cr, r, q, D) |r, q, D) 
during the second half of the clock oscillator cycle, when the Track 1 © state 
sweeps from right to left. This is done by including the following transition rule 
terms (cf. (92a)-(92c)), where q is any state in Q/{q/} (i.e. we exclude all the 
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transitions out of the QTM’s final state): 
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q 


cr 



^ 6{q, o'. 


T.qx 


T,qL 



(96a) 


(96b) 


As in the implementation of the reversible counter TM (Section 7.4), these transition 
rules implement the Reft-moving, Non-moving and Light-moving transitions 
separately. This time, since the © arrow is sweeping to the left, it is the left- 
moving transitions that are implemented in two stages, first transitioning to an 
auxiliary state q' L € Q' L , then in a second step transitioning to the corresponding q L 
state. 

Note that if the QTM head ever ends up at the end of the chain, the transition 
rules of (96) can never apply, because they require the © to be to the right of the 
QTM head. Thus the final Track 5 tape cell is not used, and the effective tape 
length available to the QTM for a chain of length L + 3 is only L rather than 1 L + 1. 
If the Track 4 QTM head ends up next to the © at the very end of the chain, this 
indicates that the head has stepped off the usable portion of the tape, regardless of 
the internal state. 

We also include the following transition rules involving a left-moving © arrow 
to cover the case where there is no QTM head present: 
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7.6.2 QTM Unitarity 

With the set of transition rule terms defined so far in (96), there is no transition 
out of the final state qf of the QTM. There is also no transition out of a QTM 
configuration in which the QTM head is at the end of the chain, and the next step 
would move the head to the right, off the end of the chain. The same is true of 

1 Two cells are always used by the ©, © markers. 
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configurations trying to move the head left when it is located at the beginning of 
the chain. 1 

Ideally, we would like it if, when the QTM reaches one of these configurations, 
the clock continued ticking but the encoded QTM did nothing for the remaining 
time steps. However, simply omitting transitions out of these configurations causes 
a subtle but important issue relating to unitarity of the quantum transition rules, 
which will be crucial to the analysis of the resulting Hamiltonian in Section 7.7. 


If we complete the partial U (associated with the Track 1-3 transitions O© 


©(§)) as defined so far by (96)-(98) to a unitary, this will add additional transitions 


out of the final state qf, and also out of configurations with the head at the beginning 
or end of the chain. But this means that, instead of the encoded QTM doing nothing 
for the remaining time steps after reaching such a configuration, it will continue to 
evolve in some arbitrary way, depending on the particular choice of completion 
of U. On the other hand, if we omit transitions out of these configurations, the 
resulting U will only be a partial isometry, not a full unitary. 

We need to use up any remaining time steps in a controlled way, such that the 
Track 5 QTM tape is left unaltered after the QTM enters q/ or runs out of tape. 
To this end, we deliberately add additional transitions out of the final state qf, and 
also out of configurations with the head at the beginning or end of the chain, which 
cause the QTM to switch over to running some other, inconsequential computation 
to use up the remaining time. This time-wasting computation will use Track 6 as 
its tape, leaving the configuration of the Track 5 QTM tape untouched. 

The only requirement on this time-wasting computation is that it must be 
guaranteed not to halt or run out of tape before it has used up all the remaining 
time steps of the clock (otherwise we face the same unitarity issue once again). 
The natural choice is simply to run the same base-f counter computation as the 
clock (but running on Tracks 4 and 6). The counter TM never halts; it increments 
the number written on its tape forever, or - when encoded on a finite chain - until 
it runs out of tape space. Furthermore, since the clock will already have “ticked” 
for some positive number of time steps before this time-wasting counter TM starts 
running, the clock counter TM is guaranteed to run out of tape (and hence the clock 
stop ticking) first. 

Transitioning from qf to the time-wasting counter TM can be accomplished 
simply by dovetailing the counter TM after the QTM, and encoding this dovetailed 
machine instead of the original QTM. Concretely, this means we must add to 
the Hamiltonian more transition rule terms of the form (96) that encode all the 
transition rules of the counter TM from (76). 2 Note that Track 6 contains additional 


Although the QTMs we will encode will in fact never do this. 

2 Since these are the transition rules of a classical reversible TM, they have a particularly simple 
form, the right hand side of each (96a) rule containing a single term. 
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subscripted variants \- q of the h tape symbol; the time-wasting counter TM transition 
rules simply ignore the subscript, and treat all of these as a I-. (I.e. the counter 
TM rules that read a i- symbol in (76) are duplicated for each \- q variant.) These 
time-wasting counter TM transition rules involve a disjoint set of Track 4 internal 
states q e P' to those q e Q used by the QTM. Thus the left- and right-hand 
sides of all time-wasting counter TM (96) transition rules are orthogonal to those 
encoding the QTM transition rules. 

We also add the following transition out of the final state qf of the original 
QTM into the initial state p a of the counter TM, which dovetails the counter TM 
after the QTM. This transition rule acts on Tracks 1, 4 and 6: 



(98) 


We must also add transitions that switch to running the time-wasting counter 
TM if the QTM runs out of tape, i.e. if it enters a configuration in which the head 
is at the beginning of the tape and would move left in the next time step, or one 
with the head at the end of the chain that would move right. (Though in the case of 
the proper machines considered in parts (iv) to (vi) of Theorem 33, the former will 
never occur.) Moving left off the beginning of the tape is particularly easy, both 
because the head is already at the beginning of the chain, and because left-moving 
QTM transitions are implemented by first transitioning to an auxiliary q' L e Q' L 
state without moving the head. The following transition rule acting on Tracks 1, 2, 
4 and 6 accomplishes this: 
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Note that, to preserve reversibility (unitarity), we must keep a record of which 
internal state q' L led to the QTM running out of tape. We record this in the \- q 
symbol written to the Track 6 time-wasting TM tape. 

Configurations in which the head step right off the end of the chain involve 
slightly more effort, as the head is at the wrong end of the chain to start running the 
time-wasting counter TM. From Section 7.6, a configuration in which the Track 4 
QTM head is at the very end of the chain indicates that the head has stepped 
right off the usable portion of the Track 5 tape. We first transition from any such 
configuration into an auxiliary Track 4 head reset state r, which moves the head 
all the way back to the beginning of the chain, before transitioning into the initial 
state of the counter TM. This requires three additional sets of transition rules. 
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We want the first new rule to act on Tracks 1 and 4 at the very end of the chain, 
and transition into the auxiliary reset state r: 



( 100 ) 


However, we have to make it compatible with the transitions already defined for 
Tracks 1-3. Hence, instead of (100), we include the following transitions acting on 
Tracks 1,2,4 and 1,2,3,4. 




( 101 ) 


The second new rule also acts on Tracks 1 and 4, and steps the auxiliary r state 
to the left along with the ©: 
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o 

r q 



( 102 ) 


The third acts on on Tracks 1,2 4 and 6 at the very beginning of the chain, and 
transitions into the initial configuration of the time-wasting counter TM: 



(103) 


Again, to preserve reversibility, we record in the \- q symbol which internal state 
caused the QTM to run out of tape. 
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Table 3: All transition rules for Tracks 4, 5 and 6. 


Track 1 and 4 rules 
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Track 1, 4 and 5 rules 
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Track 1, 4 and 6 rules 




Track 1, 2 and 4 rules 
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Track 1, 2, 4 and 6 rules 



Track 1, 2, 3 and 4 rules 
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The transition rules in Table 3 implement all the transitions that will take place 
for a properly initialised QTM. All of them have the desired form \ab) |i j) —> 

| cd) U a bcd |z j) where ab —> cd is a Track 1-3 transition and U a bcd are partial isome¬ 
tries defined by their action on a subset |i j) of the computational basis of Tracks 
4-6. We can now complete the unitaries U abc d at will. All we have to check is that 
they do indeed define a partial isometry. 

The left hand sides of each transition rule in (96)-(99) and (101)—(103) are 
manifestly orthogonal. The right hand sides of all (96b), (97)-(99) and (101)—(103) 
rules are clearly mutually orthogonal and normalised, and are also orthogonal to 
the right hand sides of all (96a) rules. 

It remains to show that (96a) defines an isometry. Recall that there are two 
sets of (96a) terms, one encoding the QTM transition rules, and the other the 
time-was ting counter TM. Furthermore, the left- and right-hand sides of all rules 
in the QTM set are orthogonal to those of the counter TM set, as they involve 
disjoint sets Q' and P' (resp.) of internal states. Thus it suffices to prove that each 
set of transition rules considered separately defines an isometry. Since classical 
reversible TMs are special cases of QTMs, we prove the result for QTMs. 

Let \\f/(p, cr)) denote the state on the RHS of (96a). Then 

\\i//(p,o-)\\ = ^j\6(p,cr,T,q R ,R)\ 2 + ^\6(p,cr,T,q N ,N)\ 2 

r,q R r,q N 

+ 5{p, cr, r, q L , L) | 2 (104) 

r,qL 

= ||<5(/?, cr)|| = 1 

by the normalisation condition of Theorem 19, so the right hand sides of (96a) are 
normalised. 

Similarly, for (puT t ) ± (p 2 , r 2 ), 

(i/s(.PuO-i)\il/(p2,o- 2 )) = ^ j S(pi,o- 1 ,T,q R ,R)*6(p2,o-2,T,q R ,R) 

r,qR 

+ V 5(pi ,<t u t, q N , N)*5(p 2 , ct 2 , r, q N , N ) 
rm (105) 

+ ^ 6(p\ ,cti,t, q L , L)*6(p 2 , cr 2 , r, q L , L ) 

r,q L 

= (6(pi,T l ),6(p 2 ,T 2 )) = 0 

by the orthogonality condition of Theorem 19, so the right hand sides of (96a) are 
mutually orthogonal. Thus the QTM and time-wasting counter TM transition rules 
from (96a) preserve orthonormality, as required. 
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7.6.3 QTM initialisation sweep 

We use penalty terms that only apply during the initialisation sweep to enforce 
the correct initial configurations of Tracks 4, 5 and 6. For Track 4, we use the 
following penalty terms, which force a single q 0 at the very beginning of the track: 



(106) 


The (j^) Track 1 states only ever occur during the initialisation sweep in the final 
K sites at the right end of the chain. We use these to force the initial configuration 
of Track 5 to contain K blank symbols at the very end, as required by Theorem 33. 1 



(107) 


Track 6 is forced to be in the all-blank configuration, except for the very first 
cell which contains a K 



(108) 


Table 4: Initialisation illegal pairs for Tracks 4-6. 


Illegal pairs 

Tracks 1 and 4 Tracks 1 and 5 Tracks 1 and 6 
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7.7 Analysis 

Lemmas 39 and 40 are the key results needed to prove the desired ground state prop¬ 
erties of the Hamiltonian we have constructed, thanks to the “Clairvoyance Lemma” 

1 These additional blank symbols will be needed later to provide for the small space overhead 
of the QTM constructed in Theorem 11 . 
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from Aharonov et al. [Aha+09, Lemma 4.2] (see also [GI09, Lemma 5.6]). 1 Since 
the Clairvoyance Lemma has previously only been stated and proven for specific 
constructions, we state and prove a very general version of the Lemma here, which 
applies to any standard-form Hamiltonian, as defined above: 

Lemma 42 (Invariant subspaces) Let H tmns and H pen define a standard-form 
Hamiltonian as defined in Definition 41. Let S = {S/} be a partition of the 
standard basis states of He into minimal subsets S, that are closed under the 
transition rules (where a transition rule \ab) CD \ fi) —» \cd) CD U a b c d \ifr) acts on He 
by restriction to (C c )® 2 , i.e. it acts as ab —» cd). 

Then H = ((J) s 7C,.) ® Hq decomposes into invariant subspaces ( K Si ® Hq of 
H = H pen + H tram where ( K Si is spanned by S 

Proof H pen is diagonal in the standard basis by definition (part (ii) of Definition 41), 
so ( Ks, ® Hq are trivially invariant under H pen . But, by the form of the transition 
rule terms (part (i) of Definition 41), the image H tmns \x) c \ip) Q of a standard basis 
state |x) c under H Udns has support only on standard basis states of He that are 
reachable by transition rules from \x) c . Closure of ( K$ i ® Hq under H lrMK is then 
immediate from the definition of S ,. □ 

Lemma 43 (Clairvoyance Lemma) Let H = H tmns + H pen be a standard-form 
Hamiltonian as defined in Definition 41, and let *7 Q be defined as in Lemma 42. 
Let AffiKs) denote the minimum eigenvalue of the restriction H\x s ? jHo °f H = 
H t ram + H pen to the invariant subspace H s <g> Hq. 

Assume that there exists a subset *W of standard basis states for He with the 
following properties: 

(i) . All legal standard basis states for He are contained in AV. 

(ii) . AV is closed with respect to the transition rules. 

(iii) . At most one transition rule applies in each direction to any state in AV. 

(iv) . For any subset S c AV that contains only legal states, there exists at least 

one state to which no backwards transition applies. 

Then each 7C, falls into one of the following categories: 

(1) . S contains only illegal states, and A(fiKs ) > 1. 

(2) . S contains both legal and illegal states, and AffiKs ) = Q(1 /\S | 3 ). 

'Note that our definition of “legal” states follows that of [GI09], which differs from the 
definition in [Aha+09]. 
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(3). S contains only legal states, and do (70) = 0. The corresponding eigenspace 
is 



H, 


trans 



= span 


1 151-1 ) 

VCT g \>)c l^t> e | 


(109) 


where \t) c are the states in S, |<^ 0 ) is any state in Hq, and \if/ t ) := U t ...U\ \^o)q 
where U, is the unitary on Hq appearing in the transition rule that takes 
1 1 — l) c to | t) c . All other states in TO ® Hq have energy at least Q(l/|51 2 ). 


To prove this, we will need Kitaev’s geometrical lemma: 

Lemma 44 (Geometrical Lemma - Lemma 14.4 in [KSV02]) Let A,B>0be 
non-negative operators such that ker A n ker B = {0}, d m i n (A| supP A) > p and 
dmin(5| supp g) > p. Then 

A + B > 2p sin 2 ^ j (110) 

where 

cos 6 := max \{ik\(p)\. (Ill) 

\ij/)ekerA 
|^)eker B 


Proof (of Clairvoyance Lemma 43) Case (1) is trivial since {x\ c {i]/\ Q H ptn \x) c \i[/) Q > 

I for any illegal standard basis state |x) c . 

Now consider cases (2) and (3). By assumption, all legal standard basis states 
of He are contained in AV, which is closed under transition rules. Thus closure 
of S (Lemma 42) implies S c AV. Consider the directed graph of states in AV 
formed by adding a directed edge between pairs of states connected by transition 
rules. By assumption, only one transition rule applies in each direction to any state 
in AV, so the graph consists of a union of disjoint paths (which could be loops in 
case (2)). Minimality of S (Lemma 42) implies that S consists of a single such 
connected path. 

Let t - 0,..., |S | - 1 denote the states in S enumerated in the order induced by 
the directed graph. H trans then acts on the subspace 70, <g> Hq as 

T l 

H ^A<K s m Q = Z 9 ® 1 + I* + 1){t + 11 ® 1 “ lt + 1)<?l ® Ut ~ m + H ® U?) ’ 

52 t=0 z 

(112) 

where U, is the unitary on Hq appearing in the transition rule that takes \t) c to 
\t + l) c . T = \S | - 1 if the path in 5 is a loop, otherwise T = \S \ - 2. 

Consider a single term H t := |r)(r| ® 1 + \t + l)(t + 1| ® 1 + \t + l)(t| ® U t + 

| t)(t + 1| ® Uf from (112). Defining the unitary W t = \t + l)(t + 1| ® uj + (1 - 

I I + 1 )(t + 1|) ® 1, we have 

W t H t wJ = (|r> -\t+ 1» «f| - <t + 1|) ® 1 > 0. (113) 
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Thus, whether the path in S forms a loop or not, we have 1 


H{n A'K s &H Q 


|S |-2 . 

> ® i + \t + i)<t + 1 | <g> t - \t + i)<t| <s> u t 

t =o 1 

— • ^path? 


( 114 ) 


with equality if the path is not a loop. 
Defining the unitary 


|S |-2 t 

w=Yj io<^i ® n u i + 1(|51 ~ i))<(|51 ~ 1)1 ® 

t -0 i -0 


we have H pat h - WH p . dl u W f = E ® 1 where 


2 

1 

2 


0 


E = 





1 

2 


1 

2 


1 


0 0 ' 
0 


0 

0 0 


0 


0 


1 

1 

‘2 


1 

2 

1 

2 ' 


(115) 


(116) 


The matrix E is the Laplacian of the random walk on a line, and it is well known 
that its eigenvalues are given by A k = 1 - cos q k where q k = kn / |S|,k = 0,...,|S|-l, 
with corresponding eigenvectors \<p k ) oc cos (q k (j + |)) |/> [ CSV02], Thus 


[ 1 151-1 ) ( 1 151-1 ) 
ker Hp ath = span | W 1 -^= ^ |f> c |^ 0 ) G j = span | -^= ^ \t) c \i// t ) Q J (117) 

where | 0 v> := U t ...U 0 \ifro) for any |^ 0 >- 

First consider case (2). Take A = 7/ path and B = H ptn \-K s <d'Ho in Lemma 44. We 
have 4 min (A| SU p pA ) = 1 - cos 41 = D(1/|S| 2 ) and T min (5| suppB ) = 1. Furthermore, 
since S contains at least one illegal state in case (2), we have ker A n ker B = {0} 
and 

cos 2 e = max {<p Q \ c (if/\ Q W + n kerB W \<p Q ) c \if/) Q < 1 - 1 /\S |, (118) 

llW 

lr This operator has the same form as Kitaev’s Hamiltonian, and the analysis from here on is 
similar to that in Kitaev, Shen, and Vyalyi [KSV02]. 
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where n kerB is the projector onto ker B and we have used the fact that B = 
UptnWs&Hq is diagonal in the standard basis. Invoking Lemma 44, we obtain 

(tftrans + //pen)|^ > //path + = ^(1/|S| 3 ), (119) 

thus MXs) = &(1/IS | 3 ) as claimed. 

Finally, consider case (3). Since S contains only legal states in this case, 
H pen \% s&Hq = 0- Furthermore, by condition (iv) the states in S cannot form a loop, 
so H[x S 0 'h q = Hu--dns\'K s »/H Q = //path - E <g> 1. The claim in the Lemma now follows 
from the form of the eigenvalues and kernel of E, given above. □ 

Note that the proof of the Clairvoyance Lemma relied crucially on unitarity of 
the operators U a bcd appearing in the transition rules (see Lemma 42). In particular, 
it is not sufficient for the U a b C d to be partial isometries. That would not allow us to 
unitarily transform H muK \- Ks «/H Q into E <g> 1 using the unitary W in (115) and (116). 
Thus, in order to apply the Clairvoyance Lemma 43, we must complete the quantum 
parts of the transition rules to a full unitary. This necessitates the time-wasting 
construction in Section 7.6.2 (or similar). 

We are finally in a position to prove Theorem 33, the main result of this section. 
As commented above, let H c = (C C )® L and Hq = (C~)* L be the Track 1-3 
and Track 4-5 Hilbert spaces, respectively, for a chain of length L, where C c - 
span{| 5 , p, cr>} ® span{©, ©} and C e - spanflg, r) © span{@, ©}} (cf. (67)). Let 
/?trans e S((C C ® C e )® 2 ) be the sum of all the transition rule terms defined in 
Tables 1 and 3 (omitting those marked in Table 1) after completing to a full unitary 
so that the resulting Hamiltonian is standard-form. Let h pen e £>((C c ® C~)® 2 ) be 
the sum of all penalty terms defined in Table 2, and h init e ®((C c ® C e )® 2 ) be the 
sum of all Track 4-6 initialisation penalty terms defined in Table 4. 

Define the standard-form (Definition 41) Hamiltonian H(L) = i7 trans (L) + 
Hpen(L) e 'H c ® ‘Hq on a chain of length L, where /7 t rans,pen = Zf=i' and 

similarly define H nut (L) := Xf=Ti' /c'i' t +l) - 

Proposition 45 (Unique ground state) The unique 0-energy eigenstate of 
(Htmns(L) + H pen (E) + H init (L)) | is the computational history state 

1 T 

—j= Wt)c > (120) 

v/ 1 t = o 

where T = Q ^L|S x Q\ ! ^, \(J/o ) Q e 'Hq is the initial configuration of the QTM with 
the unary representation of the number L - K - 3 as input, and \ e Hq is the 
sequence of states produced by evolving \ifo) under the QTM and time-wasting 
counter TM transition rules (where each such state is duplicated 0(L ) times in 
succession in the sequence). 


Ill 


Proof Let us first apply the Clairvoyance Lemma to (H tians (L) + // pen (L))| . We 
have to check that h lnms and h ptn fulfil the requirements of the Lemma. 

Take 'He to be the Hilbert space of Tracks 1-3, and 'Uq the Hilbert space of 
Tracks 4-5. Define TV to be the set of all well-formed Track 1-3 standard basis 
states. Any state that is not well-formed violates an illegal pair enforcing a regular 
expression, so all legal Track 1-3 states are well-formed. TV therefore fulfils 
condition (i) of Lemma 43. By Lemma 39, the set of well-formed states is closed 
under the transition rules and at most one transition rule applies in each direction to 
any well-formed state, so TV fulfils conditions (ii) and (iii) of Lemma 43. Finally, 
by Lemma 40, the only subset of legal Track 1-3 standard basis states that is closed 
under the transition rules is the set of clock states {| (f> t )}, which has one state | <p 0 ) 
to which no backward transition applies, and one state \(p T } to which no forward 
transition applies. So TV also fulfils condition (iv) of Lemma 43. 

All penalty terms are diagonal in the standard basis, as required by Lemmas 42 
and 43, and all transition rules are of the form required in Lemma 43, by construc¬ 
tion. h trans and h pen therefore fulfil all the requirements of Lemmas 42 and 43. 

Invoking the Clairvoyance Lemma 43, the 0-energy eigenspace of (H aans + 
# P en)ls br has the form 



H, 


trans 



= span 



( 121 ) 


where \ifr t ) Q = U t ...U\ \^o)q and \ipo)o is any state of 'Hq. 

Now we include the h init terms, and show that \i//o) Q must be as claimed in 
Proposition 45. For 0 < t < 2L - 1, the clock states \<p t ) consist of an @ 
that sweeps from left to right and back along the chain (where for the time steps 
t = L+ 1,..., L+K of this sweep, the arrow is in the states @,..., (§), respectively). 
The transition rules that apply during the initialisation sweep act trivially on 
Tracks 4-5, so |= |t (/q)q for t < 2L - 1. By construction, the h inil penalty terms 
from Table 4 give an additional energy penalty to any Track 4-5 states that are not 
in the desired initial QTM configuration when the @ sweeps past. 

Consider any Track 4-5 state |0o) e that is not the desired initial QTM configu¬ 
ration. Then there is some 0 < t < 2L - 1 for which {(f>,\ c H init \4>,) c > 0. 

Noting that the overall Hamiltonian H > 0 and the \<p t ) are mutually orthogonal, 
this immediately implies that the unique 0-energy state is that given in (120). 

The initial state | \J/q) q on Track 4 matches the regular expression ©O *q©*©. 
The transition rules U t preserve this form - which allows us to identify Tracks 4 
and 5 with the head and tape configuration of the QTM - and by construction 
implement one step of the QTM or time-wasting counter TM in each © sweep (all 
other unitaries applied to Tracks 4 and 5 during each sweep being identities). 
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Finally, the clock counter TM advances through Q(| Q x L\ l ) configurations 
before it runs out of tape, and each step of this TM requires one complete steep of 
the Track 1 arrow which takes Q(L) steps. Thus T = 0.(L\Q x L\ L ) as claimed. We 
are done. □ 

The local Hilbert space dimension in (67) manifestly depends only on the 
alphabet size and internal states of the counter TM encoded on Tracks 2 and 3, the 
QTM encoded on Tracks 4 and 5, and the time-wasting counter TM encoded on 
Tracks 4 and 6. However, the alphabet size and internal state space of both counter 
TMs depends only on the base £ we choose, and this is completely determined by 
the alphabet size and internal state space of the QTM. Thus the local Hilbert space 
dimension d depends only on the alphabet size and number of internal states of the 
QTM, as required by part (i) of Theorem 33. 

Since the local interaction h of our Hamiltonian is constructed by summing 
transition rule terms and penalty terms, and h tIans , fi pen , h- ina > 0, we have h > 0. But 
we have shown that the overall Hamiltonian has a 0-energy eigenstate, hence it is 
frustration-free, as required in part (ii) of Theorem 33. 

Proposition 45 implies that the unique ground state of H\ SbI has the computa¬ 
tional history state form claimed in part (iii) of Theorem 33. Moreover, since the 
history state encodes £l(\Q x L\ l ) complete sweeps of the Track 1 arrow, and the 
QTM is advanced by exactly one step in each right-to-left sweep, the computational 
history state encodes Q(|<2 x L \ L ) steps of the computation. If the QTM halts before 
this, it transitions by the qj rule of (98) to running the time-wasting counter TM, 
which never alters the QTM tape encoded on Track 5. Track 5 is therefore only in 
the state qj for one time step. Thus part (iv) of Theorem 33 is also satisfied. 

Now imagine that the QTM reaches a configuration in which the next step would 
move the head before the starting cell. (Though in the case of the proper machines 
considered in parts (iv) to (vi) of Theorem 33, the former will never occur.) 
Consider the corresponding | iff t ) in the computational history state. The QTM 
transition rules do not alter Track 6, so Track 6 is still in its initial configuration. 
Since the QTM has deterministic head movement, then at no point during the 
evolution is its head ever in a superposition of locations. \\J/ t ) is therefore of the 
form: 


l^) = l©> 


Z ■ ■ ■ @) Track4 ® 


# 'Track 5 


K#...#y 


Track 6 


l©>, (122) 


where \ip q ) are unnormalised Track 5 QTM tape states, and all configurations in the 
superposition would step the QTM head left in the next time-step. The qf transition 
rule from (98) therefore applies to all states in the superposition, and the next state 
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must have the form: 


Wt+\) - I©) 


I Pa© • • • )xrack 4 ® Vrack 5 ® I 1- ? # • • • ^ 


Track 6 


l@). (123) 


/ 


Thenceforth, only time-wasting counter TM transition rules apply. 

A very similar argument applies if the QTM reaches a configuration in which 
the head moves beyond cell L. Since the time-wasting counter TM transitions 
never alter Track 5, part (v) of Theorem 33 follows. 

This concludes the proof of Theorem 33, and this section. 
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8 Quasi-periodic tilings 

In order to ultimately prove the required spectral properties of our final Hamiltonian, 
we require significantly stronger properties of quasi-periodic tilings than those 
used to prove tiling results. In particular, we will need strong “rigidity” results, 
which show that the quasi-periodic structure is in a sense robust to errors. As far as 
we are aware, these rigidity results are new. All the tiling results used in the proof 
of our main results are gathered in this section. 

We exploit the very particular properties of a quasi-periodic tiling due to 
Robinson, which we briefly review in the following section. The Robinson tiling 
has a hierarchical geometric structure that lends itself to inductive proofs of the 
required rigidity results. 

8.1 Robinson’s tiling 

We now describe Robinson’s quasi-periodic tiling. It was discovered by Robinson 
in 1971 [Rob71] as a tool to simplify Berger’s proof of the undecidability of the 
tiling problem [Ber66], 


(a) (b) (c) (d) (e) 

Figure 9: The five basic tiles of Robinson’s tiling. 

Robinson’s tiling is based on the five basic tiles showed in Figure 9 and all 
those that can be obtained from them by rotation or reflection. Tile (a) is called 
a cross, whereas tiles (b)-(e) are called arms. In a valid tiling, arrow heads must 
meet arrow tails. As drawn, cross (a) is said to face up/right. An arm is said to 
point in the direction of its only complete central arrow. 


(a) (b) (c) (d) 

Figure 10: Parity markings to force the position of the crosses. 

For each one of the five basic tiles of Figure 9, we introduce two possible 
colours, red or green, to each of the side arrows (note that there is one tile without 
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colour) with the following restrictions: at most one colour may be used horizontally 
and at most one colour may be used vertically; for the cross, the same colour must 
be used in both directions; for tile (b) one colour must be used horizontally and the 
other colour vertically. Moreover, we impose the restriction that green crosses must 
appear in alternate positions in alternate rows (and maybe in other positions too). 
This can be achieved by adding parity markings to the above tiles: the additional 
arrows showed in Figure 10, which should match properly in any valid tiling, with 
the following rules: parity marking (a) is associated with green crosses, parity 
marking (b) with horizontal green arrows, and (c) with vertical green arrows. Parity 
marking (d) is associated with all tiles. It is important to remark here that parity 
markings “live in a different layer” than the five basic tiles. That is, arrows from 
the parity markings must match only those from the parity markings and arrows 
from the basic tiles must match only those from the basic tiles. 


□□□□□□□□□□□□□□DO 



□□□□□□□□□□□□□□□□ 


Figure 11: A possible Robinson tiling of the plane. 

If we now draw only the red collared lines, one possible tiling of the plane looks 
like Figure 11 , which has the crucial (for our purposes) quasi-periodic structure 
consisting of squares of increasing size. 

Let us try to explain briefly (full details con be found in [Rob7 ]) how Figure 11 
emerges, and the remaining freedom it allows in constructing valid tilings of the 
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Figure 12: Left : A 3 x 3 square of the Robinson tiling, with the choice of a right/up 
facing cross in the middle tile. Right: A 7 x 7 square of the Robinson tiling, with 
the choice of a right/down facing cross in the middle tile. (To avoid confusion, we 
have only drawn the coloured lines of the tiles in this second figure.) 


plane. Since crosses must appear in alternate rows in alternate positions, together 
with the form of the five basic tiles, if we take any cross, this completely determines 
the structure of the 3x3 square constructed in the direction it faces. For example, 
assume we start with a green cross facing left/down. Then we are forced to complete 
the 3x3 square having it as the top/right comer as shown in Figure 12, where in 
the middle tile we must have a red cross, with the only freedom of choosing where 
it faces. 

Once this is fixed, the 7x7 square obtained by growing the 3x3 square in 
the direction this red crosses faces is fixed as shown in Figure 12, where there is 
a green cross in the central position, with the only freedom of choosing where it 
faces. Keeping this procedure gives a tiling like the one shown in Figure 11. It 
consists of a quasi-periodic structure of squares (called borders) of sizes 4" (which 
means following Robinson’s notation that the distance between two of the facing 
crosses delimiting the border is 4” + 1) repeated with period 2" +1 , for all n e N. 
Such a tiling can tile the whole plane, a half-plane or a quarter-plane, depending on 
how we sequentially choose the orientation of the central crosses. Lines between 
valid half or quarter planes must consist only of arms. If tiling divides the plane 
into two half-planes, the patterns in these half-planes may be shifted relative to 
one another by an arbitrary even number of cells. In this case, we call the line 
separating the two half-planes a fault. However, if a half-plane is further divided 
into two quarter-planes, no further shift is possible between these quarter-planes. 
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8.2 Rigidity of the Robinson tiling 

We have just seen that with the original Robinson tiles, the quasi-periodic pattern 
of borders does not necessarily extend throughout the entire plane; there can be a 
fault line between two half-planes, with “slippage” between the patterns on either 
side of the fault (see Section 8.1 and [Rob71, § 8]). It will be convenient to slightly 
modify the Robinson tiles in order to prevent fault lines. To this end, wherever one 
of the five basic tiles does not have side-arrows, we add dashed side-arrows parallel 
to the solid central arrows (see Figure 13). We then extend the five basic tiles to 
the full tile set by rotation and reflection, adding parity markings, and adding red 
and green colourings to the solid side arrows, exactly as in Section 8.1. We refer to 
these as the “modified Robinson tiles”. (Similar modifications to the Robinson tile 
set appear elsewhere in the literature, e.g. in [Mi§97].) 



(a) (b) (c) (d) (e) 

Figure 13: The five basic tiles of the modified Robinson tiling. 

Since the all original tile markings are still present, any tiling using the modified 
tiles must also be a valid tiling using the original tiles if the dashed side-arrows 
are ignored. It is easy to verify that, for any tiling in which the 2”-borders repeat 
periodically throughout the plane for all n (i.e. any tiling that does not contain a 
fault line), dashed side-arrows can be added to make it into a valid tiling using the 
modified tiles. 

However, if the tiling contains a fault line, there is no way to consistently add 
dashed side-arrows. To see this, assume without loss of generality that the fault 
line is vertical. Then there must exist two back-to-back crosses on either side of 
the fault that are not aligned, so that one points up and the other down. Consider 
the cross to the left of the fault line. The dashed side-arrow on this cross that points 
to the right forces a dashed arm to extend in that direction, towards the other cross. 
Similarly, the left-pointing dashed side-arrow on the cross to the right of the fault 
forces a dashed arm to extend to the left. These dashed arms must meet somewhere 
between the two crosses. But, since one cross points up and the other down, one 
of these arms has the dashed side-arrow on top, the other has it on the bottom, so 
there is no way to consistently join up the arms. 

The crucial property we will use is that the Robinson tiling has a very “rigid” 
structure. Even if we start from a tile configuration that contains a defect (a pair of 
non-matching adjacent tiles), outside of a ball around the defect, the defect has no 
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effect on the pattern of borders that appears in the rest of the tiling. The following 
Lemmas make this rigorous. 1 

Lemma 46 (Cross rigidity) In any tiling of a connected region with Robinson 
tiles, crosses must occur in alternate columns and in alternate rows. 

Proof Immediate from the observation that, for each “parity tile” in Figure 10, its 
neighbour in each direction is uniquely defined. □ 

Lemma 47 (Robinson rigifity) Consider a connected region of a 2D square grid 
made up of square blocks of size 2" +1 x 2" +1 . Any tiling of such a region with 
modified Robinson tiles must contain a periodic pattern of 2" -borders (some of 
which may be incomplete due to boundaries), repeating horizontally and vertically 
with period 2" +1 . That is, the tiling must contain the same periodic pattern of 
2 ” -borders as a section of the tiling of the infinite plane, up to translation and 
defects that do not affect the 2 ” borders. 

Proof By Lemma 46, crosses must occur in alternate rows and columns throughout 
the region being tiled. Following Robinson [Rob71], we refer to these crosses as 
“1-squares”. (Note that there may be additional crosses in the tiling that do not 
constitute 1-squares.) Pick an orientation for one of these 1-square crosses. Its 
solid side-arrows force a solid arm to extend in the directions it faces. Unless there 
is insufficient space to the boundary of the region, there must be a cross two cells 
away in the direction of each arm. Just as in Robinson [Rob71, §3], these solid 
arms force the orientation of any two facing crosses to be mirror images. 

Meanwhile, the dashed side-arrows force a dashed arm to extend in the di¬ 
rections that the cross faces away from. Again, these arms force the orientations 
of any back-to-back crosses to be mirror images. Thus, as soon as we pick the 
orientation of a single 1-square cross, the orientations of all adjacent crosses are 
forced, and hence the orientations of all 1-square crosses throughout the connected 
region. 

Base case Consider first the case n = 1 with 4x4 blocks. Since orientations of 
adjacent pairs of 1-square crosses are mirror images of each other, they are forced 
to form groups of four facing towards the centre of a 3 x 3 block. As in Robinson 
[Rob71, §3], these blocks must be completed to form “3-squares” with 2-borders 
running around the edges and crosses in the centre. (The orientation of the central 
crosses is not fixed.) In our modified Robinson tiles, the dashed side-arrows also 
force the orientations of back-to-back crosses to be mirror images, which forces 

'Our “rigidity” results for the Robinson tiling could certainly be strengthened, but they suffice 
for our purposes. 
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adjacent 3-squares to be aligned. Thus the pattern of complete 2-borders must 
repeat horizontally and vertically with period 4 throughout the region. 

However, in general some of the 1-squares will be too close to the boundaries 
of the region to form complete 3-squares. We need to show that these still force 
partial 2-borders. If a 1-square faces a boundary and is located at that boundary, 
then there is no further section of 2-border in that direction in any case, and we 
have nothing to show. If a 1-square faces a boundary and is more than two cells 
away from it, there must be another 1-square facing it before reaching the boundary, 
and the section of 2-border between them must be completed as usual. The only 
interesting case is if a 1-square faces a boundary one cell away from it. In this case, 
the solid side-arrows still force the cell between the cross and the boundary form 
an arm extending to the boundary. The only difference is that there is no facing 
cross in that direction to force the side-arrows of the arm to point inwards; they 
could instead point towards the boundary. Thus any of the arm tiles with solid 
side-arrows (in a suitable orientation) can be used to form this piece of the 2-border. 
For our purposes, we consider any of these to still constitute a section of 2-border. 
So partial 2-borders are still forced at the boundaries. 

This completes the proof for the base case n = 1. It will be important for later 
to note that, if a partial 2-border is large enough to contain the cell that would 
normally be at the centre of the border, then the central cross is still forced in that 
cell. To see this, note that because our region is made up of adjacent square blocks, 
there must be arms on at least two sides of the central cell (not necessarily opposite 
sides, if the partial 2-square is located in a corner of the region). Regardless of 
which tile is used to form these arms, the central arrows perpendicular to the arms 
must face outwards as usual, because the solid side-arrows on the arms must be 
on the side away from the central cell. But the only tile that has outward-facing 
central arrows on adjacent sides is the cross, so it is forced. 

n = 2 case Before proving the general case, it is instructive to consider the next 
case n = 2 (with 8x8 blocks). We can always divide each 8x8 block into four 4x4 
blocks and apply the previous argument, so the 2-borders must repeat horizontally 
and vertically with period 4, throughout the region, with a cross at the centre of 
each 2-border (whether partial or complete). Each 8x8 block must therefore 
contain at least one complete 2-border, and all 2-borders must be aligned (hence 
their central crosses are also aligned). 

Pick an orientation for one of these central crosses. The side-arrows of this cross 
force solid arms extending in the directions that the cross faces, and dashed arms 
extending in the directions it faces away from. Now, the only tile at which an arm 
can terminate is the cross. But there can be no crosses between adjacent 2-borders. 
Thus the arms must extend all the way from a central cross to its neighbours, again 
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forcing the orientation of both facing and back-to-back crosses to be mirror images 
of each other, so that they form groups of four facing towards the centre of a 
5x5 block. As in Robinson [Rob71, §3], these blocks must be completed to form 
“5-squares” with 4-borders running around the edges and crosses in the centre. As 
before, the dashed arms force adjacent 5-squares to be aligned. So the pattern of 
complete 4-borders repeats horizontally and vertically throughout the region, with 
period 8. 

It remains to consider crosses that are too close to the boundaries to form 
complete 4-borders. If the central cross of a 2-border faces a boundary (rather 
than another cross), its solid side-arrows still force an arm that extends towards 
the boundary. The only tile that can terminate an arm is a cross. However, we 
cannot place a cross anywhere within the 3-square surrounding the central cross, 
as it would break the tiling within the 3-square. The only remaining possibility 
is to place a cross in the one-cell-wide corridor that runs along each side of a 
3-square. However, where an arm intersects a 2-border it must point outwards, 
which prevents placing a cross in the corridor. Thus the arm must be extended all 
the way to the boundary to form a partial 4-border (with some additional freedom 
in which arm tiles are used to form this section of 4-border). 

This proves the Lemma for n = 2. As in [Rob71, § 3], a cross is forced at the 
centre of all complete 4-borders. If a partial 4-border is large enough to surround 
the cell that would be at the centre of the complete border, a cross is forced there 
too. To see this, observe that for partial borders of this size, the centre cell has 
partial 4-borders on at least two sides. The central arrows perpendicular to a 
border necessarily have arrow tails on the inside of the border. These force an arm 
extending back towards the central cell. That cell therefore has arrow tails on at 
least two sides, and the central cross is forced. 

General case The argument for general n is very similar to the argument for the 
n = 2 case, and proceeds by induction on n. Assume for induction that the pattern 
of 2”“ 1 -borders repeats horizontally and vertically with period 2” in any connected 
region made up of 2" x 2" blocks, and that a cross is forced at the centre of each of 
those 2 n ~'-borders (whether partial or complete). 

Now consider a connected region made up of 2 n+l x 2" +1 blocks. We can 
always divide each block into four 2" x 2" blocks. So, by assumption, the pattern 
of 2” -1 -borders must repeat horizontally and vertically with period 2" throughout 
the region, with a cross at the centre of each (partial or complete) 2”“ 1 -border. 
Each 2” +1 x 2” +1 block therefore contains at least one complete 2 n ~ 1 -border, and 
adjacent borders (hence also their central crosses) are aligned. The solid and 
dashed side-arrows of these central crosses force arms extending between adjacent 
crosses, which in turn force the orientations of adjacent crosses to be mirror images. 
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Thus the arms extending between facing crosses form 2"-borders, and adjacent 
2"-borders are aligned. The pattern of 2”-borders therefore repeats horizontally 
and vertically throughout the region, with period 2 " +1 . 

For central crosses that face a boundary instead of an adjacent cross, the solid 
side-arrows on the cross still force an arm extending towards the boundary. To 
show that this arm necessarily extends all the way to the boundary, we must 
show that it is impossible to place a cross anywhere along the arm. Within the 
(2 n - l)-square surrounding the central cross, we cannot place a cross along the 
arm without breaking the tiling within 2” - 1-square (cf. [Rob71, § 3]). Thus the 
arm must extend to the boundary of this (2” - l)-square. Similarly, we cannot place 
a cross along the arm anywhere within the adjacent (2" - l)-square. So, if the arm 
reaches the adjacent (2” - l)-square, it must continue onwards to the boundary. 
The only remaining possibility is to place a cross in the one-cell-wide corridor 
between the adjacent (2” - l)-squares. However, the arm necessarily points out of 
the (2” - l)-square, which prevents placing a cross in the corridor. 

Finally, for the induction step we must show that a cross is forced at the 
centre of each (complete or partial) 2"-border. The argument applies to any border 
that extends sufficiently far to surround its centre cell (which includes complete 
borders). The central cell has (complete or partial) 2”-borders on at least two 
sides. The central arrows perpendicular to a border necessarily have arrow tails 
pointing inwards, which force arms extending back towards the central cell. That 
cell therefore has arrow tails on at least two sides, forcing the central cross. 

Together with the base case n - 1, this completes the proof. □ 

We refer to the top edge of a 2”-border as a 2 n -segment} 

Lemma 48 (Segment bound) The minimum number of 2 n -segments in a tiling 
of an L x H rectangle (width L, height H) using modified Robinson tiles is 
lH/2 ,,+ '\({L/2 n+ '\ - 1), and this minimum can be attained simultaneously for 
all n. 

Proof Lemma 47 implies that the tiling must contain the usual periodic pattern of 
2"-borders, up to translation. Since the borders repeat vertically with period 2 n+l , 
the minimum number of rows of 2”-segments is [H/2 n+l ] (obtained by translating 
the pattern vertically so that there is a row of incomplete borders along the top edge 
of the region). Similarly, since the borders repeat horizontally with period 2" +l , the 
minimum number of 2"-segments in each row of segments is |_T/2" +l J - 1 (obtained 
by translating the pattern horizontally so that there are incomplete segments at the 
left and/or right edges of the region). 

*It is on these 2"-segments that the QTMs will ultimately “run” in the final construction. 
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Since any 2"-border in the Robinson tiling contains 2" _1 borders inside it, the 
translations can be chosen such that there are incomplete segments of all sizes at 
the boundaries of the region. The result is now trivial. □ 

The following Lemma is the key rigidity result that we will need later. It shows 
that a defect in the tiling (i.e. a non-matching pair of adjacent tiles) cannot affect 
the pattern of 2”-borders in the tiling outside a ball of size 2” +1 centred on the 
defect. 

Lemma 49 (Segment rigidity) In any tiling of an LxH rectangle (width L, height 
H) with d defects using modified Robinson tiles, the total number of 2" -segments is 
at least |H/2” +1 J(|L/2” +1 J - 1) - 2d. 

Proof Divide the LxH rectangle into contiguous square blocks of size 2” +1 x 2' i+1 , 
with |_L/2" +l J blocks in each row and |_77/2" +1 J in each column. (Any leftover 
cells play no further role in the argument.) Delete any block that contains at 
least one defect (arbitrarily assigning defects that occur at block boundaries to the 
bottom/left block). Consider any horizontal row of blocks. If the row contains dj 
deleted blocks, those missing blocks divide it into at most d, + 1 strips of height 
2«+i. 

Since each strip is a connected, defect-free region made up of 2 n+1 x2 ,!+l blocks, 
Lemma 47 implies that it must be tiled with the usual periodic pattern of 2"-borders, 
up to translations. The strip may be connected to strips in adjacent rows, in which 
case we are not free to translate the tiling in each strip independently. Nonetheless, 
we can lower-bound the number of 2"-segments by allowing the tiling within each 
strip to be translated independently, and minimising over all such translations. 

Each strip has height 2" +l , so, regardless of how the pattern is translated 
vertically, it must contain exactly one row of 2"-segments. The minimum number 
of segments within a strip of block-length / is / - 1 (obtained by translating the 
pattern horizontally so that there are incomplete segments at the beginning and 
end of the strip). In other words, each block in the strip effectively contributes one 
segment to the total, except for the block at the right end of the strip, which does 
not contribute (we assign the segment between adjacent blocks to the left one). 

If the row contains d, deleted blocks, those missing blocks divide it into at most 
di + 1 strips. 1 Thus the row contains |_L/2' i+l J - d t non-deleted blocks, with at most 
di + 1 of them located at ends of strips. The whole row therefore contains at least 
|_L/2" +l J - di - (di + 1) = |L/2 n+1 J - 2d, - 1 segments. 

Summing over all [H/ 2” +1 J rows, and noting that the total number of deleted 
blocks cannot be more than the total number of defects, we obtain the claimed 
lower bound on the total number of 2”-segments. □ 

1 For di > m/2, this substantially over-counts the number of strips, but it is sufficient for our 
purposes. 
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9 Putting it all together 

9.1 Undecidability of the g.s. energy density 

To prove our first main result, Theorem 6, we will prove a sequence of lemmas 
which allow us to combine together the Hamiltonian constructions from the previ¬ 
ous sections, and progressively build up the final Hamiltonian. 

We will repeatedly need to refer to ID translationally-invariant Hamiltonians 
with a particular set of properties. For conciseness, we will call these “Gottesman- 
Irani Hamiltonians”, captured in the following definition: 

Definition 50 (Gottesman-Irani Hamiltonian) Let C e be a finite-dimensional 
Hilbert space with two distinguished orthogonal states labelled |@), |@). A 
Gottesman-Irani Hamiltonian is a ID, translationally-invariant, nearest-neighbour 
Hamiltonian H q (r) on a chain of length r with local interaction h q e S(C~ ® C e ), 
which satisfies the following properties: 

(i) . h q > 0. 

(ii) . [M©X@I®I@X©I] = [fc„l©X@l®l©X©|] = [M©X©I®I©X©|] = 

[M©X©l®l©X©l] = o. 

(iii) . A()(r) := A( } (H q (r)\ Shr ) < 1, where Sb r is the subspace of states with fixed 

boundary conditions |©), |@) at the left and right ends of the chain, respec¬ 
tively. 

(iv) . Vne N : A 0 (4 n ) > 0, and £“ =1 A 0 ( 4") < 1/2. 

(v) . Ao(H q (r)\s < ) = A()(H q (r)\ s ) = 0, where S< and S > are the subspaces of 

states with, respectively, a | < ) at the left end of the chain or a |©) at the 
right end of the chain. 

Lemma 51 (Tiling + quantum layers) Let /iJ' 0W , lfi° l e !B( <D C ® C c ) be the local 
interactions of a 2D tiling Hamiltonian H c , with two distinguished states (tiles) 

| L ), | R) € C c . Let h q € S(C e ® C e ) be the local interaction of a Gottesman-Irani 
Hamiltonian H q (r), as in Definition 50. Then there is a Hamiltonian on a 2D 
square lattice with nearest-neighbour interactions /? row , h co1 e !B( C c+e+1 ® C c+e+1 ) 
with the following properties: For any region of the lattice, the restriction of the 
Hamiltonian to that region has an eigenbasis of the form \T) C ® \fi) q , where \ T) c is 
a product state representing a classical configuration of tiles. Furthermore, for 
any given \T) C , the lowest energy choice for \fi) q consists of ground states of H q (r) 
on segments between sites in which \T) C contains an \L) and an \R), a 0-energy 
eigenstate on segments between an \L) or |/?) and the boundary of the region, and 
|0) ’s everywhere else. 
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Figure 14: Intuition behind construction in Lemma 51: we “sandwich” together a 
Gottesman-Irani Hamiltonian and the Robinson tiling Hamiltonian, so that com¬ 
putational history states are forced to appear in the “top” layer above every every 
square border in the Robinson tiling of the “bottom” layer. 
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Proof The idea is to sandwich the two Hamiltonians H c and H q together in two 
“layers”, so that the overall Hamiltonian acts as H c on the c-layer, with constraints 
between the layers that force low-energy configurations of the c/-layer to be in the 
auxiliary |0> “blank” state, except between pairs of | L) and \R) states appearing in 
the same row of the c-layer, where the < 7 -layer acts like H q on that line segment. 
(See Figure 14.) 

To this end, define the local Hilbert space to be f H := *K C ® ( f H e © f H q ) - 
C c ®(|0)©C~). The Hamiltonian H is defined in terms of the two-body interactions 
as follows: 



(124a) 

(124b) 

(124c) 

(124d) 

(124e) 

(124f) 

(124g) 

(124h) 

(124i) 

(124j) 

(124k) 


=hT ® lg ® ^ 

+ 1® ® l< f+1) ® h q 

+ (1 C -|L)<L| C )®®|©)(©| ( °®1^ 1) 

+ 1% ® \R){R\t x) ® (l eq - l©X©l) (i+1) 

+ 1 ® ® (l c - \R)(R\ c ) {i+l) ® |©)<©| ( ' +1) 

+ 1® ® |0)<0|® ® \R)(R\ ( ‘ +1) © if 1} 

+ |L><L|® ® 1® ® l< i+1) ® |0>(0|f x) 

+ 1® ® |0)<0|® ® (l c - |L)(L| c ) (i ' +1) ® (leg ~ |0)<0| e ) (,+1) 

+ (l c - | R)<R\ C )® ® deg - |0)<0|g)® ® If 1} ® |0)<0|f 1} , 


where l c , t eq and t ceq are the identity operators on the corresponding Hilbert 
spaces. The Hamiltonian can be understood as follows. (124d) and (124e) force a 
|©) in the g-layer whenever there is an | L) in the c-layer. (124f) and (124g) do the 
same with |@) and \R). (124h) and (124i) force non-blank to the left and right of 
an | R) or | L), respectively. Finally, (124j) and (124k) force a non-blank to the left 
and right of any other non-blank in the <y-layer, except when a non-blank coincides 
with an |L> or \R) in the c-layer. 

Since h c is a tiling Hamiltonian, (124b) is by assumption diagonal in the 
canonical product basis on "Hf 2 (and acts trivially on ( ( H q © ‘K ( ,) 02 ). Meanwhile, 
(124c) is block-diagonal with respect to the four one-dimensional subspaces of < H® 1 
spanned by products of ©) and |@), together with the orthogonal complement 
thereof. ((124c) acts trivially on ‘Hf 1 .) The remaining terms are manifestly block- 
diagonal with respect to both of these decompositions simultaneously. The overall 
Hamiltonian is therefore block-diagonal w.r.t. the product basis on the c-layer. 
Hence there is a basis of eigenstates of H of the form \T) C \ip) q , where | T) c is a 
product state in the canonical basis of the c-layer. This proves the first claim of the 
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Lemma. 


For a given classical tile configuration | T) c on the c-layer, let X denote the set 
of all horizontal line segments £ that lie between an | L) and an | R) (inclusive) in the 
classical configuration \T\. Let Xl denote the set of all horizontal line segments 
between an L and the right boundary of the region, and similarly Xa> the segments 
between the left boundary and an R. 

Consider first a state \ijJo) q consisting of the ground state of H q (£) in the c/-layer 
for each £ e X, a 0 -energy eigenstate of H q (£) in the c/-layer for each £ e£ s U £ R , 
and |0) everywhere else in the c/-layer. The associated energy (L| c (i(/o\ q \H\ \T) C \4>o) q 
is clearly (T\H C \T) + Zfex Indeed (T\H C \T) is the contribution of (124a) 

and (124b), and Xcex / to(KI) the contribution of (124c), the rest of the terms being 0 
as | T) c \4>o) q satisfies all constraints imposed by (124d) to (124k). 

To see that this is the lowest energy state for a given classical configuration on 
the c-layer, we define a signature cr for each state in the canonical basis of the q- 
layer. The local Hilbert space at site i in the g-layer is { H e ® t H q , so we assign cr, = 0 
for states in f H e and cr, = 1 for states in 9-( q . By collecting computational basis 
states with the same signature, we decompose *H q into a direct sum of subspaces 
with given signature, 
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The overall Hamiltonian H is block-diagonal with respect to this decomposition, so 
all eigenstates | T) c \ij/) q have a \if/) q part with well-defined signature. We distinguish 
two cases: 

Case 1: cr-, = 1 for all i e £ e X. Consider a g-layer state \4>, T ) q with this 
signature; that is, a state supported on the corresponding Hilbert space H^. Since 
all Hamiltonian terms (124d) to (124k) are > 0, and for all £ e X the (124c) 
contribution is > X 0 (|X|), we easily get that (T\ c (ip ir \ q H \T) C \^ (r ) q > (T\H C \T) + 
Zcex^od^l), with equality iff \if/ a ) q = \if/ 0 ) q . 

Case 2: There exists i e £ e X such that cr, = 0. Let j be the rightmost position 
in £ such that cr ; = 0. If this is the right or left end of £, then the state picks up an 
energy contribution of 1 from (124d) or (124f). If it is the position next to the right 
or left end, it picks up a contribution of 1 from from (124h) or (124i). Finally, if it is 
not one of the above cases, it acquires a contribution of 1 from (124k). In all cases, 
the contribution is > Aq{\£\) since To(l^l) < 1 by assumption for a Gottesman-Irani 
Hamiltonian (see Definition 50). As all terms in the Hamiltonian apart from (124a) 
and (124b) are positive-semidefinite, and the contribution from (124a) and (124b) 
is (T\H C \T), we have that <T| C (^\ q H \T) C \ifr <T ) q > < T\H C \T) + Xcex ^o(l^l), which 
completes the proof of the Lemma. □ 


127 




With this, we can now prove the following: 


Lemma 52 (Robinson + Gottesman-Irani Hamiltonian) 

Let h r c ov, ,h c c ° 1 e S(G C <g> <D C ) be the local interactions of the tiling Hamiltonian 
associated with the modified Robinson tiles. For a given ground state configuration 
(tiling) ofH c , let £. denote the set of all horizontal line segments of the lattice that 
lie between down/right-facing and down/left-facing red crosses (inclusive). Let 
h q € £>(C~ ® C e ) be the local interaction of a Gottesman-Irani Hamiltonian H q (r), 
as in Definition 50. 

Then there is a Hamiltonian on a 2D square lattice of width L and height H 
with nearest-neighbour interactions h row , h co] e !B( C c+<2+1 ® C c+e+1 ) such that, for 
any L, H, the ground state energy 


A () (H ML >) = 


min 

XcA(L) 


X^d, 

£e£ 
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where the minimisation is over all valid tilings of the Lx H rectangle. 

Proof Construct the Hamiltonian H as in Lemma 51, with the red down/right- and 
down/left-facing crosses from the modified Robinson tile set as the designated |L) 
and | R) states. 

The tiling and QTM Hamiltonians satisfy the requirements of Lemma 51 , 
implying that the lowest energy for a given c-layer configuration is attained by 
putting |0)’s in the g-layer everywhere except in the segments between an |L) and 
an | R), inclusive. In the modified Robinson tiling, these are exactly the 4"-segments. 
Therefore, to bound the ground state energy, we can restrict our attention to classical 
configurations of Robinson tiles (not necessarily valid tilings) with an eigenstate 
of h q along each 4"-segment. 

By Lemma 48, the minimum energy of an eigenstate with a defect-free Robin¬ 
son tiling on the L x H rectangle is 


Llog 4 (L/2)J 

L (0 defects) = X ^o(l^l) = X 

£e£ n =1 

where the final sum is over all segment lengths 4". (Recall that we only use the red 
4”-segments to define line-segments on which h q acts non-trivially.) 

On the other hand, since each defect in the classical tile configuration con¬ 
tributes energy at least 1 from the h c term, Lemma 49 implies that the energy of an 
eigenstate with d defects on the L x H rectangle is at least 



Llog 4 (L/2)j, rr j \ 

L(u! defects) = J+XWI)>^+X ( ^TT “ 1 )“ 2 A 

£e£ n -1 ' ' 


Mo(4"). (128) 
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Since f] r A 0 (r) < 1/2 by assumption for a Gottesman-Irani Hamiltonian (see 
Definition 50), we have for all d > 0 


Llog 4 (i/2)J 

E(d defects) - £(0 defects) > d (1 - 2 ^ d 0 (4”)) > 0. (129) 

n -1 

Therefore, the lowest energy eigenstate always contains a defect-free tiling, and 
the Lemma follows. □ 

We can now apply this Lemma to construct a Hamiltonian h u with ground state 
energy that is undecidable even with a constant promise on the energy gap. 

Proposition 53 (Diverging g.s. energy) 

There exists a family of interactions h’f w (n), h c f l (n) e £((C U <E> (C u ) and h£\ri) e 
!B(<C U ) with operator norm < \ and algebraic matrix entries, and strictly pos¬ 
itive functions a[(n),ao(n),d 2 (n),a u l (n),di(n) (where the a functions are com¬ 
putable and the 6 functions are uncomputable), such that either AfHf'fn)) = 
-La\ (n) + afn), or AfH^in)) = Lfdoin) - L[n“(n) + 8\{ri)\ but determining which 
is undecidable. 

Moreover, the interactions can be taken to have the following form: h„\n) = 
a 2 (n)t with a 2 ( 11 ) an algebraic computable number, h r f" (n) {0,/3}-valued and 
independent of n and 

h c u °\n) =fi(A + e i7up B + e m2 ~ h 'c) + h.c.. 

where A e !B(<C U ® C u ) is independent of n and has coefficients in 7L + 

B,C e £>(<C U ® <C U ) are independent ofn and have coefficients in 7L, and/3 e Q 
is independent ofn. Recall that </? is defined as the rational number whose binary 
decimal expansion contains the digits ofn in reverse order after the decimal. 

Proof Let h q o be the Hamiltonian obtained by applying Theorem 33 with K = 3 
to the QTM from Theorem 11 with a proper reversible universal TM dovetailed 
after it. The Hamiltonian h q (n) in Lemma 52 will then be h q (n) = h q0 (n) + |T)(T| ® 
1 + 1 <S> \q f){c/ f\, where |T) is the halting state of the universal TM. h q has clearly 
the form given in part (vi) of Theorem 33. We claim that this Hamiltonian is a 
Gottesman-Irani Hamiltonian according to Definition 50. 

The only requirements in Definition 50 that are not immediate are those con¬ 
cerning the minimum eigenvalues restricted to various subspaces. Recall the 
construction of the computational history state Hamiltonian from Theorem 33. In 
the initialisation sweep, the @ normally sweeps once from left-to-right along the 
chain, turns around at the © to become a (f), which sweeps once right-to-left 
back along the chain. The superpositions over the standard basis states in these 
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sequences contain no illegal pairs as long as the other tracks are initialised correctly. 
However, if the © is missing, then when the @ reaches the right end of the chain, 
there is no forward transition out of the resulting state. Thus the uniform superpo¬ 
sition over the first part of the initialisation sweep, involving just the left-to-right 
sweep, is an eigenstate of h q . Furthermore, it has 0 energy if the rest of the tracks 
are correctly initialised. Similarly, if the © is missing, there is no further forward 
transition once the © reaches the left end of the chain, and the superposition over 
the full initialisation sweep is a 0-energy eigenstate of h q . Thus condition (v) of 
Definition 50 is satisfied. 

Let |0) = -jj: Yj T t =\ \<Pt ) I0r) be the normalised computational history state for 
the QTM, where T = Q(|X x Q \ r ) and | 0 ,) is the state encoding the f’th step of the 
computation. Note that 1 0 ) is a zero-energy eigenstate of H q0 , and at most one | 0 f ) 
can have support on the state |T) that represents the halting state of the universal 
TM, by Theorem 33. For r > 2, we have 


A 0 (r) < <0|/L(r)|0> = <01 


Y h T l) w + ® l <+i + 1 ‘ ® i t >< t i !+ i 

i 


10 ) 


T l ( 

Y Tj 0 1 '+1 + 1 «- ® l T >< T Ul 


t= 1 


I0r> < O 


1 


I SxglT 
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thus YjZ=\ do(4"') < 1/2. The remaining conditions (iii) and (iv) of Definition 50 
are therefore also satisfied. Hence h q (n ) is a Gottesman-Irani Hamiltonian, as 
claimed. 

Let ^ ow (n), h c u olumn (n) be the local interactions obtained by applying Lemma 52 
to h q (n). Let N(n) := max{p TOM, (n)||, ||/i CY,/ (n)||), and take a rational number^ < 
for all n. Such [3 exists by the form of h q guaranteed by part (vi) in Theorem 33 
and the definition of /?™ w (n), h c u olumn (n) based on h q . Define the normalised local 
interactions h T u ow (n) := / 3 / 0 ow (n), ^ olumn (n) := phf amn (n). 

For any r > \n\ + 6 , the QTM from Theorem 11 has sufficient tape and time to 
finish, and we can be sure that the reversible universal TM starts. If the universal 
TM does not halt on input n, then for all r > \n\ + 6 we have that A 0 (r) = 0. By 
Lemma 48, the minimum number of r-segments in any tiling of an L x L square 
(for r = 4 m , me N) is |_T/2rJ (LT/2rJ - 1) (and this minimum can be attained for 
all r simultaneously). Hence, as long as we take L > L (J (n) where L (J (n) is the 
minimal L such that the modified Robinson tiling of A(L) necessarily contains a 
4 m -segment of size 4 "' > \n\ + 6 , then Lemma 52 gives a ground state energy for 
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H Aa> (n) of 


min J>(M) =p Y l^l(ld- 1 )« r) <131a) 

le£ l<r<M+6 


l<r<|rc|+6 
r=4 m , me N 


= l- 


f I 


l<r<|«|+6 
r=4 m , me N 


jofr) 

4r 2 


- L 




l<r<|rc|+6 
r=4 m , raeN 


+ 


P Y ' 1 o(r )frac (d( frac (d + ‘) 


l<r<|rc|+6 
r-4 m , me N 


(131b) 


=: L 2 or 2 (« ) - La\(n, L) + ao(n, L ), 


(131c) 


where frac(x) := x - U'J denotes the fractional part of x. Note that the number of 
terms in the sums are finite, and for all finite r the quantity A(r) is an eigenvalue 
of a finite-dimensional matrix. Therefore, o^fn) is always a computable algebraic 
number. We also have 

a[(n):=(3 J] ^ < «i (n,L)<fi J] = : a »> (132) 

l<r<|«|+6 l<r<|rc|+6 

r=4 m , me N r=4 m , me N 

0<a 0 (n,L)</3 ^ =: a 0 (n). (133) 

l<r<|ft|+6 
r-4 m , raeN 

If the universal TM does halt on input n, then for any r larger than the size of 
tape needed to halt it is clear that A 0 (r) > 0. This follows immediately from the 
fact that the computational history state encoding the evolution (which necessarily 
has support on |T}) is the unique ground state of h q(h and h q o > 0. Let r\(ri) be the 
minimal such r of the form 4 m . Then by Lemma 52, the ground state energy of 
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H a<l > is 

A 0 (H ML) )={3 min Y^dfl) 

jcA(L) 




-p £ 

l<r<\n\+6 
r= 4 m , meN 


L \{ 

L 

_2rJ V 

.2r. 


L ! L 


r>r\(ri) 
r- 4 m , meN 


= L 


«2 (n)+p ^ 


Mr) 

4 r 2 

r>r\{n) 

r= 4 m , meN / 


(134a) 
jMr) 
(134b) 

Mr) (.„ (L 


“' < ") + ' , £ -^ 2 ( 2frac ( 5 7)+ 1 ) 


r>ri(«) 

r=4 m , meN 


+ 


Mn) +/3 ^ Mr) frac ( ^ j (frac ( ^ ) + 1) 


r>ri(«) 
r=4 m , meN 


(134c) 


=: L 2 [or 2 (n) + d 2 (ra)] - L[ari(n, L) + di(n, L)] + or 0 (n, L) + 6 0 (n, L). 

(134d) 


Note that b 2 (n) > 0, since in the halting case 4 0 (r) > 0 for all r > /-|(«), and 




34 0 (r) 


r>r\ ( n ) 
r=4 m , meN 


2 r 


r>ri(«) 

r=4 m , meN 


2r 


(135) 


0 < 8o(n, L) < f3 £ ^ ~:=S 0 (n ). (136) 

r>ri00 

r=4 m , meN 


We now modify /?(«) by adding the 1-body term = -a 2 (n)t acting at each 
site. The ground state energy of H AtL) is simply shifted by exactly -L 2 a 2 (n), so in 
the non-halting case we have 

A 0 (H A(L) ) = -La x (n, L ) + a 0 (n, L) < -La\(n) + a 0 (n). (137) 

In the halting case, we have 

A q (H A{L) ) = L 2 8 2 (n) - L[ ai (n, L ) + 6 An, L)\ + a 0 (n, L) + S 0 (n, L) 

n (138) 

> L 2 6 2 (n) - L[a\(n) + 8An)\. 

The Proposition follows from undecidability of the Halting Problem. □ 


The following corollary follows immediately by letting L(n) be the minimal L 
in Proposition 53 such that -La\{n) + ao(n) < 0 and L 2 5 2 {n) - L\a l \(n) + 6\(n)\ > 1. 
This will be useful shortly. 
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Corollary 54 (Undecidability of g.s. energy with promise) There exists a fam¬ 
ily of interactions h" m ’(n), h c u ol (n ) e £(<C U ®C U ) and h^\ri) e £>(< 0 ^) with operator 
norm < \ and algebraic matrix entries, and an (uncomputable) function L(n), such 
that either AfH^^in)) < 0 for all L > L{n), or XfH'^in)) > 1 for all L > L(n), 
but determining which is undecidable. Moreover, the interactions can be taken to 
have the same form as in Proposition 53. 

Undecidability of the ground state energy density is now immediate from 
Proposition 53 by the definition E p := lini/^oo Aq(H a<l> )/L 2 of the ground state 
energy density. We restate this result here for convenience: 

Theorem 6 (restated) Let d e N be sufficiently large but fixed, and consider 
translationally-invariant nearest-neighbour Hamiltonians on a 2D square lattice 
with open boundary conditions, local Hilbert space dimension d, algebraic matrix 
entries, and local interaction strengths bounded by 1. Then determining whether 
E p = 0 or E p > 0 is an undecidable problem. □ 

9.2 From g.s. energy density to spectral gap 

We are finally in a position to prove our main result: undecidability of the spectral 
gap, which we restate here for convenience. Recall that for each natural number 
n, we define if = <fi(n) to be the rational number whose binary decimal expansion 
contains the digits of n in reverse order after the decimal. 

Theorem 3 (restated) We construct explicitly a dimension d, d 2 x d 2 matrices 
A, B,C,D and a rational number (3 so that 

(i) . A is Hermitian and with coefficients in 7L + [3h + 

(ii) . B, C have integer coefficients, 

(iii) . D is Hermitian and with coefficients in {0,!,/?}. 

For each natural number n, define: 

hi(n) = a 2 (n)t (a^in) an algebraic number) 

^row(ft) = D (independent of n) 

h co \(ri) = A +f(e inv B + e^B' + e in2 ~ W C + . 


Then: 

(i). The local interaction strength is < 1. I.e. all terms h\(n), h„ m (n), h CQ \(n) have 
operator norm bounded by 1. 
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(ii) . If UTM halts on input n, then the associated family of Hamiltonians \ H ML> (n)\ 

is gapped in the strong sense of Definition 1 and, moreover, the gap y > 1. 

(iii) . If UTM does not halt on input n, then the associated family of Hamiltonians 

\H ML> (n)\ is gapless in the strong sense of Definition 2. □ 


Proof We assign a Hilbert space < H (,) := |0> ® 'H„ ® H( d - C © C 6 ' <g> C° to each 
site i e A, where |0) denotes the generating element of C. Let h„ be the two-body 
interaction obtained in Corollary 54 (keeping the one-body term h^in) separate, 
for later convenience). 

Let h d be the normalised (so that \\h d \\ < \) two-body interaction of any nearest- 
neighbour Hamiltonian H d with 0 ground state energy whose spectrum becomes 
dense in the thermodynamic limit (cf. Definition 2): lim L ^ txl spec H y iL) — > [0, oo). 
(For example [LSM61], we can take D = 2 and h d to be the critical XY-model 
h row = (T x ®cr x + (Ty®(Ty + CT Z ®1 + t®(T z along the rows, where (T xy ~ are the 
Pauli matrices, with no interactions along the columns.) 

Define the Hamiltonian H(ri) in terms of its two-body and one-body interactions 
h(n) as follows: 

h (Ui \n) := |0><0| ( ° ® nj/j (139a) 

+ (139b) 

+ 1®®1 (139c) 

h®(n) := - a 2 (n) n®. (139d) 


Clearly, both have norm bounded by 1. 

As in Lemma 51, t u , t d denote the identity operators on 'H„, 'H d , and Il ud 
denotes the projection of onto its 'H u ® < H d subspace. We decompose the global 
Hamiltonian in the square A(L) as H A(L} =: Hq + H u +H d -a^lud-, where Ho, H u , H d 
are defined by taking the sum over sites separately for the expressions in (139a) 
to (139c), respectively. Once again, we assign a signature cr € {0,1} L to every 
state in the computational basis, by assigning a 0 to the |0> state, and a 1 to any 
state in the ud subspace, and decompose the Hilbert space into subspaces of given 
signature, 




(140) 


iG A 


cr 


Since (139) is manifestly block-diagonal with respect to this direct sum decompo¬ 
sition, we have 


spec if A(L) = M specify 

cr 


(141) 


Therefore, it is sufficient to analyse the spectra for each signature separately. We 
distinguish three cases: 
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(i) . cr = (0,0): In this sector the Hamiltonian is 0. 

(ii) . cr = (1 ,..., 1): In this sector H 0 = 0 and H u + H d - a 2 (n) Zi = H„ (L) <g> 

1 + 1 ®H^ l \ where H‘^ L> = ZijeA hu’ j \n) - a 2 (n ) Zi 1« ' is the Hamiltonian 
from Corollary 54. Thus the spectrum stemming from this subspace is 
spec/^ (L) + specHj (L) . 

It only remains to analyse the case of mixed signatures: 

(iii) . cr ^ (0 ,..., 0) but cr, = 0 for some i. 

Let | iff) be any state with this signature. Since all cr, = 0 sites must be in the 
state |0), the overall state is of the form |^) ( , :cr . =1| (^) . 7T =0 |0) y -. We must show that 
any such state has higher energy than either the cr = (0 ,..., 0) or cr = (1,..., 1) 
state (whichever is the lower). 

Intuitively, 0-signature regions of the lattice reduce the number of QTMs 
“running” in the lattice in two ways: by reducing the total area that contains 
QTMs, and by breaking QTMs at boundaries between signature-1 and signature-0 
regions. The former reduces the energy contribution from the H u (139b) by 
something that scales as a 2 {n ) x area. But it also misses out on the negative 
-a 2 {n) x area contribution from the -a 2 (n)t ud (139d). So the net effect of the 
reduction in area occupied by QTMs cancels out. Meanwhile, breaking the QTMs 
at the boundaries reduces the H u energy contribution by something that scales 
as Zr 4 (j(r) x length, since each length of boundary can break one QTM of each 
size r. But the boundaries also pick up a positive 1 x length energy penalty from 
Ho ((139a) term). Since Zr^o(f) < 1, the energy penalty outweighs the energy 
reduction. Thus, overall, signature-0 sites within signature-1 region always increase 
the energy, and the ground state is never a mixed-signature eigenstate. 

We now make this rigorous, analysing the halting and non-halting cases sepa¬ 
rately. 1 

Non-halting case: Let M denote the total number of signature-0 (“marked”) 
sites, U the total number of signature-1 (“unmarked”) sites, and B the total number 
of boundaries between signature-0 and signature-1 sites. Now divide the L x L 
lattice into blocks of size 2r x 2r with [L/2r\ blocks in each row and column, 2 
where r = 4 m , m e N. We classify blocks according to the signature within the 
block: unmarked blocks, whose signature contains only l’s; completely marked 
blocks, whose signature contains only 0’s; and partially marked blocks, whose 
signature contains both 0’s and l’s. 

'The argument is reminiscent of Lemma 49. 

2 The remaining L - 2r[L/2r] columns will play essentially no role. 
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We first analyse the contribution of the H u + H d terms. Consider any horizontal 
row i of blocks, and let c, and p, denote, respectively, the number of completely and 
partially marked blocks in that row. The number of unmarked blocks in the row is 
then \_L/2r\ - c, - p t . Let /?, denote the number of unmarked blocks immediately to 
the left of a marked block (complete or partial). The marked blocks divide the row 
into strips of contiguous unmarked blocks. Within each such strip, the Hamiltonian 
acts as Hu" v + H^ np > H^, np . By Lemma 52, the lowest energy eigenstate of H^ np 
has energy min£ CA( x) YjC<e£ d () (|f|), where the minimisation is over all valid tilings 
of the strip. By Lemma 48, the minimum number of r-segments in a strip of 2r x 2r 
blocks of length / is / - 1; in other words, each unmarked block contributes one 
r-segment, except for blocks at the ends of strips. Thus the energy contribution 
from the r-segments in the row is at least 


( 


L 
.2r . 



(142) 


Now, each boundary between an unmarked and a marked block (complete or 
partial) contributes at least one additional vertical boundary between a signature-0 
and a signature-1 site. Meanwhile, each partially marked block contributes at least 
one additional 0/1 signature boundary. Thus, summing over all rows, we have 


J]pi + bi<B. (143) 


Also, each completely marked block must contain at least 4r 2 signature-0 sites, 
thus 


I 



(144) 


Summing (142) over all \_L/2r\ rows of blocks and over all segment lengths 
r, we obtain the following lower bound on the total contribution of the H u + H d 
terms: 


m h u + hm) > Yj Z {[t\ - Ci - Pi ~ bi ~ l ) io(r) 


r=4 m i 
me N 




Mr) 

4 r 2 


(145a) 

(145b) 


Now consider the energy contribution of the Hq - a 2 (n)t uq terms. Ho gives an 
energy contribution of 1 to every boundary between signature-0 and signature-1 
sites, whereas -a(n)t uq gives a negative energy contribution of - a 2 (n ) to every 
signature-1 site. Thus the total energy contribution from these terms is 

MH 0 - a 2 (n) l uq \i/s) = B - a 2 (n)(L 2 - M ). (146) 
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Putting (145) and (146) together, the total energy of any mixed-signature state 
is at least 


{m HL) w) = m ht + h* l) +# 0 A(L) - a 2 (n) m 

+ B- a 2 (n)(L 2 - M ) 

n / 

L 


(147a) 


Sis 


-”(? 


([^\- [ ) A o(r)-L 2 a 2 (n) +B 1 - £ A 0 (r) 
^ J V r 

Mr) 


4 r 2 


- a 2 (n) 


>A 0 (tf A(t) (n)) + | - M J] 


r>\n\+6 


Mr) 

4 r 2 


(147b) 


(147c) 


In the final inequality we have used the fact that a 2 (ri) := Yj r <\n\+6 Ao(r)/4r 2 and 
A 0 (r) < 1 /2, from Proposition 53 and Definition 50. 

In the non-halting case, Ao(r) = 0 for all r > |«| + 6. Furthermore, for any 
mixed-signature state, the total number of 0/1 signature boundaries B >2. Thus 
(147) implies the bound (if/\H K(L) \if/) > Ao(H^ iL> ) + 1. 


Halting case: The argument is similar to the non-halting case, but we count 
things slightly differently. Let U denote the total number of signature-1 (“un¬ 
marked”) sites, and B the total number of boundaries between signature-0 and 
signature-1 sites. We again divide the Lx L lattice into blocks of size 2 r x 2 r, 
with |_T/2rJ blocks in each row and column. Let r R r be the set of blocked rows 
that contain only signature-1 sites, with R r = \R r \ denoting the size of this set. As 
before, we categorise the blocks into unmarked, completely marked and partially 
marked, letting P r denote the total number of partially marked blocks. 

We first analyse the contribution of the H u + Hj terms. Consider any horizontal 
row i of blocks. Let n, denote the number of unmarked blocks in row i, and b, 
the number of unmarked blocks which do not have an unmarked block to their 
right (either because they are adjacent to a completely- or partially-marked block, 
or because they are the last block in the row). By the same argument as above, 
Lemmas 48 and 52 imply that the total number of r-segments in the i-th row must 
be at least n, - /?,. Thus the energy contribution from the r-segments in the row is 
at least 

(u t - bj)Ao(r). (148) 

Now we count the number of 0/1 signature boundaries. An unmarked block 
with a marked block on its right contributes at least one boundary between a 
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signature-1 site to the left of a signature-0 site. An unmarked block at the end of a 
row does not contribute this boundary. However, in this case, if the row contains 
at least one marked block, there must exist at least one unmarked block that has 
a marked block on its left, which contributes an additional boundary between a 
signature-1 site to the right of a signature-0 site. Therefore, in any row i that 
contains at least one marked block, there are at least b, 0/1 signature boundaries. 
Finally, in a completely unmarked row we have b, = 1 from the unmarked block at 
the end of the row. 

The number of vertical 0/1 signature boundaries is 0 if a row belongs to r R r and 
at least 1 if it does not. Thus rows in ( R r contribute one less 0/1 boundary than the 
others. Meanwhile, each of the P r partially marked blocks contributes at least one 
additional 0/1 boundary. Therefore, summing over all rows i, we have 

P r + ^ b, - R r < B. (149) 


Each unmarked block contains 4r 2 signature-1 sites, and each of the P, partially 
marked blocks contains at most 4r 2 - 1, out of a total of U. Thus 




U - 4 r 2 P r 
4 r 2 



(150) 


Finally, each row in % contains 2rL signature-1 sites, so 


R r < 


U 

2rL 


(151) 


Summing (148) over all rows i and over all segment lengths r, and using (149)- 
(151), we obtain the following lower bound on the total contribution of the H u + H d 
terms: 


+ H d \f) > 


> 


> 


2 - bi ) A <> (r ) 

(152a) 

r=4 m i 
me N 


r 

(152b) 


(152c) 


r 


Now consider the energy contribution of the H 0 - a 2 (n)t uq terms. As before, 
Ho gives an energy contribution of 1 to every boundary between signature -0 and 
signature -1 sites, whereas -a(n) 2 t uq gives a negative energy contribution of -a 2 (n) 
to every signature-1 site. So the total energy contribution from these terms is 

<tA|// 0 - a 2 {n)l uq \\J/) = B - a 2 (n)U. (153) 
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Putting (152) and (153) together, the total energy of any mixed-signature state 
is at least 


= <<AI H A{L) + H^ l) + H A(L> - a 2 (n)m (154a) 

£ Z( + (154b) 

r 

= u[Z ^- --(">) + b[i-- Z *«)^ 

V r / V r ) r 

>t/U(n)-i J + f. (154d) 


In the final inequality, we have used a 2 (n) := Zi< r <|„|+6^o( r )/4r 2 , 

S 2 (n) := 2r>i«i+6 doO')/4r 2 and do(r) < 1 /2, all from Proposition 53 and Defini¬ 
tion 50. 

For any mixed-signature state, the total number of 0/1 signature boundaries 
B > 2. Moreover, in the halting case, S 2 (n) > 0 independent of L. Thus, for 
sufficiently large lattice sizes L, (154) implies the bound > 1. 


Denoting the part of the spectrum which corresponds to these mixed signatures 
by 5 c R, we obtain that 

spec H ML) = {0} U j spec tf A(L) + spec H^ (L) ] U S. (155) 

Recall that the spectrum of H d becomes dense in the thermodynamic limit: 

lim spec H^ L> -» [0, oo). (156) 

L—>oo 

By taking L > L(n) and using Corollary 54, one of two things may happen. If the 
universal TM does not halt on input n, then S > + 1, A 0 (H^' (i) ) < 0, and 

spec H ML} becomes dense in [Ao(H'^ L) ), co) as L goes to infinity. If the universal 
TM does halt, then S > 1 and AoiH^ 11 ^) > 1, thus A(H a(l> ) > 1. Determining 
which is clearly undecidable, thus concluding the proof of our main Theorem. □ 

9.3 Periodic boundary conditions 

The previous section proves Theorem 3 for open boundary conditions. This is 
arguably the most important type of boundary conditions in the context of the 
spectral gap problem, both physically, because periodic boundary conditions rarely 
occur in real physical systems, and mathematically, because the thermodynamic 
limit is better behaved. Nonetheless, our result can also be extended to other types 
of boundary condition. Extending it to fixed boundary conditions is trivial, so 
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we omit the argument. In this section, we consider the more interesting case of 
periodic boundary conditions. 

For this, we must first address an issue that does not occur for open boundary 
conditions. For any local interaction h > 0 on a lattice with open boundary 
conditions, it is easy to see 1 that the ground state energy is non-decreasing with 
the lattice size: for any L 2 > L u A 0 (H(Li)) < A 0 (H(L 2 ). However, this is no 
longer necessarily the case for other types of boundary conditions. In particular, 
for periodic boundary conditions, there are trivial examples in which the ground 
state energy oscillates periodically with the lattice size. A simple example is the 
classical tiling Hamiltonian on a line, with N tile-types and the constraint that tile 
type n must be to the left of type n + 1 mod N. For periodic boundary conditions, 
this has ground state energy 0 for lines of length L = 0 mod N, but has ground 
state energy 1 otherwise. 

Using a very similar construction to that in the proof of Theorem 3 (see Sec¬ 
tion 9.2), this oscillatory behaviour of the ground state energy translates into 
oscillatory behaviour of the spectral gap. If we restrict the lattice size L to the 
subsequence with L mod N = 0, then the system is gapless. If we restrict to any 
subsequence with L mod N ^ 0, then the system is gapped with spectral gap 1. 
Our Definitions 1 and 2 of gapped and gapless systems are therefore not adequate 
for periodic boundary condition. 

Since our Definition 2 of gapless systems expresses convergence of the se¬ 
quence of spectra to a dense subset, the most natural generalisation appropriate to 
other boundary conditions is to relax this requirement to a convergent subsequence 
of lattice sizes: 

Definition 55 (Gapless - generalised) We say that a family [H A(L> \ of Hamiltoni¬ 
ans characterises a gapless system if there is a constant c > 0 such that for all 
s > 0 there is an Lq e I I and a subsequence {L,} of N such that for all L, > Lq 
every point in [Ao(H A<L ‘ } ), Af H MLl) ) + c] is within distance s from spec H A(L,) . 

With this generalised definition of gapless systems, Theorem 3 proving undecid¬ 
ability of the spectral gap generalises to the case of periodic boundary conditions. 
The proof involves slightly modifying the Robinson tiling part of our construction, 
using a tiling trick from [GI09]. We first restate and prove this trick here (in a 
slightly generalised form). 

The construction works for any (unweighted) tile set. First, duplicate all the tiles 
to create black- and white-coloured variants of each tile, and impose the constraint 
that adjacent tiles must have different colours. Turn this into an equivalent weighted 

'If If) is a ground state with energy E 0 = 'L<ij>eML- l )xML 2) Tr A(L2) (h il ' J} \d)(^\) for the lattice of 
size L 2 with open boundary conditions, then the mixed state p = Tr^ A (L 0 l0 r )(<AI on the lattice of size 

Li < L 2 has energy Tr A(Ll) (Z< ! j> eA (z. 1 )xA(i 1 ) h (iJ) p) = Tr A(Ll) (Z< 1 ;>€A(l 1 )xA(l 1 ) h (ui) I^X<Al) < E 0 . 
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tiling by defining the weight of adjacent pairs of matching tiles to be 0, and that 
of non-matching pairs to be +5. (Note that adjacent pairs can fail to match either 
because the original tile markings fail to match, or because their colours are the 
same.) Finally, add to the tile set the three new tile types shown in Figure 15 (which 
do not have additional colour markings). 


Figure 15: Boundary tiles for periodic boundary condition construction. 

The weights for adjacent pairs of tiles in this extended tile set are given in the 
following table, where □ and | denote any tiles of the given colour 1 (the weight 
for adjacent pairs of opposite-colour tiles is 0 if the original tile markings match, 
+5 if they do not match): 



□ 

Right tile 

■ B 

□ 

BB 



□ 

Tile above 

■ B 

□ 

ffl 

JD 

□ 

5 

0/5 

5 

0 

5 

£ 

o 

□ 

5 

0/5 

0 

5 

5 

• i—H 
+-> 

■4—> 

■ 

0/5 

5 

5 

0 

5 

13 

X) 

■ 

0/5 

5 

0 

5 

5 

CD 

B 

5 

5 

0 

5 

1 


B 

0 

0 

5 

5 

5 


□ 

0 

0 

5 

5 

5 


□ 

5 

5 

5 

0 

1 


ffl 

5 

5 

1 

5 

5 


B 

5 

5 

5 

1 

5 


(157) 


The following Lemma encapsulates the construction given in [GI09, § 4.2], 
It shows that these weights effectively turn the weighted tiling problem for an 
odd-sized square region with periodic boundary conditions (a torus), into the 
unweighted tiling problem for the original tile set on an even-sized square region 
with open boundary conditions. 

Lemma 56 (Periodic to open b.c.) Given any unweighted tile set and any odd 
L £ N, all minimum weight tilings of an Lx L torus using the extended tile set 
described above, with the weights given in (157), consist of a single m tile at an 
arbitrary location, a complete row of 3 t ^ es extending horizontally from the B 
a complete column of □ tiles extending vertically from the B and a valid tiling 
of the remaining (L - l)x(L- 1) region with the original tile set using alternating 
colours. 


'The choice of weights is somewhat arbitrary. Here, we choose them slightly more symmetri¬ 
cally than in [GI09], which makes the proof of Lemma 56 marginally simpler to describe. 
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Proof The tilings described in the Lemma have total weight +4 (coming from the 
□ and Q tiles adjacent to the ^). We must show that all other tilings have 
higher weight. Since L is odd, the toms cannot be tiled consistently with alternating 
colours, thus any weight < 5 tiling must contain at least one □.ra or ffl. 

If the tiling contains any tile other than R or ffl to the left or right of a R, 
it has weight > 5. Thus, in any weight < 5 tiling, R tiles can only appear in 
complete rows of 0 and ffl tiles. The analogous argument for columns implies 
that □ tiles can only appear in complete rows of □ and EB tiles. A EE adjacent 
to a black or white tile has weight +5, so EB tiles can only appear where such rows 
and columns meet. Moreover, adjacent R and □ tiles have weight +5, so all 
such rows and columns must meet at m tiles. 

Therefore, any tiling with weight < 5 must consist of some non-zero number 
of rows of Q til es and columns of □ tiles meeting at EB tiles, with coloured 
tiles everywhere else. Each intersection of a R row with a □ column at a EB 
contributes weight +4, so the weight is minimised by having just one such row and 
column. 

Finally, the remaining region to be filled with coloured tiles is a square of size 
(L-l)x(L- 1). L— lis even, so this can be tiled consistently with alternating 
colours. It contributes weight 0 iff, in addition, the tiles (ignoring the colour 
markings) form a valid tiling for the original tile set. □ 

Applying Lemma 56 to the modified Robinson tiles of Section 8, we obtain a 
weighted tile set whose minimum weight tilings on a torus of size L have exactly 
the same quasi-periodic structure as the original Robinson tilings on a square of 
size L— 1, for any odd lattice size L. Using this tile set in Lemma 51, the proofs of 
Lemmas 51 and 52, Proposition 53, and Corollary 54 go through unchanged for 
odd-sized square lattices with periodic boundary conditions. Undecidability of the 
ground state energy density (Theorem 6), and undecidability of the spectral gap 
(Theorem 3) under the appropriate Definition 55 of gapless, follow immediately. 
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